In orbital mechanics, Kepler's equation relates various geometric properties of the orbit of a body subject to a central force. It was derived by Johannes Kepler in 1609 in Chapter 60 of his Astronomia nova, and in book V of his Epitome of Copernican Astronomy (1621) Kepler proposed an iterative solution to the equation. This equation and its solution, however, first appeared in a 9th-century work by Habash al-Hasib al-Marwazi, which dealt with problems of parallax. The equation has played an important role in the history of both physics and mathematics, particularly classical celestial mechanics.
Equation
Kepler's equation is
where M {\displaystyle M} is the mean anomaly, E {\displaystyle E} is the eccentric anomaly, and e {\displaystyle e} is the eccentricity. The 'eccentric anomaly' E {\displaystyle E} is useful to compute the position of a point moving in a Keplerian orbit. As for instance, if the body passes the periastron at coordinates x = a ( 1 − e ) {\displaystyle x=a(1-e)} , y = 0 {\displaystyle y=0} , at time t = t 0 {\displaystyle t=t_{0}} , then to find out the position of the body at any time, you first calculate the mean anomaly M {\displaystyle M} from the time and the mean motion n {\displaystyle n} by the formula M = n ( t − t 0 ) {\displaystyle M=n(t-t_{0})} , then solve the Kepler equation above to get E {\displaystyle E} , then get the coordinates relative to the central gravitational body from:
where a {\displaystyle a} is the semi-major axis, b {\displaystyle b} the semi-minor axis. Kepler's equation is a transcendental equation because sine is a transcendental function, and it cannot be solved for E {\displaystyle E} algebraically. Numerical analysis and series expansions are generally required to evaluate E {\displaystyle E} .
Alternate forms There are several forms of Kepler's equation. Each form is associated with a specific type of orbit. The standard Kepler equation is used for elliptic orbits ( 0 ≤ e < 1 {\displaystyle 0\leq e<1} ). The hyperbolic Kepler equation is used for hyperbolic trajectories ( e > 1 {\displaystyle e>1} ). The radial Kepler equation is used for linear (radial) trajectories ( e = 1 {\displaystyle e=1} ). Barker's equation is used for parabolic trajectories (for which e = 1 {\displaystyle e=1} ). With the parabolic orbit, unlike the elliptical or hyperbolic orbits, it is possible to solve Barker's equation and find a closed-form expression for the position as a function of time. When e = 0 {\displaystyle e=0} , the orbit is circular. Increasing e {\displaystyle e} causes the circle to become elliptical. When e = 1 {\displaystyle e=1} , there are four possibilities:
a parabolic trajectory, a trajectory that goes back and forth along a line segment from the centre of attraction to a point at some distance away, a trajectory going in or out along an infinite ray emanating from the centre of attraction, with its speed going to zero with distance or a trajectory along a ray, but with speed not going to zero with distance. A value of e {\displaystyle e} slightly above 1 results in a hyperbolic orbit with a turning angle of just under 180 degrees. Further increases reduce the turning angle, and as e {\displaystyle e} goes to infinity, the orbit becomes a straight line of infinite length.
Hyperbolic Kepler equation The Hyperbolic Kepler equation is:
where H {\displaystyle H} is the hyperbolic eccentric anomaly. This equation is derived by redefining M to be the square root of −1 times the right-hand side of the elliptical equation:
M = i ( E − e sin E ) {\displaystyle M=i\left(E-e\sin E\right)}
(in which E {\displaystyle E} is now imaginary) and then replacing E {\displaystyle E} by i H {\displaystyle iH} .
Radial Kepler equations The Radial Kepler equation for the case where the object does not have enough energy to escape is:
where t {\displaystyle t} is proportional to time and x {\displaystyle x} is proportional to the distance from the centre of attraction along the ray and attains the value 1 at the maximum distance. This equation is derived by multiplying Kepler's equation by 1/2 and setting e {\displaystyle e} to 1:
t ( x ) = 1 2 [ E − sin E ] . {\displaystyle t(x)={\frac {1}{2}}\left[E-\sin E\right].}
and then making the substitution
… excerpt ends here. Continue reading the full article.



