Abstract
An efficient algorithm is presented for the solution of Kepler's equationf(E)=E−M−e sinE=0, wheree is the eccentricity,M the mean anomaly andE the eccentric anomaly. This algorithm is based on simple initial approximations that are cubics inM, and an iterative scheme that is a slight generalization of the Newton-Raphson method. Extensive testing of this algorithm has been performed on the UNIVAC 1108 computer. Solutions for 20 000 pairs of values ofe andM show that for single precision (∼10−8) 42.0% of the cases require one iteration, 57.8% two and 0.2% three. For double precision (∼10−18) one additional iteration is required. Single- and double-precision FORTRAN subroutines are available from the author.
Similar content being viewed by others
References
Brouwer, D. and Clemence, G.: 1961,Methods of Celestial Mechanics, Academic Press, N.Y.
Danby, J.: 1962,Fundamentals of Celestial Mechanics, McMillan, N.Y.
Smith, O. K.: 1961, ‘Terminating the Iterative Solution of Kepler's Equation’,Am. Rocket Soc. J. 31, 1598.
Traub, J.: 1964,Iterative Methods for the Solution of Equations, Prentice-Hall, N.J.
‘Handbook of Mathematical Functions’, National Bureau of Stds Appl. Math. Series, No. 55, p. 17, 1964.
Hart, J. et al.:Computer Approximations, J. Wiley, N.Y.
Ng, E.: 1971, ‘A General Algorithm for the Solution of Kepler's Equation for Elliptic Orbit’, JPL Tech. Mem. 33-496.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Ng, E.W. A general algorithm for the solution of Kepler's equation for elliptic orbits. Celestial Mechanics 20, 243–249 (1979). https://doi.org/10.1007/BF01371365
Received:
Accepted:
Issue Date:
DOI: https://doi.org/10.1007/BF01371365