Elliptic Integrals, Elliptic Functions and

Theta Functions
Reading Problems
Outline
Background . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
Theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
Elliptic Integral . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
Elliptic Function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
Theta Function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
Assigned Problems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
1
Background
This chapter deals with the Legendre elliptic integrals, the Theta functions and the Jaco-
bian elliptic functions. These elliptic integrals and functions nd many applications in the
theory of numbers, algebra, geometry, linear and non-linear ordinary and partial dierential
equations, dynamics, mechanics, electrostatics, conduction and eld theory.
An elliptic integral is any integral of the general form
f(x) =
_
A(x) + B(x)
C(x) + D(x)
_
S(x)
dx
where A(x), B(x), C(x) and D(x) are polynomials in x and S(x) is a polynomial of
degree 3 or 4. Elliptic integrals can be viewed as generalizations of the inverse trigonometric
functions. Within the scope of this course we will examine elliptic integrals of the rst and
second kind which take the following forms:
First Kind
If we let the modulus k satisfy 0 k
2
< 1 (this is sometimes written in terms of the
parameter m k
2
or modular angle sin
1
k). The incomplete elliptic integral of the
rst kind is written as
F(, k) =
_
sin
0
dt
_
(1 t
2
)(1 k
2
t
2
)
, 0 k
2
1 and 0 sin 1
if we let t = sin and dt = cos d =

1 t
2
d, then
F(, k) =
_

0
d
_
1 k
2
sin
2

, 0 k
2
1 and 0 /2
This is referred to as the incomplete Legendre elliptic integral. The complete elliptic integral
can be obtained by setting the upper bound of the integral to its maximum range, i.e.
sin = 1 or = /2 to give
K(k) =
_
1
0
dt
_
(1 t
2
)(1 k
2
t
2
)
=
_
/2
0
d
_
1 k
2
sin
2

2
Second Kind
E(, k) =
_
sin
0

1 k
2
t
2

1 t
2
dt
=
_

0
_
1 k
2
sin
2
d
Similarly, the complete elliptic integral can be obtained by setting the upper bound of inte-
gration to the maximum value to get
E(k) =
_
1
0

1 k
2
t
2

