Laplace limit

In mathematics, the Laplace limit is the maximum value of the eccentricity for which a solution to Kepler's equation, in terms of a power series in the eccentricity, converges. It is approximately

0.66274 34193 49181 58097 47420 97109 25290.

Kepler's equation M = E − ε sin E relates the mean anomaly M with the eccentric anomaly E for a body moving in an ellipse with eccentricity ε. This equation cannot be solved for E in terms of elementary functions, but the Lagrange reversion theorem gives the solution as a power series in ε:

$E=M+\sin(M)\,\varepsilon +{\tfrac {1}{2}}\sin(2M)\,\varepsilon ^{2}+\left({\tfrac {3}{8}}\sin(3M)-{\tfrac {1}{8}}\sin(M)\right)\,\varepsilon ^{3}+\cdots$

or in general

$E=M\;+\;\sum _{n=1}^{\infty }{\frac {\varepsilon ^{n}}{2^{n-1}\,n!}}\sum _{k=0}^{\lfloor n/2\rfloor }(-1)^{k}\,{\binom {n}{k}}\,(n-2k)^{n-1}\,\sin((n-2k)\,M)$

Laplace realized that this series converges for small values of the eccentricity, but diverges for any value of M other than a multiple of π if the eccentricity exceeds a certain value that does not depend on M. The Laplace limit is this value. It is the radius of convergence of the power series.

It is the unique real solution of the transcendental equation

${\frac {x\exp({\sqrt {1+x^{2}}})}{1+{\sqrt {1+x^{2}}}}}=1$

A closed-form expression in terms of r-Lambert special function and an infinite series representation were given by István Mező.

History

Laplace calculated the value 0.66195 in 1827. The Italian astronomer Francesco Carlini found the limit 0.66 five years before Laplace. Cauchy in the 1829 gave the precise value 0.66274.

See also

• Orbital eccentricity

References

• Finch, Steven R. (2003), "Laplace limit constant", Mathematical constants, Cambridge University Press, ISBN 978-0-521-81805-6.
• Moulton, Forest R. (1914), "V. The Problem of Two Bodies", An Introduction to Celestial Mechanics (2d ed.), MacMillan.

External links

• Weisstein, Eric W. "Laplace Limit". MathWorld.
• OEIS sequence A033259 (Decimal expansion of Laplace's limit constant)