A bicubic polynomial approximation to Kepler's equation for elliptic orbits is shown to provide accurate starting values for efficient numerical solution of this equation for eccentric anomaly. In the approximate equation, a cubic in mean anomaly is set equal to a cubic in eccentric anomaly. The coefficients in the two cubics are obtained as functions of eccentricity by specifying values of function and slope at the midpoint and the endpoint of the complete interval (0 to pi). The initial estimate of eccentric anomaly to use in an iteration formula is obtained by evaluating the cubic in mean anomaly and finding the single real root of the cubic in eccentric anomaly. Numerical results are presented which indicate that the estimate accuracy of this method is roughly an order of magnitude better than that of other recently-reported formulas.
Efficient inverse solution of Kepler's equation
1986-12-01
Sonstige
Keine Angabe
Englisch
Solution of Kepler's equation as a differential equation
AIAA | 1984
|Bounds on the solution to a universal Kepler's equation
AIAA | 1978
|Efficient Solutions of Kepler’s Equation via Hybrid and Digital Approaches
Online Contents | 2015
|