1 t
2
dt
=
_
/2
0
_
1 k
2
sin
2
t dt
Another very useful class of functions can be obtained by inverting the elliptic integrals. As
an example of the Jacobian elliptic function sn we can write
u(x = sin , k) = F(, k) =
_
sin
0
dt
_
(1 t
2
)(1 k
2
t
2
)
If we wish to nd the inverse of the elliptic integral
x = sin = sn(u, k)
or
u =
_
sn
0
dt
_
(1 t
2
)(1 k
2
t
2
)
While there are 12 dierent types of Jacobian elliptic functions based on the number of poles
and the upper limit on the elliptic integral, the three most popular are the copolar trio of sine
amplitude, sn(u, k), cosine amplitude, cn(u, k) and the delta amplitude elliptic function,
dn(u, k) where
3
sn
2
+ cn
2
= 1 and k
2
sn
2
+ dn
2
= 1
Historical Perspective
The rst reported study of elliptical integrals was in 1655 when John Wallis began to study
the arc length of an ellipse. Both John Wallis (1616-1703) and Isaac Newton (1643-1727)
published an innite series expansion for the arc length of the ellipse. But it was not until the
late 1700s that Legendre began to use elliptic functions for problems such as the movement
of a simple pendulum and the deection of a thin elastic bar that these types of problems
could be dened in terms of simple functions.
Adrien-Marie Legendre (1752-1833), a French mathematician, is remembered mainly for the
Legendre symbol and Legendre functions which bear his name but he spent more than forty
years of his life working on elliptic functions, including the classication of elliptic integrals.
His rst published writings on elliptic integrals consisted of two papers in the Memoires de
lAcadmie Francaise in 1786 based on elliptic arcs. In 1792 he presented to the Acadmie a
memoir on elliptic transcendents.
Legendres major work on elliptic functions appeared in 3 volumes
5
in 1811-1816. In the
rst volume Legendre introduced basic properties of elliptic integrals as well as properties
for beta and gamma functions. More results on beta and gamma functions appeared in the
second volume together with applications of his results to mechanics, the rotation of the
Earth, the attraction of ellipsoids and other problems. The third volume contained the very
elaborate and now well-known tables of elliptic integrals which were calculated by Legendre
himself, with an account of the mode of their construction. He then repeated much of this
work again in a three volume set
6
in 1825-1830.
Despite forty years of dedication to elliptic functions, Legendres work went essentially unno-
ticed by his contemporaries until 1827 when two young and as yet unknown mathematicians
Abel and Jacobi placed the subject on a new basis, and revolutionized it completely.
In 1825, the Norwegian government funded Abel on a scholarly visit to France and Germany.
Abel then traveled to Paris, where he gave an important paper revealing the double period-
icity of the elliptic functions. Among his other accomplishments, Abel wrote a monumental
work on elliptic functions
7
which unfortunately was not discovered until after his death.
Jacobi wrote the classic treatise
8
on elliptic functions, of great importance in mathematical
physics, because of the need to integrate second order kinetic energy equations. The motion
5
Exercises du Calcul Intgral
6
Trait des Fonctions Elliptiques
7
Abel, N.H. Recherches sur les fonctions elliptiques. J. reine angew. Math. 3, 160-190, 1828.
8
Jacobi, C.G.J. Fundamentia Nova Theoriae Functionum Ellipticarum. Regiomonti, Sumtibus fratrum
Borntraeger, 1829.
4
equations in rotational form are integrable only for the three cases of the pendulum, the
symmetric top in a gravitational eld, and a freely spinning body, wherein solutions are in
terms of elliptic functions.
Jacobi was also the rst mathematician to apply elliptic functions to number theory, for
example, proving the polygonal number theorem of Pierre de Fermat. The Jacobi theta
functions, frequently applied in the study of hypergeometric series, were named in his honor.
In developments of the theory of elliptic functions, modern authors mostly follow Karl Weier-
strass. The notations of Weierstrasss elliptic functions based on his p-function are conve-
nient, and any elliptic function can be expressed in terms of these. The elliptic functions
introduced by Carl Jacobi, and the auxiliary theta functions (not doubly-periodic), are more
complex but important both for the history and for general theory.
5
Theory
1. Elliptic Integrals
There are three basic forms of Legendre elliptic integrals that will be examined here; rst,
second and third kind. In their most general form, elliptic integrals are presented in a form
referred to as incomplete integrals where the bounds of the integral representation range
from
0 sin 1 or 0 /2.
a) First Kind: The incomplete elliptic integral can be written as
F(sin , k) =
_
sin
0
dt
_
(1 t
2
)(1 k
2
t
2
)
, 0 k
2
1 (3.1)
0 sin 1
by letting t = sin , Eq. 3.1 becomes
F(, k) =
_

0
d
_
(1 k
2
sin
2
)
, 0 k
2
< 1 (3.2)
0 <

2
The parameter k is called the modulus of the elliptic integral and is the amplitude
angle.
The complete elliptic integral is obtained by setting the amplitude = /2 or
sin = 1, the maximum range on the upper bound of integration for the elliptic
integral.
F
_
=

2
, k
_
= F(sin = 1, k) = K(k) = K (3.3)
A complementary form of the elliptical integral can be obtained by letting the modulus
be
(k

)
2
= 1 k
2
(3.4)
6
If we let v = tan and in turn dv = sec
2
d = (1 + v
2
) d, then
F(, k) =
_
tan
0
dv
(1 + v
2
)

_
1 k
2
_
v
2
1 + v
2
__
=
_
tan
0
dv

1 + v
2
_
(1 + v
2
k
2
v
2
=
_
tan
0
dv
_
(1 + v
2
)(1 + k

v
2
)
(3.5)
The complementary, complete elliptic integral can then be written as
F
_
=

2
, k

_
= F(sin = 1, k

) = K(k

) = K

(3.6)
b) Second Kind:
E(, k) =
_
sin
0

1 k
2
t
2
1 t
2
dt, 0 k
2
1 (3.7)
or its equivalent
E(, k) =
_

0
_
1 k
2
sin
2
d, 0 k
2
1 (3.8)
0

2
And similarly, the complete elliptic integral of the second kind can be written as
E
_
=

2
, k
_
= E(sin = 1, k) = E(k) = E (3.9)
and the complementary complete integral of the second kind
E
_
=

2
, k

_
= E(sin = 1, k

) = E(k

) = E

(3.10)
7
c) Third Kind:
(, n, k) =
_
sin
0
dt
(1 + nt)
2
_
(1 t
2
)(1 k
2
t
2
)
, 0 k
2
1
(3.11)
or its equivalent
(, n, k) =
_

