Local finite Taylor series expansion is shown to be an efficient numerical method for the study of the n-body problem. Application is made to the accurate determination of the motion of the perihelion of Mercury according to Newtonian celestial mechanics.