Laplace limit

From HandWiki

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 ε:

[math]\displaystyle{ E = M + \sin(M) \, \varepsilon + \tfrac12 \sin(2M) \, \varepsilon^2 + \left( \tfrac38 \sin(3M) - \tfrac18 \sin(M) \right) \, \varepsilon^3 + \cdots }[/math]

or in general[1][2]

[math]\displaystyle{ 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) }[/math]

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 given by the solution to the transcendental equation

[math]\displaystyle{ \frac{x\exp(\sqrt{1+x^2})}{1+\sqrt{1+x^2}}=1. }[/math] [3]

No closed-form expression or infinite series is known for the Laplace limit.

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.[4]

See also

References

  1. Finch (2003), §4.8
  2. Moulton (1914), §99
  3. Math World
  4. Francesco Carlini: Kepler's equation and the asymptotic solution to singular differential equations. Historia Mathematica Volume 53, November 2020, Pages 1-32. https://arxiv.org/pdf/2002.02679.pdf

External links