What's New
About the Project
NIST
19 Elliptic IntegralsSymmetric Integrals

§19.29 Reduction of General Elliptic Integrals

Contents

§19.29(i) Reduction Theorems

These theorems reduce integrals over a real interval (y,x) of certain integrands containing the square root of a quartic or cubic polynomial to symmetric integrals over (0,) containing the square root of a cubic polynomial (compare §19.16(i)). Let

19.29.1 Xα =aα+bαx,
Yα =aα+bαy,
x>y, 1α5,
19.29.2 dαβ=aαbβ-aβbα,
dαβ0 if αβ,

and assume that the line segment with endpoints aα+bαx and aα+bαy lies in \(-,0) for 1α4. If

19.29.3 s(t)=α=14aα+bαt

and α,β,γ,δ is any permutation of the numbers 1,2,3,4, then

19.29.4 yxdts(t)=2RF(U122,U132,U232),

where

19.29.5 Uαβ =(XαXβYγYδ+YαYβXγXδ)/(x-y),
Uαβ =Uβα=Uγδ=Uδγ,
Uαβ2-Uαγ2 =dαδdβγ.

There are only three distinct U’s with subscripts 4, and at most one of them can be 0 because the d’s are nonzero. Then

19.29.6 Uαβ=bαbβYγYδ+YαYβbγbδ,
x=,
Uαβ =XαXβ-bγ-bδ+-bα-bβXγXδ,
y=-.
19.29.7 yxaα+bαtaδ+bδtdts(t)=23dαβdαγRD(Uαβ2,Uαγ2,Uαδ2)+2XαYαXδYδUαδ,
Uαδ0.
19.29.8 yxaα+bαta5+b5tdts(t)=23dαβdαγdαδdα5RJ(U122,U132,U232,Uα52)+2RC(Sα52,Qα52),
Sα52\(-,0),

where

19.29.9 Uα52 =Uαβ2-dαγdαδdβ5dα5=Uβγ2-dαβdαγdδ5dα50,
Sα5 =1x-y(XβXγXδXαY52+YβYγYδYαX52),
Qα5 =X5Y5XαYαUα50,
Sα52-Qα52 =dβ5dγ5dδ5dα5.

The Cauchy principal value is taken when Uα52 or Qα52 is real and negative. Cubic cases of these formulas are obtained by setting one of the factors in (19.29.3) equal to 1.

The advantages of symmetric integrals for tables of integrals and symbolic integration are illustrated by (19.29.4) and its cubic case, which replace the 8+8+12=28 formulas in Gradshteyn and Ryzhik (2000, 3.147, 3.131, 3.152) after taking x2 as the variable of integration in 3.152. Moreover, the requirement that one limit of integration be a branch point of the integrand is eliminated without doubling the number of standard integrals in the result. (19.29.7) subsumes all 72 formulas in Gradshteyn and Ryzhik (2000, 3.168), and its cubic cases similarly replace the 18+36+18=72 formulas in Gradshteyn and Ryzhik (2000, 3.133, 3.142, and 3.141(1-18)). For example, 3.142(2) is included as

19.29.10 uba-t(b-t)(t-c)3dt=-23(a-b)(b-u)3/2RD+2b-c(a-u)(b-u)u-c,
a>b>u>c,

where the arguments of the RD function are, in order, (a-b)(u-c), (b-c)(a-u), (a-b)(b-c).

§19.29(ii) Reduction to Basic Integrals

(19.2.3) can be written

19.29.11 I(m)=yxα=1h(aα+bαt)-1/2j=1n(aj+bjt)mjdt,

where x>y, h=3 or 4, nh, and mj is an integer. Define

19.29.12 m=(m1,,mn)=j=1nmjej,

where ej is an n-tuple with 1 in the jth position and 0’s elsewhere. Define also 0=(0,,0) and retain the notation and conditions associated with (19.29.1) and (19.29.2). The integrals in (19.29.4), (19.29.7), and (19.29.8) are I(0), I(eα-eδ), and I(eα-e5), respectively.

The only cases of I(m) that are integrals of the first kind are the two (h=3 or 4) with m=0. The only cases that are integrals of the third kind are those in which at least one mj with j>h is a negative integer and those in which h=4 and j=1nmj is a positive integer. All other cases are integrals of the second kind.

I(m) can be reduced to a linear combination of basic integrals and algebraic functions. In the cubic case (h=3) the basic integrals are

19.29.13 I(0);
I(-ej),
1jn.

In the quartic case (h=4) the basic integrals are

19.29.14 I(0);
I(-ej),
1jn;
I(eα),
1α4.

Basic integrals of type I(-ej), 1jh, are not linearly independent, nor are those of type I(ej), 1j4.

The reduction of I(m) is carried out by a relation derived from partial fractions and by use of two recurrence relations. These are given in Carlson (1999, (2.19), (3.5), (3.11)) and simplified in Carlson (2002, (1.10), (1.7), (1.8)) by means of modified definitions. Partial fractions provide a reduction to integrals in which m has at most one nonzero component, and these are then reduced to basic integrals by the recurrence relations. A special case of Carlson (1999, (2.19)) is given by

