Rapidly convergent series for ellipse perimeter

The previous post looked at two simple approximations for the perimeter of an ellipse. Approximations are necessary since the perimeter of an ellipse cannot be expressed as an elementary function of the axes.

Kepler noted in 1609 that you could approximate the perimeter of an ellipse as the perimeter of a circle whose radius is the mean of the semi-axes of the ellipse, where the mean could be either the arithmetic mean or the geometric mean. The previous post showed that the arithmetic mean is more accurate, and that it under-estimates the perimeter. This post will explain both of these facts.

There are several series for calculating the perimeter of an ellipse. In 1798 James Ivory came up with a series that converges more rapidly than previously discovered series. Ivory's series is

where

If you're not familiar with the !! notation, see this post on multifactorials.

The 0th order approximation using Ivory's series, dropping all the infinite series terms, corresponds to Kepler's approximation using the arithmetic mean of the semi-axes a and b. By convention the semi-major axis is labeled a and the semi-minor axis b, but the distinction is unnecessary here since Ivory's series is symmetric in a and b.

Note that h >= 0 and h = 0 only if the ellipse is a circle. So the terms in the series are positive, which explains why Kepler's approximation under-estimates the perimeter.

Using just one term in Ivory's series gives a very good approximation

The approximation error increases as the ratio of a to b increases, but even for a highly eccentric ellipse like the orbit of the Mars Orbital Mission, the ratio of a to b is only 2, and the relative approximation error is about 1/500, about 12 times more accurate than Kepler's approximation.

The post Rapidly convergent series for ellipse perimeter first appeared on John D. Cook.