Kepler equation Newton iteration

What you are seeing: the orbit of a planet in a Keplerian ellipse (left) and the Newton-iteration history for the Kepler equation $M = E - e \sin E$ (right). At each rAF frame the mean anomaly $M = 2\pi t / T$ ticks forward; the solver iterates on $E$ until $|E_{n+1} - E_n| \lt 10^{-12}$ and the planet is placed at $(a(\cos E - e), b \sin E)$ on the ellipse.

The right panel shows successive $|E_n - E_\infty|$ on a log axis. Newton's method has quadratic local convergence: each iteration roughly squares the remaining error. For moderate eccentricity ($e \lesssim 0.9$) the residual drops to machine precision in 4-6 iterations; near $e = 0.99$ it takes 10-15.