Definitions

Elliptic integral

In integral calculus, elliptic integrals originally arose in connection with the problem of giving the arc length of an ellipse. They were first studied by Giulio Fagnano and Leonhard Euler. Modern mathematics defines an elliptic integral as any function f which can be expressed in the form

$f\left(x\right) = int_\left\{c\right\}^\left\{x\right\} R\left(t,P\left(t\right)\right) dt ,!$

where R is a rational function of its two arguments, P is the square root of a polynomial of degree 3 or 4 with no repeated roots, and c is a constant.

In general, elliptic integrals cannot be expressed in terms of elementary functions. Exceptions to this general rule are when P has repeated roots, or when R(x,y) contains no odd powers of y. However, with the appropriate reduction formula, every elliptic integral can be brought into a form that involves integrals over rational functions and the three canonical forms (i.e. the elliptic integrals of the first, second and third kind).

Besides the forms given below, the elliptic integrals may also be expressed in Legendre form and 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, and this one in particular: we have F (sn(z;k);k) = z where sn is one of Jacobi's elliptic functions.

Notation

Elliptic integrals are often expressed as functions for a variety of different arguments. These different arguments are completely equivalent (they give the same elliptic integral), but they can be confusing due to their different appearance. Most texts adhere to a canonical naming scheme. Before defining the integrals, we review the naming conventions for the arguments:

• $o!varepsilon,,!$ the modular angle;
• $k=sin o!varepsilon,,!$ the elliptic modulus;
• $m=k^2=\left[sin\left(o!varepsilon\right)\right]^2,,!$ the parameter;

Note that the above three conventions are completely determined by one another. Specifying one is the same as specifying another. The elliptic integrals also depend on another argument, which can be specified in a number of different ways:

• $phi,!$ the amplitude
• x where $x=sin phi= textrm\left\{sn\right\} ; u,!$
• u, where x = sn u and sn is one of the Jacobian elliptic functions

Specifying any one of these arguments determines the others. Thus, they may be used interchangeably in the notation. Note that u also depends on m. Some additional relationships involving u include

$cos phi = textrm\left\{cn\right\}; u,!$

and

$sqrt\left\{1-msin^2 phi\right\} = textrm\left\{dn\right\}; u.,!$

The latter is sometimes called the delta amplitude and written as $Delta\left(phi\right)=textrm\left\{dn\right\}; u,!$. 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.

Incomplete elliptic integral of the first kind

The incomplete elliptic integral of the first kind F is defined as

$F\left(phisetminus o!varepsilon \right) = F\left(phi|m\right) =$
int_0^phifrac{dtheta}{sqrt{1-(sinthetasin o!varepsilon)^2}}.,!

Equivalently, using notation in Jacobi's form, one sets $x=sin phi ~,~ t=sin theta;!$; then

$F\left(phisetminus o!varepsilon \right) = F\left(x;k\right) =$
int_{0}^{x} frac{dt}{sqrt{(1-t^2)(1-k^2 t^2)} },! where it is understood that when there is a vertical bar used, the argument following the vertical bar is the parameter (as defined above); and, when a backslash is used, it is followed by the modular angle. In this sense, $F\left(sinphi;sin o!varepsilon\right) = F\left(phi|sin \left(o!varepsilon\right)^2\right) = F\left(phisetminus o!varepsilon \right)~ ,!$, with the notations directly borrowed from the reference book of standards, Abramowitz and Stegun. The use of the delimiters ; | is traditional in elliptic integrals.

However, there remain different conventions for the notation of elliptic integrals! The differences can be very confusing, especially to a novice . The functions that evaluate the elliptic integrals do not have standard and unique names and meanings (like sqrt, sin and erf have). Even the literatures on the subject use differentiated notations. Gradstein, Ryzhik , Eq.(8.111)] and the Wikipedia article "Legendre form" use $F\left(phi,k\right) ,!$. The notation is equivalent to our $F\left(phi|k^2\right)~ ,!$; also $E\left(phi,k\right)=E\left(phi|k^2\right)~ ,!$ below.

Accordingly, if one translates the code from the Mathematica language into the language used by Maple, one should replace the argument of the EllipticK function with its square root. Correspondingly, in the translation from Maple to Mathematica, the argument should be replaced by its square. EllipticK(x) in Maple is almost equivalent to EllipticK[x^2] in Mathematica; one may expect to get the same result in both systems, at least while 0Note that

$F\left(x;k\right) = u ,!$
with u as defined above: thus, the Jacobian elliptic functions are inverses to the elliptic integrals.

Incomplete elliptic integral of the second kind

The incomplete elliptic integral of the second kind E is

$E\left(phisetminus o!varepsilon\right) = E\left(phi|m\right) =$
int_0^phi!E'(theta) dtheta = int_0^phisqrt{1-(sinthetasin o!varepsilon)^2} dtheta.,!

Equivalently, using an alternate notation (substituting $t=sintheta,!$),

$E\left(x;k\right) = int_\left\{0\right\}^\left\{x\right\} frac\left\{sqrt\left\{1-k^2 t^2\right\} \right\}\left\{sqrt\left\{1-t^2\right\}\right\} dt. ,!$