0
d
(1 + nsin
2
)
_
(1 k
2
sin
2
)
, 0 k
2
1 (3.12)
0

2
8
Computer Algebra Systems
More than any other special function, you need to be very careful about the arguments
you give to elliptic integrals and elliptic functions. There are several conventions in
common use in competing Computer Algebra Systems and it is important that you
check the documentation provided with the CAS to determine which convection is
incorporated.
Function Mathematica
Elliptic Integral of the rst kind, F[|m] EllipticF[,m]
- amplitude and modulus m = k
2
Complete Elliptic Integral EllipticK[m]
of the rst kind, K(m)
Elliptic Integral of the second kind, E[|m] EllipticE[,m]
- amplitude and modulus m = k
2
Complete Elliptic Integral EllipticE[m]
of the second kind, E(m)
Elliptic Integral of the third kind, [n; |k] EllipticPi[n, ,m]
- amplitude and modulus m = k
2
Complete Elliptic Integral EllipticPi[n,m]
of the third kind, (n|m)
Potential Applications
Determining the arc length of a circle is easily achieved using trigonometric functions,
however elliptic integrals must be used to nd the arc length of an ellipse.
Tracing the arc of a pendulum can be achieved for small angles using trigonometric
functions but to determine the full path of the pendulum elliptic integrals must be
used.
9
Relations and Selected Values of Elliptic Integrals
Complete Elliptic Integrals of the First and Second Kind, K, K

, E, E

The four elliptic integrals K, K

, E, and E

, satisfy the following identity attributed

to Legendre
KE

+ K

E KK

2
(3.13)
The elliptic integrals K and E as functions of the modulus k are connected by means
of the following equations:
dE
dk
=
1
k
(E K) (3.14)
dK
dk
=
1
k(k

)
2
[E (k

)
2
K] (3.15)
Incomplete Elliptic Integrals of the First and Second Kind, F(, k), E(, k)
It is convenient to introduce another frequently encountered elliptic integral which is
related to E and F.
D(, k) =
_

0
sin
2

d =
F E
k
2
(3.16)
where
=
_
1 k
2
sin
2
(3.17)
Therefore
F = E + k
2
D (3.18)
10
Other incomplete integrals expressed by D, E, and F are
_

0
cos
2

d = F D (3.19)
_

0
tan
2

d =
tan E
(k

)
2
(3.20)
_

0
d
cos
2

=
tan + k
2
(DF)
(k

)
2
(3.21)
_

0
sin
2

2
d =
F D
(k

)
2

sin cos
(k

)
2

(3.22)
_

0
cos
2

2
d = D +
sin cos

(3.23)
_

0
tan
2
d = tan + F 2E (3.24)
Special Values of the Elliptic Integrals
E(0, k) = 0 (3.25)
F(0, k) = 0 (3.26)
(0,
2
, k) = 0 (
2
= n) (3.27)
E(, k) = (3.28)
F(, k) = (3.29)
(,
2
, 0) = (if n = 0) (3.30)
=
arctan
_
_
(1
2
) tan
_

1
2
, if
2
< 1 (3.31)
=
arctanh
_
_
(
2
1) tan
_

2
1
, if
2
> 1 (3.32)
11
K(0) = K

(1) = /2 (3.33)
E(0) = E

(1) = /2 (3.34)
E(, 1) = sin (3.35)
F(, 1) = ln(tan + sec ) (3.36)
(,
2
, 1) =
1
1
2
_
ln(tan + sec ) ln

1 + sin
1 sin
_
(3.37)
if
2
> 0,
2
= 1
=
1
1
2
[ln(tan + sec ) +||arctan(|| sin )]
if
2
< 0
Dierentiation and Integration with Respect to the Modulus k
F
k
=
k
(k

)
2
_
F D
sin cos

_
(3.38)
E
k
= kD (3.39)
D
k
=
1
k(k

)
2
_
F(, k) D(, k)
sin cos

D(, k)
k
_
(3.40)
_
F k dk = E(, k) (k

)F(, k) (1 )cotan (3.41)

_
D k dk = E(, k) (3.42)
_
E k dk =
1
3
(1 + k
2
)E(, k) (k

)
2
F(, k) (1 )cotan (3.43)
12
Jacobis Nome q and Complementary Nome q
1
The nome q and the complementary nome q
1
are dened by
q = q(k) = exp[K

/K] (3.44)
and
q
1
= q(k

) = exp[K/K

] (3.45)
Therefore
ln q ln q
1
=
2
(3.46)
To compute the nome q we set k = sin and introduce the parameter dened by
2 =
1