19.29.15 bjI(el-ej)=dljI(-ej)+blI(0),
j,l=1,2,,n,

which shows how to express the basic integral I(-ej) in terms of symmetric integrals by using (19.29.4) and either (19.29.7) or (19.29.8). The first choice gives a formula that includes the 18+9+18 = 45 formulas in Gradshteyn and Ryzhik (2000, 3.133, 3.156, 3.158), and the second choice includes the 8+8+8+12 = 36 formulas in Gradshteyn and Ryzhik (2000, 3.151, 3.149, 3.137, 3.157) (after setting x2=t in some cases).

If h=3, then the recurrence relation (Carlson (1999, (3.5))) has the special case

19.29.16 bβbγI(eα)=dαβdαγI(-eα)+2bα(s(x)aα+bαx-s(y)aα+bαy),

where α,β,γ is any permutation of the numbers 1,2,3, and

19.29.17 s(t)=α=13aα+bαt.

(This shows why I(eα) is not needed as a basic integral in the cubic case.) In the quartic case this recurrence relation has an extra term in I(2eα), and hence I(eα), 1α4, is a basic integral. It can be expressed in terms of symmetric integrals by setting a5=1 and b5=0 in (19.29.8).

The other recurrence relation is

19.29.18 bjqI(qel)=r=0q(qr)blrdljq-rI(rej),
j,l=1,2,,n;

see Carlson (1999, (3.11)). An example that uses (19.29.15)–(19.29.18) is given in §19.34.

For an implementation by James FitzSimons of the method for reducing I(m) to basic integrals and extensive tables of such reductions, see Carlson (1999) and Carlson and FitzSimons (2000).

Another method of reduction is given in Gray (2002). It depends primarily on multivariate recurrence relations that replace one integral by two or more.

§19.29(iii) Examples

The first formula replaces (19.14.4)–(19.14.10). Define Qj(t)=aj+bjt2, j=1,2, and assume both Q’s are positive for 0y<t<x. Then

19.29.19 yxdtQ1(t)Q2(t)=RF(U2+a1b2,U2+a2b1,U2),
19.29.20 yxt2dtQ1(t)Q2(t)=13a1a2RD(U2+a1b2,U2+a2b1,U2)+(xy/U),

and

19.29.21 yxdtt2Q1(t)Q2(t)=13b1b2RD(U2+a1b2,U2+a2b1,U2)+(xyU)-1,

where

19.29.22 (x2-y2)U=xQ1(y)Q2(y)+yQ1(x)Q2(x).

If both square roots in (19.29.22) are 0, then the indeterminacy in the two preceding equations can be removed by using (19.27.8) to evaluate the integral as RG(a1b2,a2b1,0) multiplied either by -2/(b1b2) or by -2/(a1a2) in the cases of (19.29.20) or (19.29.21), respectively. If x=, then U is found by taking the limit. For example,

19.29.23 ydt(t2+a2)(t2-b2)=RF(y2+a2,y2-b2,y2).

Next, for j=1,2, define Qj(t)=fj+gjt+hjt2, and assume both Q’s are positive for y<t<x. If each has real zeros, then (19.29.4) may be simpler than

19.29.24 yxdtQ1(t)Q2(t)=4RF(U,U+D12+V,U+D12-V),

where

19.29.25 (x-y)2U =S1S2,
Sj =(Qj(x)+Qj(y))2-hj(x-y)2,
Djl =2fjhl+2hjfl-gjgl,
V =D122-D11D22.

(The variables of RF are real and nonnegative unless both Q’s have real zeros and those of Q1 interlace those of Q2.) If Q1(t)=(a1+b1t)(a2+b2t), where both linear factors are positive for y<t<x, and Q2(t)=f2+g2t+h2t2, then (19.29.25) is modified so that

19.29.26 S1 =(X1Y2+Y1X2)2,
Xj =aj+bjx,
Yj =aj+bjy,
D12 =2a1a2h2+2b1b2f2-(a1b2+a2b1)g2,
D11 =-(a1b2-a2b1)2=-d122,

with other quantities remaining as in (19.29.25). In the cubic case, in which a2=1, b2=0, (19.29.26) reduces further to

19.29.27 S1 =(X1+Y1)2,
D12 =2a1h2-b1g2,
D11 =-b12.

For example, because t3-a3=(t-a)(t2+at+a2), we find that when 0ay<x

19.29.28 yxdtt3-a3=4RF(U,U-3a+23a,U-3a-23a),

where

19.29.29 (x-y)2U =(x-a+y-a)2((ξ+η)2-(x-y)2),
ξ =x2+ax+a2,
η =y2+ay+a2.

Lastly, define Q(t2)=f+gt2+ht4 and assume Q(t2) is positive and monotonic for y<t<x. Then

19.29.30 yxdtQ(t2)=2RF(U,U-g+2fh,U-g-2fh),

where

19.29.31 (x-y)2U=(Q(x2)+Q(y2))2-h(x2-y2)2.

For example, if 0yx and a40, then

19.29.32 yxdtt4+a4=2RF(U,U+2a2,U-2a2),

where

19.29.33 (x-y)2U=(x4+a4+y4+a4)2-(x2-y2)2.