[go: up one dir, main page]

19 Elliptic IntegralsSymmetric Integrals

§19.22 Quadratic Transformations

Contents
  1. §19.22(i) Complete Integrals
  2. §19.22(ii) Gauss’s Arithmetic-Geometric Mean (AGM)
  3. §19.22(iii) Incomplete Integrals

§19.22(i) Complete Integrals

Let ℜ⁡x>0, ℜ⁡y>0, a=(x+y)/2, and p≠0. Then

19.22.1 RF⁡(0,x2,y2) =RF⁡(0,x⁢y,a2),
19.22.2 2⁢RG⁡(0,x2,y2) =4⁢RG⁡(0,x⁢y,a2)−x⁢y⁢RF⁡(0,x⁢y,a2),
19.22.3 2⁢y2⁢RD⁡(0,x2,y2) =14⁢(y2−x2)⁢RD⁡(0,x⁢y,a2)+3⁢RF⁡(0,x⁢y,a2).
19.22.4 (p±2−p∓2)⁢RJ⁡(0,x2,y2,p2)=2⁢(p±2−a2)⁢RJ⁡(0,x⁢y,a2,p±2)−3⁢RF⁡(0,x⁢y,a2)+3⁢π/(2⁢p),

where

19.22.5 2⁢p±=(p+x)⁢(p+y)±(p−x)⁢(p−y),

and hence

19.22.6 p+⁢p− =p⁢a,
p+2+p−2 =p2+x⁢y,
p+2−p−2 =(p2−x2)⁢(p2−y2),
4⁢(p±2−a2) =(p2−x2±p2−y2)2.

Bartky’s Transformation

19.22.7 2⁢p2⁢RJ⁡(0,x2,y2,p2)=v+⁢v−⁢RJ⁡(0,x⁢y,a2,v+2)+3⁢RF⁡(0,x⁢y,a2),
v±=(p2±x⁢y)/(2⁢p).

If p=y, then (19.22.7) reduces to (19.22.3), but if p=x or p=y, then both sides of (19.22.4) are 0 by (19.20.9). If x<p<y or y<p<x, then p+ and p− are complex conjugates.

§19.22(ii) Gauss’s Arithmetic-Geometric Mean (AGM)

The AGM, M⁡(a0,g0), of two positive numbers a0 and g0 is defined in §19.8(i). Again, we assume that a0≥g0 (except in (19.22.10)), and define cn=an2−gn2. Then

19.22.9 4π⁢RG⁡(0,a02,g02)=1M⁡(a0,g0)⁢(a02−∑n=0∞2n−1⁢cn2)=1M⁡(a0,g0)⁢(a12−∑n=2∞2n−1⁢cn2),

and

where

19.22.11 Q0 =1,
Qn+1 =12⁢Qn⁢an−gnan+gn.

Qn has the same sign as a0−g0 for n≥1.

where p0>0 and

19.22.13 pn+1 =pn2+an⁢gn2⁢pn,
εn =pn2−an⁢gnpn2+an⁢gn,
Q0 =1,
Qn+1 =12⁢Qn⁢εn.

(If p0=a0, then pn=an and (19.22.13) reduces to (19.22.11).) As n→∞, pn and εn converge quadratically to M⁡(a0,g0) and 0, respectively, and Qn converges to 0 faster than quadratically. If the last variable of RJ is negative, then the Cauchy principal value is

19.22.14 RJ⁡(0,g02,a02,−q02)=−3⁢π4⁢M⁡(a0,g0)⁢(q02+a02)⁢(2+a02−g02q02+g02⁢∑n=0∞Qn),

and (19.22.13) still applies, provided that

19.22.15 p02=a02⁢(q02+g02)/(q02+a02).

§19.22(iii) Incomplete Integrals

Let x, y, and z have positive real parts, assume p≠0, and retain (19.22.5) and (19.22.6). Define

19.22.16 a =(x+y)/2,
2⁢z± =(z+x)⁢(z+y)±(z−x)⁢(z−y),

so that

19.22.17 z+⁢z− =z⁢a,
z+2+z−2 =z2+x⁢y,
z+2−z−2 =(z2−x2)⁢(z2−y2),
4⁢(z±2−a2) =(z2−x2±z2−y2)2.

Then

19.22.18 RF⁡(x2,y2,z2)=RF⁡(a2,z−2,z+2),
19.22.19 (z±2−z∓2)⁢RD⁡(x2,y2,z2)=2⁢(z±2−a2)⁢RD⁡(a2,z∓2,z±2)−3⁢RF⁡(x2,y2,z2)+(3/z),
19.22.20 (p±2−p∓2)⁢RJ⁡(x2,y2,z2,p2)=2⁢(p±2−a2)⁢RJ⁡(a2,z+2,z−2,p±2)−3⁢RF⁡(x2,y2,z2)+3⁢RC⁡(z2,p2),
19.22.21 2⁢RG⁡(x2,y2,z2)=4⁢RG⁡(a2,z+2,z−2)−x⁢y⁢RF⁡(x2,y2,z2)−z,
19.22.22 RC⁡(x2,y2)=RC⁡(a2,a⁢y).

If x,y,z are real and positive, then (19.22.18)–(19.22.21) are ascending Landen transformations when x,y<z (implying a<z−<z+), and descending Gauss transformations when z<x,y (implying z+<z−<a). Ascent and descent correspond respectively to increase and decrease of k in Legendre’s notation. Descending Gauss transformations include, as special cases, transformations of complete integrals into complete integrals; ascending Landen transformations do not.

If p=x or p=y, then (19.22.20) reduces to 0=0 by (19.20.13), and if z=x or z=y then (19.22.19) reduces to 0=0 by (19.20.20) and (19.22.22). If x<z<y or y<z<x, then z+ and z− are complex conjugates. However, if x and y are complex conjugates and z and p are real, then the right-hand sides of all transformations in §§19.22(i) and 19.22(iii)—except (19.22.3) and (19.22.22)—are free of complex numbers and p±2−p∓2=±|p2−x2|≠0.

The transformations inverse to the ones just described are the descending Landen transformations and the ascending Gauss transformations. The equations inverse to (19.22.5) and (19.22.16) are given by

19.22.23 x+y =2⁢a,
x−y =(2/a)⁢(a2−z+2)⁢(a2−z−2),
z =z+⁢z−/a,

and the corresponding equations with z, z+, and z− replaced by p, p+, and p−, respectively. These relations need to be used with caution because y is negative when 0<a<z+⁢z−⁢(z+2+z−2)−1/2.