cos
1 +

cos
(3.47)
by means of which we have
q = + 2
5
+ 15
9
+ 150
13
+ 1707
17
(3.48)
The above series converges rapidly provided 0 /4 or 0 k 1/

2. For
/4 < /2 or 1/

2 k 1, let
2
1
=
1

sin
1 +

sin
(3.49)
and compute the complementary nome q
1
by means of
q
1
=
1
+ 2
5
1
+ 15
9
1
+ 150
13
1
+ 1707
17
1
(3.50)
13
2. Elliptic Functions
There are several types of elliptic functions including the Weierstrass elliptic functions
as well as related theta functions but the most common elliptic functions are the
Jacobian elliptic functions, based on the inverses of the three types of elliptic integrals.
1. Jacobi elliptic functions: The three standard forms of Jacobi elliptic integrals
are denoted as sn(u, k), cn(u, k) and dn(u, k) and are the sine, cosine and
delta amplitude elliptic functions, respectively. These functions are obtained by
inverting the elliptic integral of the rst kind where
u = F(, k) =
_

0
d
_
1 k
2
sin
2

(3.51)
where 0 < k
2
< 1 and k is referred to as the elliptic modulus of u and , the
upper bound on the elliptic integral is referred to as the Jacobi amplitude (amp).
The inversion of the elliptic integral gives
= F
1
(u, k) = amp(u, k) (3.52)
and from this we can write
sin = sin(amp(u, k)) = sn(u, k) (3.53)
cos = cos(amp(u, k)) = cn(u, k) (3.54)
_
1 k
2
sin
2
=
_
1 k
2
sin
2
(amp(u, k)) = dn(u, k) (3.55)
These functions are doubly periodic generalizations of the trigonometric functions
satisfying
sn(u, 0) = sin u (3.56)
cn(u, 0) = cos u (3.57)
dn(u, 0) = 1 (3.58)
14
0 1 2 3 4 5 6 7
u
1
0.5
0
0.5
1
s
n

dn u
snu
cnu
0
K 2 K
3 K
4 K
Figure 3.1: Plot of Jacobian Elliptic Functions sn(u), cn(u) and dn(u) for k = 1/2
In total there are 12 Jacobian elliptic functions, where the remaining 9 can be
related to the 3 we have already seen as follows:
cd(u) =
cn(u)
dn(u)
dc(u) =
dn(u)
cn(u)
ns(u) =
1
sn(u)
sd(u) =
sn(u)
dn(u)
nc(u) =
1
cn(u)
ds(u) =
dn(u)
sn(u)
nd(u) =
1
dn(u)
sc(u) =
sn(u)
cn(u)
cs(u) =
cn(u)
sn(u)
2. Weierstrass elliptic functions: The principal dierence between the Jacobi
and the Weierstrass elliptic integrals is in the number of poles in each fundamental
cell. While the Jacobi elliptic functions has two simple poles per cell and can be
considers as a solution to the dierential equation
d
2
x
dt
2
= A + Bx + Cx
2
+ Dx
3
the Weierstrass elliptic function has just one double pole and is a solution to
d
2
x
dt
2
= A + Bx + Cx
2
15
We will focus primarily on the Jacobi elliptic function in this course but you should be
aware of the Weierstrass form of the function.
3. Theta Functions
Theta functions are the elliptic analogs of the exponential function and are typically
written as,
a
(u, q) where a ranges from 1 to 4 to represent the fours variations of
the theta function, u is the argument of the function and q is the Nome, given as
q = e
it
= e
K

/K
(3.59)
where
t = i
K

(k)
K(k)
16
Assigned Problems
Problem Set for Elliptic Integrals and Functions
1. Determine the perimeter of the ellipse
4x
2
+ 9y
2
= 36
2. Obtain the expression for the length of arc of the ellipse
x
2
+
y
2
4
= 1
between (0, 2) and (1/2,

3). Note that b > a in this problem. Compute the arc

length to three decimal places.
3. Compute to four decimal places the following integrals
i)
_
/4
0
dt
_
1
1
2
sin
2
t
ii)
_
/3
0
_
1
3
4
sin
2
t dt
iii)
_
/2
0
dt
_
1
1
15
sin
2
t
iv)
_
/2
0
_
1
80
81
sin
2
t dt
4. By means of the substitution x = 2 sin , show that
_
2
0
dx
_
(4 x
2
)(9 x
2
)
=
1
3
K
_
2
3
_
5. Prove that
_
/2
0
dx

