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

${\displaystyle 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[1][2]

${\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)}$

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

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