Elliptic integral

From Wikipedia, the free encyclopedia

Jump to: navigation, search

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(x) = \int_{c}^{x} R(t,P(t))\ 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. In particular, we have F (sn(z;k);k) = z, where sn is one of Jacobi's elliptic functions.

Contents

[edit] Notation

Elliptic integrals are a function of two arguments. These arguments are expressed in a variety of different but completely equivalent ways (they give the same elliptic integral). Most texts adhere to a canonical naming scheme, using the following naming conventions.

For expressing one

Each of the above three quantities is completely determined by any of the others (given that they are nonnegative). Thus, they can be used interchangeably.

The other argument can likewise be expressed in a number of different ways:

Specifying the value of any one of these quantities determines the others. Note that u also depends on m. Some additional relationships involving u include

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

and

\sqrt{1-m\sin^2 \phi} = \textrm{dn}\; u.\,\!

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

[edit] Incomplete elliptic integral of the first kind

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

 F(\phi\setminus o\!\varepsilon ) = F(\phi|m) =
\int_0^\phi\frac{d\theta}{\sqrt{1-(\sin\theta\sin o\!\varepsilon)^2}}.\,\!

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

 F(\phi\setminus o\!\varepsilon ) = F(x;k) =
\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(\sin\phi;\sin o\!\varepsilon) = F(\phi|\sin (o\!\varepsilon)^2) = F(\phi\setminus o\!\varepsilon )~ \,\!, 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[citation needed]. 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[1] and the Wikipedia article "Legendre form" use  F(\phi,k) \,\!. The notation is equivalent to our  F(\phi|k^2)~ \,\!; also  E(\phi,k)=E(\phi|k^2)~ \,\! 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 0 < x < 1.

Note that

F(x;k) = u \,\!

with u as defined above: thus, the Jacobian elliptic functions are inverses to the elliptic integrals.

[edit] Incomplete elliptic integral of the second kind

The incomplete elliptic integral of the second kind E is

 E(x;k) = \int_0^x \sqrt{1-k^2 \sin^2\theta}  \ d\theta. \,\!

Equivalently, substituting t=\sin\theta:\,\!

 E(x;k) = \int_0^x \frac{\sqrt{1-k^2 t^2} }{\sqrt{1-t^2}}\ dt. \,\!

Equivalently, using an alternate notation:

 E(\phi\setminus o\!\varepsilon) = E(\phi|m) =
\int_0^\phi\!E'(\theta)\ d\theta = \int_0^\phi\sqrt{1-(\sin\theta\sin o\!\varepsilon)^2}\ d\theta.\,\!

Additional relations include

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

[edit] Incomplete elliptic integral of the third kind

The incomplete elliptic integral of the third kind \Pi\,\! is

 \Pi(n; \phi|m) = \int_0^\phi  \frac{1}{1-n\sin^2 \theta}
\frac {d\theta}{\sqrt{1-(\sin\theta\sin o\!\varepsilon)^2}},\,\!

or

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

or

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

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

[edit] 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(k) = \int_0^{\pi/2} \frac{d\theta}{\sqrt{1-k^2 \sin^2\theta}}

or

K(k) = \int_0^1 \frac{dt}{\sqrt{(1-t^2)(1-k^2 t^2)}}.\!

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

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

The special case can be expressed as a power series

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

which is equivalent to

K(k) = \frac{\pi}{2}\left\{1 + \left(\frac{1}{2}\right)^2 k^{2} + \left(\frac{1 \cdot 3}{2 \cdot 4}\right)^2 k^{4} + \cdots + \left[\frac{\left(2n - 1\right)!!}{\left(2n\right)!!}\right]^2 k^{2n} + \cdots \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(k) = \frac{\pi}{2} \,_2F_1 \left(\frac{1}{2}, \frac{1}{2}; 1; k^2\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.

[edit] Special values

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

[edit] The derivative of the complete elliptic integral of the first kind

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

[edit] Complete elliptic integral of the second kind

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

E(k) = \int_0^{\pi/2}\sqrt {1-k^2 \sin^2\theta}\ d\theta\!

or

 E(k) = \int_0^1 \frac{\sqrt{1-k^2 t^2}}{\sqrt{1-t^2}}\ dt.\!

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

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

that can be expressed as a power series

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

which is

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

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

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

[edit] Special values

E(0) = \frac \pi 2
E(1) = 1 \,
E\left(\frac\sqrt{2} 2\right) = \pi^{\frac 3 2} \Gamma\left(\frac 1 4\right)^{-2}+\frac{\Gamma\left(\frac 1 4\right)^2}{8 \sqrt \pi}
E\left(\tfrac{1}{4}(\sqrt{6} - \sqrt{2})\right) = 2^{\frac 1 3} 3^{-\frac 3 4} \pi^2 \Gamma\left(\tfrac 1 3\right)^{-3} + 2^{-\frac {10} 3} 3^{-\frac {1} 4} (\sqrt3 + 1)\pi^{-1} \Gamma\left(\tfrac 1 3\right)^3 \,
E\left(\tfrac{1}{4}(\sqrt{6} + \sqrt{2})\right) = 2^{\frac 1 3} 3^{-\frac 1 4} \pi^2 \Gamma\left(\tfrac 1 3\right)^{-3} + 2^{-\frac {10} 3} 3^{\frac 1 4} (\sqrt3 - 1)\pi^{-1} \Gamma\left(\tfrac 1 3\right)^3 \,

[edit] The derivative of the complete elliptic integral of the second kind

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

[edit] Complete elliptic integral of the third kind

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

\Pi(n,k) = \int_0^{\pi/2}\frac{\ d\theta}{(1-n\sin^2\theta)\sqrt {1-k^2 \sin^2\theta}}

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

\Pi'(n,k) = \int_0^{\frac{\pi}{2}}\frac{\ d\theta}{(1+n\sin^2\theta)\sqrt {1-k^2 \sin^2\theta}}.

[edit] The partial derivatives of the complete elliptic integral of the third kind

\frac{\partial\Pi(n,k)}{\partial n}=
\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{\partial\Pi(n,k)}{\partial k}=
\frac{k}{n-k^2}\left(\frac{E(k)}{k^2-1}+\Pi(n,k)\right)

[edit] See also

[edit] References

  1. ^ ISBN 0122947576 Eq.(8.111)