sin x
=
_
/2
0
dx

cos x
=

2 K
_
1

2
_
17
6. Given q = 1/2, compute k, K and K

to six decimal places.

7. Compute to three decimal places the area of the ellipsoid with semi-axes 3, 2, and 1.
8. The thermal constriction resistance on an isothermal, elliptical disk (a > b) on an
insulated isotropic half-space of thermal conductivity k is
R

= kaR = K
_
_

1
_
b
a
_
2
_
_
Compute R

to four decimal places for the following values of a/b : 1, 1.5, 2, 3, 4 and
5. Use the arithmetic-geometric mean method of Gauss.
9. Derivatives of the elliptic integrals. Show that
i)
dE
dk
=
E K
k
ii)
d
2
E
dk
2
=
1
k
dK
dk
=
E (k

)
2
K
k
2
(k

)
2
10. Integrals of the elliptic integrals. Show that
i)
_
K dk =
k
2
_
1 +

n=1
[(2n)!]
2
k
2n
(2n + 1)2
4n
(n!)
4
_
ii)
_
kKL dk = E (k

)
2
K
iii)
_
k E dk =
1
3
[(1 + k
2
)E (k

)
2
K]
18
11. Show that
i)
_
1
0
kE dk
k

=
_
1
0
E(u

) du =

2
8
ii)
_
1
0
K dk
1 + k
=

2
8
12. Derivatives of elliptic integrals with respect to the argument
i)
d
d
F(, k) =
1
_
1 k
2
sin
2

ii)
d
d
E(, k) =
_
1 k
2
sin
2

19
References
1. Abramowitz, M. and Stegun, I.A., (Eds.), Gamma (Factorial) Function and
Incomplete Gamma Function. 6.1 and 6.5 in Handbook of Mathematical Functions
and Formulas, Graphs and Mathematical Tables, 9th printing, Dover, New York, 1972,
pp. 255-258 and pp 260-263.
2. Anderson, D., Vamanamurthy, K. and Vuorinen, M., Functional Inequalities
for Complete Elliptic Integrals, SIAM J. Math. Anal., 21, 1990, pp. 536-549.
3. Anderson, D., Vamanamurthy, K. and Vuorinen, M., Functional inequalities
for hypergeometric functions and complete elliptic integrals, SIAM J. Math. Anal.,
23, 1992, pp. 512-524.
4. Arfken, G., Mathematical Methods for Physicists, 3rd ed., Section 5.8, Elliptic
Integrals, Orlando, FL, Academic Press, pp.321-327, 1985.
5. Borwein, J.M. and Borwein, P.B. Pi, A Study in Analytic Number Theory and
Computational Complexity, Wiley, New York, 1987.
6. Byrd, P.F. and Friedman, M.D., Handbook of Elliptic Integrals for Engineers and
Scientists, Second ed., Springer-Verlag, New York, 1971.
7. Carlson, B.C., Some Inequalities for Hypergeometric Functions, Proc. Amer.
Math. Soc., 17, 1966, pp. 32-39.
8. Carlson, B.C., Inequalities for a Symmetric Elliptic Integral, Proc. Amer. Math.
Soc., 25, 1970, pp. 698-703.
9. Carlson, B.C. and Gustafson, J.L., Asymptotic Expansion of the First Elliptic
Integral, SIAM J. Math. Anal., 16, 1985, pp. 1072-1092.
10. Cayley, A., Elliptic Functions, Second ed., Dover Publications, New York, 1961.
11. Hancock H., Elliptic Integrals, Dover Publications, New York, 1917.
12. Karman, T. von and Biot, M.A., Mathematical Methods in Engineering: An
Introduction to the Mathematical Treatment of Engineering Problems, McGraw-Hill,
New York, 1940.
13. King, L.V., The Direct Numerical Calculation of Elliptic Functions and Integrals,
Cambridge University Press, London, 1924.
14. Prasolov, V. and Solovyev, Y., Elliptic Functions and Elliptic Integrals, Amer.
Math. Soc., Providence, RI, 1997.
15. Weisstein, E.W., Books about Elliptic Integrals,
http://www.ericweisstein.com/encyclopedias/books/EllipticIntegrals.html.
16. Whittaker, E.T. and Watson, G.N. A Course in Modern Analysis, 4th ed., Cam-
bridge University Press, Cambridge, England, 1990.
20