In this paper we find the exact solution of the nonlinear differential equation describing the trajectory of mercury considering relativistic effects. Instead of the classical perturbation method we use the Jacobi elliptic functions to obtain the exact solution. The deviation of the perihelion of Mercury and other planets is calculated. The results may be of great interest ot astronomers and for those interested in the study of nonlinear differential equations.