$E\left(phi|m\right) = int_0^u textrm\left\{dn\right\}^2 w ;dw =$
u-mint_0^u textrm{sn}^2 w ;dw = (1-m)u+mint_0^u textrm{cn}^2 w ;dw.,!

Incomplete elliptic integral of the third kind

The incomplete elliptic integral of the third kind $Pi,!$ is
$Pi\left(n; phi|m\right) = int_0^phi frac\left\{1\right\}\left\{1-nsin^2 theta\right\}$
frac {dtheta}{sqrt{1-(sinthetasin o!varepsilon)^2}},,!

or

$Pi\left(n; phi|m\right) = int_\left\{0\right\}^\left\{x\right\} frac\left\{1\right\}\left\{1-nt^2\right\}$
frac{dt}{sqrt{(1-k^2 t^2)(1-t^2) }},,!

or

$Pi\left(n; phi|m\right) = int_0^u frac\left\{dw\right\}\left\{1-n textrm\left\{sn\right\}^2 \left(w|m\right)\right\}. ; ,!$

The number n is called the characteristic and can take on any value, independently of the other arguments. Note though that the value $Pi\left(1;pi/2|m\right),!$ is infinite, for any $m,!$.

Complete elliptic integral of the first kind

Elliptic Integrals are said to be 'complete' when the amplitude is pi/2 and thus x=1. The complete elliptic integral of the first kind K may be defined as

$K\left(k\right) = int_0^\left\{frac\left\{pi\right\}\left\{2\right\}\right\} frac\left\{dtheta\right\}\left\{sqrt\left\{1-k^2 sin^2theta\right\}\right\}$

or

$K\left(k\right) = int_\left\{0\right\}^\left\{1\right\} frac\left\{dt\right\}\left\{sqrt\left\{\left(1-t^2\right)\left(1-k^2 t^2\right)\right\}\right\}.!$

It is a special case of the incomplete elliptic integral of the first kind:

$K\left(k\right) = F\left(1;,k\right) = F\left(frac\left\{pi\right\}\left\{2\right\},|,k^2\right)!$

The special case can be expressed as a power series

$K\left(k\right) = frac\left\{pi\right\}\left\{2\right\} sum_\left\{n=0\right\}^\left\{infty\right\} left\left[frac\left\{\left(2n\right)!\right\}\left\{2^\left\{2 n\right\} n!^2\right\}right\right]^2 k^\left\{2n\right\}!$

which is equivalent to

$K\left(k\right) = frac\left\{pi\right\}\left\{2\right\}left\left\{1 + left\left(frac\left\{1\right\}\left\{2\right\}right\right)^2 k^\left\{2\right\} + left\left(frac\left\{1 cdot 3\right\}\left\{2 cdot 4\right\}right\right)^2 k^\left\{4\right\} + cdots + left\left[frac\left\{left\left(2n - 1right\right)!!\right\}\left\{left\left(2nright\right)!!\right\}right\right]^2 k^\left\{2n\right\} + cdots right\right\}.!$

where $n!!$ denotes the Double Factorial. In terms of the Gauss hypergeometric function, the complete elliptic integral of the first kind can be expressed as

$K\left(k\right) = frac\left\{pi\right\}\left\{2\right\} ,_2F_1 left\left(frac\left\{1\right\}\left\{2\right\}, frac\left\{1\right\}\left\{2\right\}; 1; k^2right\right).,!$

The complete elliptic integral of the first kind is sometimes called the quarter period. It can be computed in terms of the arithmetic-geometric mean.

Special values

$K\left(0\right) = frac 1 2 pi!$
$K\left(1\right) = infty!$
$Kleft\left(fracsqrt\left\{2\right\} 2right\right) = frac\left\{Gammaleft\left(frac 1 4right\right)^2\right\}\left\{4 sqrt\left\{pi\right\}\right\}!$
$Kleft\left(frac\left\{sqrt\left\{6\right\} - sqrt\left\{2\right\}\right\}\left\{4\right\}right\right) = frac\left\{2^\left\{-frac 7 3\right\} 3^\left\{frac 1 4\right\} Gammaleft\left(frac 1 3right\right)^3\right\}\left\{pi\right\}!$
$Kleft\left(frac\left\{sqrt\left\{6\right\} + sqrt\left\{2\right\}\right\}\left\{4\right\}right\right) = frac\left\{2^\left\{-frac 7 3\right\} 3^\left\{frac 3 4\right\} Gammaleft\left(frac 1 3right\right)^3\right\}\left\{pi\right\}!$

The derivative of the complete elliptic integral of the first kind

$frac\left\{mathrm\left\{d\right\}K\left(k\right)\right\}\left\{mathrm\left\{d\right\}k\right\} = frac\left\{E\left(k\right)\right\}\left\{k\left(1-k^2\right)\right\}-frac\left\{K\left(k\right)\right\}\left\{k\right\}$

Complete elliptic integral of the second kind

The complete elliptic integral of the second kind E may be defined as

$E\left(k\right) = int_0^\left\{frac\left\{pi\right\}\left\{2\right\}\right\}sqrt \left\{1-k^2 sin^2theta\right\} dtheta!$

