Elliptic integral


In integral calculus, an elliptic integral is one of a number of related functions defined as the value of certain integrals, which were first studied by Giulio Fagnano and Leonhard Euler. Their name originates from their connection with the problem of finding the arc length of an ellipse.
Modern mathematics defines an "elliptic integral" as any function which can be expressed in the form
where is a rational function of its two arguments, is a polynomial of degree 3 or 4 with no repeated roots, and is a constant.
In general, integrals in this form cannot be expressed in terms of elementary functions. Exceptions to this general rule are when has repeated roots, when contains no odd powers of, and when the integral is pseudo-elliptic. However, with the appropriate reduction formula, every elliptic integral can be brought into a form that involves integrals over rational functions and the three Legendre canonical forms, also known as the elliptic integrals of the first, second and third kind.
Besides the Legendre form given below, the elliptic integrals may also be expressed in Carlson symmetric form. Additional insight into the theory of the elliptic integral may be gained through the study of the Schwarz–Christoffel mapping. Historically, elliptic functions were discovered as inverse functions of elliptic integrals.

Argument notation

Incomplete elliptic integrals are functions of two arguments; complete elliptic integrals are functions of a single argument. These arguments are expressed in a variety of different but equivalent ways as they give the same elliptic integral. Most texts adhere to a canonical naming scheme, using the following naming conventions.
For expressing one argument:
Each of the above three quantities is completely determined by any of the others. Thus, they can be used interchangeably.
The other argument can likewise be expressed as, the amplitude, or as or, where and is one of the Jacobian elliptic functions.
Specifying the value of any one of these quantities determines the others. Note that also depends on. Some additional relationships involving include
The latter is sometimes called the delta amplitude and written as. Sometimes the literature also refers to the complementary parameter, the complementary modulus, or the complementary modular angle. These are further defined in the article on quarter periods.
In this notation, the use of a vertical bar as delimiter indicates that the argument following it is the "parameter", while the backslash indicates that it is the modular angle. The use of a semicolon implies that the argument preceding it is the sine of the amplitude:
This potentially confusing use of different argument delimiters is traditional in elliptic integrals and much of the notation is compatible with that used in the reference book by Abramowitz and Stegun and that used in the integral tables by Gradshteyn and Ryzhik.
There are still other conventions for the notation of elliptic integrals employed in the literature. The notation with interchanged arguments,, is often encountered; and similarly for the integral of the second kind. Abramowitz and Stegun substitute the integral of the first kind,, for the argument in their definition of the integrals of the second and third kinds, unless this argument is followed by a vertical bar: i.e. for. Moreover, their complete integrals employ the parameter as argument in place of the modulus, i.e. rather than. And the integral of the third kind defined by Gradshteyn and Ryzhik,, puts the amplitude first and not the "characteristic".
Thus one must be careful with the notation when using these functions, because various reputable references and software packages use different conventions in the definitions of the elliptic functions. For example, Wolfram's Mathematica software and Wolfram Alpha define the complete elliptic integral of the first kind in terms of the parameter, instead of the elliptic modulus.

Incomplete elliptic integral of the first kind

The incomplete elliptic integral of the first kind is defined as
This is Legendre's trigonometric form of the elliptic integral; substituting and, one obtains Jacobi's algebraic form:
Equivalently, in terms of the amplitude and modular angle one has:
With one has:
demonstrating that this Jacobian elliptic function is a simple inverse of the incomplete elliptic integral of the first kind.
The incomplete elliptic integral of the first kind has following addition theorem:
The elliptic modulus can be transformed that way:

Incomplete elliptic integral of the second kind

The incomplete elliptic integral of the second kind in Legendre's trigonometric form is
Substituting and, one obtains Jacobi's algebraic form:
Equivalently, in terms of the amplitude and modular angle:
Relations with the Jacobi elliptic functions include
The meridian arc length from the equator to latitude is written in terms of :
where is the semi-major axis, and is the eccentricity.
The incomplete elliptic integral of the second kind has following addition theorem:
The elliptic modulus can be transformed that way:

Incomplete elliptic integral of the third kind

The incomplete elliptic integral of the third kind is
or
The number is called the characteristic and can take on any value, independently of the other arguments. Note though that the value is infinite, for any.
A relation with the Jacobian elliptic functions is
The meridian arc length from the equator to latitude is also related to a special case of :

Complete elliptic integral of the first kind

Elliptic Integrals are said to be 'complete' when the amplitude and therefore. The complete elliptic integral of the first kind may thus be defined as
or more compactly in terms of the incomplete integral of the first kind as
It can be expressed as a power series
where is the Legendre polynomials, which is equivalent to
where denotes the double factorial. In terms of the Gauss hypergeometric function, the complete elliptic integral of the first kind can be expressed as
The complete elliptic integral of the first kind is sometimes called the quarter period. It can be computed very efficiently in terms of the arithmetic–geometric mean:
Therefore, the modulus can be transformed as:
This expression is valid for all and :

Relation to the gamma function

If and , then is expressible in closed form in terms of the gamma function. For example,, and give, respectively,
and
and
More generally, the condition that
be in an imaginary quadratic field is sufficient. For instance, if, then and
The second formula above, written as, can be completed by 5 equations showing that is a period for all even divisors of :

Asymptotic expressions

This approximation has a relative precision better than for. Keeping only the first two terms is correct to 0.01 precision for.

Differential equation

The differential equation for the elliptic integral of the first kind is
A second solution to this equation is. This solution satisfies the relation

Continued fraction

A continued fraction expansion is:
where the nome is in its definition.

Inverting the period ratio

Here, we use the complete elliptic integral of the first kind with the parameter instead, because the squaring function introduces problems when inverting in the complex plane. So let
and let
be the theta functions.
The equation
can then be solved by
which is in fact the modular lambda function.
For the purposes of computation, the error analysis is given by
where and.
Also
where.

Complete elliptic integral of the second kind

The complete elliptic integral of the second kind is defined as
or more compactly in terms of the incomplete integral of the second kind as
For an ellipse with semi-major axis and semi-minor axis and eccentricity, the complete elliptic integral of the second kind is equal to one quarter of the circumference of the ellipse measured in units of the semi-major axis. In other words:
The complete elliptic integral of the second kind can be expressed as a power series
which is equivalent to
In terms of the Gauss hypergeometric function, the complete elliptic integral of the second kind can be expressed as
The modulus can be transformed that way:

Computation

Like the integral of the first kind, the complete elliptic integral of the second kind can be computed very efficiently using the arithmetic–geometric mean.
Define sequences and, where, and the recurrence relations, hold. Furthermore, define
By definition,
Also
Then
In practice, the arithmetic-geometric mean would simply be computed up to some limit. This formula converges quadratically for all. To speed up computation further, the relation can be used.
Furthermore, if and , then is expressible in closed form in terms of
and hence can be computed without the need for the infinite summation term. For example,, and give, respectively,
and
and

Derivative and differential equation

A second solution to this equation is.

Complete elliptic integral of the third kind

The complete elliptic integral of the third kind can be defined as
Note that sometimes the elliptic integral of the third kind is defined with an inverse sign for the characteristic,
Just like the complete elliptic integrals of the first and second kind, the complete elliptic integral of the third kind can be computed very efficiently using the arithmetic-geometric mean.

Partial derivatives

Jacobi zeta function

In 1829, Jacobi defined the Jacobi zeta function:
It is periodic in with minimal period. It is related to the Jacobi zn function by. In the literature, sometimes means Wikipedia's. Some authors use for both Wikipedia's and.