Abstract:A fundamental relation in celestial mechanics is Kepler's equation, linking an orbit's mean anomaly to its eccentric anomaly and eccentricity. Being transcendental, the equation cannot be directly solved for eccentric anomaly by conventional treatments; much work has been devoted to approximate methods. Here, we give an explicit integral solution, utilizing methods recently applied to the "geometric goat problem" and to the dynamics of spherical collapse. The solution is given as a ratio of contour integrals; these can be efficiently computed via numerical integration for arbitrary eccentricities. The method is found to be highly accurate in practice, with our C++ implementation outperforming conventional root-finding and series approaches by a factor greater than two.
Submission history
From: Oliver Henry Edward Philcox [view email]
[v1]
Mon, 29 Mar 2021 18:00:00 UTC (1,746 KB)
[v2]
Tue, 4 May 2021 18:20:01 UTC (2,559 KB)
[v3]
Thu, 13 May 2021 23:55:47 UTC (2,561 KB)