or

$E\left(k\right) = int_\left\{0\right\}^\left\{1\right\} frac\left\{sqrt\left\{1-k^2 t^2\right\}\right\}\left\{sqrt\left\{1-t^2\right\}\right\} dt.!$

It is a special case of the incomplete elliptic integral of the second kind:

$E\left(k\right) = E\left(1;,k\right) = E\left(frac\left\{pi\right\}\left\{2\right\},|,k^2\right)!$

that can be expressed as a power series

$E\left(k\right) = frac\left\{pi\right\}\left\{2\right\} sum_\left\{n=0\right\}^\left\{infty\right\} left\left[frac\left\{\left(2n\right)!\right\}\left\{2^\left\{2 n\right\} n!^2\right\}right\right]^2 frac\left\{k^\left\{2n\right\}\right\}\left\{1-2 n\right\}!$

which is

$E\left(k\right) = frac\left\{pi\right\}\left\{2\right\}left\left\{1 - left\left(frac\left\{1\right\}\left\{2\right\}right\right)^2 frac\left\{k^2\right\}\left\{1\right\} - left\left(frac\left\{1 cdot 3\right\}\left\{2 cdot 4\right\}right\right)^2 frac\left\{k^4\right\}\left\{3\right\} - cdots - left\left[frac\left\{left\left(2n - 1right\right)!!\right\}\left\{left\left(2nright\right)!!\right\}right\right]^2 frac\left\{k^\left\{2n\right\}\right\}\left\{2 n-1\right\} - cdots right\right\}.!$

In terms of the Gauss hypergeometric function, the complete elliptic integral of the second kind can be expressed as

$E\left(k\right) = frac\left\{pi\right\}\left\{2\right\} ,_2F_1 left\left(frac\left\{1\right\}\left\{2\right\}, -frac\left\{1\right\}\left\{2\right\}; 1; k^2right\right).,!$

Special values

$E\left(0\right) = frac pi 2!$

$E\left(1\right) = 1!$

$Eleft\left(fracsqrt\left\{2\right\} 2right\right) = pi^\left\{frac 3 2\right\} Gammaleft\left(frac 1 4right\right)^\left\{-2\right\}+frac\left\{Gammaleft\left(frac 1 4right\right)^2\right\}\left\{8 sqrt pi\right\}!$

$Eleft\left(frac\left\{sqrt\left\{6\right\} - sqrt\left\{2\right\}\right\}\left\{4\right\}right\right) = 2^\left\{frac 1 3\right\} 3^\left\{-frac 3 4\right\} pi^2 Gammaleft\left(frac 1 3right\right)^\left\{-3\right\} + 2^\left\{-frac \left\{10\right\} 3\right\} 3^\left\{-frac \left\{1\right\} 4\right\} frac\left\{sqrt3 + 1\right\}\left\{pi\right\} Gammaleft\left(frac 1 3right\right)^3!$

$Eleft\left(frac\left\{sqrt\left\{6\right\} + sqrt\left\{2\right\}\right\}\left\{4\right\}right\right) = 2^\left\{frac 1 3\right\} 3^\left\{-frac 1 4\right\} pi^2 Gammaleft\left(frac 1 3right\right)^\left\{-3\right\} + 2^\left\{-frac \left\{10\right\} 3\right\} 3^\left\{frac 1 4\right\} frac\left\{sqrt3 - 1\right\}\left\{pi\right\} Gammaleft\left(frac 1 3right\right)^3!$

The derivative of the complete elliptic integral of the second kind

$frac\left\{mathrm\left\{d\right\}E\left(k\right)\right\}\left\{mathrm\left\{d\right\}k\right\}=frac\left\{E\left(k\right)-K\left(k\right)\right\}\left\{k\right\}$

Complete elliptic integral of the third kind

The complete elliptic integral of the third kind $Pi$ can be defined as

$Pi\left(n,k\right) = int_0^\left\{frac\left\{pi\right\}\left\{2\right\}\right\}frac\left\{ dtheta\right\}\left\{\left(1-nsin^2theta\right)sqrt \left\{1-k^2 sin^2theta\right\}\right\}$

Note that sometimes the elliptic integral of the third kind is defined with an inverse sign in $n$, i.e.

$Pi\text{'}\left(n,k\right) = int_0^\left\{frac\left\{pi\right\}\left\{2\right\}\right\}frac\left\{ dtheta\right\}\left\{\left(1+nsin^2theta\right)sqrt \left\{1-k^2 sin^2theta\right\}\right\}.$

The partial derivatives of the complete elliptic integral of the third kind

$frac\left\{partialPi\left(n,k\right)\right\}\left\{partial n\right\}=$
frac{1}{2(k^2-n)(n-1)}left(E(k)+frac{(k^2-n)K(k)}{n}+frac{(n^2-k^2)Pi(n,k)}{n}right)

$frac\left\{partialPi\left(n,k\right)\right\}\left\{partial k\right\}=$
frac{k}{n-k^2}left(frac{E(k)}{k^2-1}+Pi(n,k)right)