Digital Library of Mathematical Functions
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,\infty) containing the square root of a cubic polynomial (compare §19.16(i)). Let

19.29.1
X_{\alpha}=\sqrt{a_{\alpha}+b_{\alpha}x},
Y_{\alpha}=\sqrt{a_{\alpha}+b_{\alpha}y},x>y, 1\leq\alpha\leq 5,
19.29.2d_{{\alpha\beta}}=a_{\alpha}b_{\beta}-a_{\beta}b_{\alpha},d_{{\alpha\beta}}\neq 0 if \alpha\neq\beta,

and assume that the line segment with endpoints a_{\alpha}+b_{\alpha}x and a_{\alpha}+b_{\alpha}y lies in \Complex\setminus(-\infty,0) for 1\leq\alpha\leq 4. If

19.29.3s(t)=\prod_{{\alpha=1}}^{4}\sqrt{a_{\alpha}+b_{\alpha}t}

and \alpha,\beta,\gamma,\delta is any permutation of the numbers 1,2,3,4, then

where

19.29.5
U_{{\alpha\beta}}=(X_{\alpha}X_{\beta}Y_{\gamma}Y_{\delta}+Y_{\alpha}Y_{\beta}%
X_{\gamma}X_{\delta})/(x-y),
U_{{\alpha\beta}}=U_{{\beta\alpha}}=U_{{\gamma\delta}}=U_{{\delta\gamma}},
U_{{\alpha\beta}}^{2}-U_{{\alpha\gamma}}^{2}=d_{{\alpha\delta}}d_{{\beta\gamma%
}}.

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

19.29.6
{U_{{\alpha\beta}}=\sqrt{b_{\alpha}}\sqrt{b_{\beta}}Y_{\gamma}Y_{\delta}+Y_{%
\alpha}Y_{\beta}\sqrt{b_{\gamma}}\sqrt{b_{\delta}},}x=\infty,
U_{{\alpha\beta}}=X_{\alpha}X_{\beta}\sqrt{-b_{\gamma}}\sqrt{-b_{\delta}}+%
\sqrt{-b_{\alpha}}\sqrt{-b_{\beta}}X_{\gamma}X_{\delta},y=-\infty.

where

19.29.9
U_{{\alpha 5}}^{2}=U_{{\alpha\beta}}^{2}-\frac{d_{{\alpha\gamma}}d_{{\alpha%
\delta}}d_{{\beta 5}}}{d_{{\alpha 5}}}=U_{{\beta\gamma}}^{2}-\frac{d_{{\alpha%
\beta}}d_{{\alpha\gamma}}d_{{\delta 5}}}{d_{{\alpha 5}}}\neq 0,
S_{{\alpha 5}}=\frac{1}{x-y}\left(\frac{X_{\beta}X_{\gamma}X_{\delta}}{X_{%
\alpha}}Y_{5}^{2}+\frac{Y_{\beta}Y_{\gamma}Y_{\delta}}{Y_{\alpha}}X_{5}^{2}%
\right),
Q_{{\alpha 5}}=\frac{X_{5}Y_{5}}{X_{\alpha}Y_{\alpha}}U_{{\alpha 5}}\neq 0,
S_{{\alpha 5}}^{2}-Q_{{\alpha 5}}^{2}=\frac{d_{{\beta 5}}d_{{\gamma 5}}d_{{%
\delta 5}}}{d_{{\alpha 5}}}.

The Cauchy principal value is taken when U_{{\alpha 5}}^{2} or Q_{{\alpha 5}}^{2} 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 x^{2} 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

where the arguments of the \mathop{R_{D}\/}\nolimits 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

where x>y, h=3 or 4, n\geq h, and m_{j} is an integer. Define

