An elliptic integral is an integral of the form
(1)
|
or
(2)
|
where ,
,
, and
are polynomials in
, and
is a polynomial of degree
3 or 4. Stated more simply, an elliptic integral is an integral of
the form
(3)
|
where
is a rational function of
and
,
is a function of
that is cubic or quartic
in
,
contains at least one odd power
of
,
and
has no repeated factors (Abramowitz and Stegun 1972, p. 589).
Elliptic integrals can be viewed as generalizations of the inverse trigonometric functions and provide solutions to a wider class of problems. For instance, while the arc length of a circle is given as a simple function of the parameter, computing the arc length of an ellipse requires an elliptic integral. Similarly, the position of a pendulum is given by a trigonometric function as a function of time for small angle oscillations, but the full solution for arbitrarily large displacements requires the use of elliptic integrals. Many other problems in electromagnetism and gravitation are solved by elliptic integrals.
A very useful class of functions known as elliptic functions is obtained by inverting elliptic integrals to obtain generalizations of the trigonometric functions. Elliptic functions (among which the Jacobi elliptic functions and Weierstrass elliptic function are the two most common forms) provide a powerful tool for analyzing many deep problems in number theory, as well as other areas of mathematics.
All elliptic integrals can be written in terms of three "standard" types. To see this, write
(4)
| |||
(5)
|
But since ,
(6)
| |||
(7)
|
then
(8)
| |||
(9)
| |||
(10)
|
so
(11)
| |||
(12)
|
But any function
can be evaluated in terms of elementary functions,
so the only portion that need be considered is
(13)
|
Now, any quartic can be expressed as where
(14)
| |||
(15)
|
The coefficients here are real, since pairs of complex roots are complex conjugates
(16)
| |||
(17)
|
If all four roots are real, they must be arranged so as not to interleave (Whittaker and Watson 1990, p. 514). Now define a quantity
such that
(18)
|
is a square number and
(19)
|
(20)
|
Call the roots of this equation and
, then
(21)
| |||
(22)
| |||
(23)
| |||
(24)
| |||
(25)
| |||
(26)
|
(27)
| |||
(28)
|
Solving gives
(29)
| |||
(30)
| |||
(31)
| |||
(32)
|
so we have
(33)
|
Now let
(34)
| |||
(35)
| |||
(36)
| |||
(37)
|
so
(38)
| |||
(39)
|
and
(40)
| |||
(41)
| |||
(42)
|
Now let
(43)
|
so
(44)
|
Rewriting the even and odd parts
(45)
| |||
(46)
|
gives
(47)
| |||
(48)
|
so we have
(49)
|
Letting
(50)
| |||
(51)
|
reduces the second integral to
(52)
|
which can be evaluated using elementary functions. The first integral can then be reduced by integration
by parts to one of the three Legendre elliptic integrals (also called Legendre-Jacobi
elliptic integrals), known as incomplete elliptic
integrals of the first, second,
and third kinds, denoted ,
, and
, respectively (von Kármán and Biot
1940, Whittaker and Watson 1990, p. 515). If
, then the integrals are called complete elliptic integrals
and are denoted
,
,
.
Incomplete elliptic integrals are denoted using a elliptic modulus ,
parameter
, or modular angle
. An elliptic integral
is written
when the parameter is used,
when the elliptic modulus
is used, and
when the modular angle is used. Complete elliptic
integrals are defined when
and can be expressed using the expansion
(53)
|
An elliptic integral in standard form
(54)
|
where
(55)
|
can be computed analytically (Whittaker and Watson 1990, p. 453) in terms of the Weierstrass elliptic function with invariants
(56)
| |||
(57)
|
If
is a root of
,
then the solution is
(58)
|
For an arbitrary lower bound,
(59)
|
where
is a Weierstrass elliptic function
(Whittaker and Watson 1990, p. 454).
A generalized elliptic integral can be defined by the function
(60)
| |||
(61)
|
(Borwein and Borwein 1987). Now let
(62)
| |||
(63)
|
But
(64)
|
so
(65)
| |||
(66)
| |||
(67)
| |||
(68)
|
and
(69)
|
and the equation becomes
(70)
| |||
(71)
|
Now we make the further substitution . The differential becomes
(72)
|
but ,
so
(73)
|
(74)
|
and
(75)
|
However, the left side is always positive, so
(76)
|
and the differential is
(77)
|
We need to take some care with the limits of integration. Write (◇) as
(78)
|
Now change the limits to those appropriate for the integration
(79)
|
so we have picked up a factor of 2 which must be included. Using this fact and plugging (◇) in (◇) therefore gives
(80)
|
Now note that
(81)
| |||
(82)
| |||
(83)
|
Plug (◇) into (◇) to obtain
(84)
| |||
(85)
|
But
(86)
| |
(87)
| |
(88)
| |
(89)
|
so
(90)
|
and (◇) becomes
(91)
| |||
(92)
|
We have therefore demonstrated that
(93)
|
We can thus iterate
(94)
| |||
(95)
|
as many times as we wish, without changing the value of the integral. But this iteration is the same as and therefore converges to the arithmetic-geometric
mean, so the iteration terminates at , and we have
(96)
| |||
(97)
| |||
(98)
| |||
(99)
| |||
(100)
|
Complete elliptic integrals arise in finding the arc length of an ellipse and the period of a pendulum. They also arise in a natural way from the theory of theta functions. Complete elliptic integrals can be computed using a procedure involving the arithmetic-geometric mean. Note that
(101)
| |||
(102)
| |||
(103)
|
So we have
(104)
| |||
(105)
|
where
is the complete elliptic integral
of the first kind. We are free to let
and
, so
(106)
|
since ,
so
(107)
|
But the arithmetic-geometric mean is defined by
(108)
| |||
(109)
| |||
(110)
|
where
(111)
|
so we have
(112)
|
where
is the value to which
converges. Similarly, taking instead
and
gives
(113)
|
Borwein and Borwein (1987) also show that defining
(114)
| |||
(115)
|
leads to
(116)
|
so
(117)
|
for
and
,
and
(118)
|
The elliptic integrals satisfy a large number of identities. The complementary functions and moduli are defined by
(119)
|
Use the identity of generalized elliptic integrals
(120)
|
to write
(121)
| |||
(122)
| |||
(123)
|
(124)
|
Define
(125)
|
and use
(126)
|
so
(127)
|
Now letting
gives
(128)
|
(129)
|
(130)
| |||
(131)
| |||
(132)
| |||
(133)
|
and
(134)
| |||
(135)
| |||
(136)
|
Writing
instead of
,
(137)
|
Similarly, from Borwein and Borwein (1987),
(138)
|
(139)
|
Expressions in terms of the complementary function can be derived from interchanging the moduli and their complements in (◇), (◇), (◇), and (◇).
(140)
| |||
(141)
| |||
(142)
| |||
(143)
| |||
(144)
| |||
(145)
|
and
(146)
|
(147)
|
Taking the ratios
(148)
|
gives the modular equation of degree 2. It is also true that
(149)
|