Laplace limit

Template:Short description 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)\epsilon +\frac{1}{2}\sin (2M){\epsilon }^2+(\frac{3}{8}\sin (3M)-\frac{1}{8}\sin (M)){\epsilon }^3+\cdots$

or in general

$E=M+\sum _{n=1}^{\mathrm{\infty }}\frac{{\epsilon }^n}{2^{n-1}n!}\sum _{k=0}^{\lfloor n/2\rfloor }(-1{)}^k(\genfrac{}{}{0ex}{}{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

• Template:Citation.
• Template:Citation.

External links

• Template:MathWorld
• Template:OEIS el

Template:Mathanalysis-stub Template:Physics-stub