where \mathbf{e}_{j} is an n-tuple with 1 in the jth position and 0’s elsewhere. Define also \boldsymbol{{0}}=(0,\dots,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(\boldsymbol{{0}}), I(\mathbf{e}_{\alpha}-\mathbf{e}_{\delta}), and I(\mathbf{e}_{\alpha}-\mathbf{e}_{5}), respectively.

The only cases of I(\mathbf{m}) that are integrals of the first kind are the two (h=3 or 4) with \mathbf{m}=\boldsymbol{{0}}. The only cases that are integrals of the third kind are those in which at least one m_{j} with j>h is a negative integer and those in which h=4 and \sum_{{j=1}}^{n}m_{j} is a positive integer. All other cases are integrals of the second kind.

I(\mathbf{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(\boldsymbol{{0}});
I(-\mathbf{e}_{j}),1\leq j\leq n.

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

19.29.14
I(\boldsymbol{{0}});
I(-\mathbf{e}_{j}),1\leq j\leq n;
I(\mathbf{e}_{\alpha}),1\leq\alpha\leq 4.

Basic integrals of type I(-\mathbf{e}_{j}), 1\leq j\leq h, are not linearly independent, nor are those of type I(\mathbf{e}_{j}), 1\leq j\leq 4.

The reduction of I(\mathbf{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 \mathbf{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

which shows how to express the basic integral I(-\mathbf{e}_{j}) 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 x^{2}=t in some cases).

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

19.29.16b_{\beta}b_{\gamma}I(\mathbf{e}_{\alpha})=d_{{\alpha\beta}}d_{{\alpha\gamma}}I%
(-\mathbf{e}_{\alpha})+2b_{\alpha}\left(\frac{s(x)}{a_{\alpha}+b_{\alpha}x}-%
\frac{s(y)}{a_{\alpha}+b_{\alpha}y}\right),

where \alpha,\beta,\gamma is any permutation of the numbers 1,2,3, and

19.29.17s(t)=\prod_{{\alpha=1}}^{3}\sqrt{a_{\alpha}+b_{\alpha}t}.

(This shows why I(\mathbf{e}_{\alpha}) is not needed as a basic integral in the cubic case.) In the quartic case this recurrence relation has an extra term in I(2\mathbf{e}_{\alpha}), and hence I(\mathbf{e}_{\alpha}), 1\leq\alpha\leq 4, is a basic integral. It can be expressed in terms of symmetric integrals by setting a_{5}=1 and b_{5}=0 in (19.29.8).

The other recurrence relation is

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(\mathbf{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 Q_{j}(t)=a_{j}+b_{j}t^{2}, j=1,2, and assume both Q’s are positive for 0\leq y<t<x. Then

and

where

19.29.22(x^{2}-y^{2})U=x\sqrt{Q_{1}(y)Q_{2}(y)}+y\sqrt{Q_{1}(x)Q_{2}(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 \mathop{R_{G}\/}\nolimits\!\left(a_{1}b_{2},a_{2}b_{1},0\right) multiplied either by -2/(b_{1}b_{2}) or by -2/(a_{1}a_{2}) in the cases of (19.29.20) or (19.29.21), respectively. If x=\infty, then U is found by taking the limit. For example,

Next, for j=1,2, define Q_{j}(t)=f_{j}+g_{j}t+h_{j}t^{2}, and assume both Q’s are positive for y<t<x. If each has real zeros, then (19.29.4) may be simpler than

where

19.29.25
(x-y)^{2}U=S_{1}S_{2},
S_{j}=\left(\sqrt{Q_{j}(x)}+\sqrt{Q_{j}(y)}\right)^{2}-h_{j}(x-y)^{2},
D_{{jl}}=2f_{j}h_{l}+2h_{j}f_{l}-g_{j}g_{l},
V=\sqrt{D_{{12}}^{2}-D_{{11}}D_{{22}}}.

(The variables of \mathop{R_{F}\/}\nolimits are real and nonnegative unless both Q’s have real zeros and those of Q_{1} interlace those of Q_{2}.) If Q_{1}(t)=(a_{1}+b_{1}t)(a_{2}+b_{2}t), where both linear factors are positive for y<t<x, and Q_{2}(t)=f_{2}+g_{2}t+h_{2}t^{2}, then (19.29.25) is modified so that

19.29.26
S_{1}=(X_{1}Y_{2}+Y_{1}X_{2})^{2},
X_{j}=\sqrt{a_{j}+b_{j}x},
Y_{j}=\sqrt{a_{j}+b_{j}y},
D_{{12}}=2a_{1}a_{2}h_{2}+2b_{1}b_{2}f_{2}-(a_{1}b_{2}+a_{2}b_{1})g_{2},
D_{{11}}=-(a_{1}b_{2}-a_{2}b_{1})^{2}=-d_{{12}}^{2},

with other quantities remaining as in (19.29.25). In the cubic case, in which a_{2}=1, b_{2}=0, (19.29.26) reduces further to

19.29.27
S_{1}=(X_{1}+Y_{1})^{2},
D_{{12}}=2a_{1}h_{2}-b_{1}g_{2},
D_{{11}}=-b_{1}^{2}.

For example, because t^{3}-a^{3}=(t-a)(t^{2}+at+a^{2}), we find that when 0\leq a\leq y<x

where

19.29.29
(x-y)^{2}U=(\sqrt{x-a}+\sqrt{y-a})^{2}\left((\xi+\eta)^{2}-(x-y)^{2}\right),
\xi=\sqrt{x^{2}+ax+a^{2}},
\eta=\sqrt{y^{2}+ay+a^{2}}.

Lastly, define Q(t^{2})=f+gt^{2}+ht^{4} and assume Q(t^{2}) is positive and monotonic for y<t<x. Then

where

19.29.31(x-y)^{2}U=\left(\sqrt{Q(x^{2})}+\sqrt{Q(y^{2})}\right)^{2}-h(x^{2}-y^{2})^{2}.

For example, if 0\leq y\leq x and a^{4}\geq 0, then

where

19.29.33(x-y)^{2}U=\left(\sqrt{x^{4}+a^{4}}+\sqrt{y^{4}+a^{4}}\right)^{2}-(x^{2}-y^{2}%
)^{2}.