About the Project
19 Elliptic IntegralsLegendre’s Integrals

§19.8 Quadratic Transformations

Contents
  1. §19.8(i) Gauss’s Arithmetic-Geometric Mean (AGM)
  2. §19.8(ii) Landen Transformations
  3. §19.8(iii) Gauss Transformation

§19.8(i) Gauss’s Arithmetic-Geometric Mean (AGM)

When a_{0} and g_{0} are positive numbers, define

19.8.1
a_{n+1}=\frac{a_{n}+g_{n}}{2},
g_{n+1}=\sqrt{a_{n}g_{n}},n=0,1,2,\dots.

As n\to\infty, a_{n} and g_{n} converge to a common limit M\left(a_{0},g_{0}\right) called the AGM (Arithmetic-Geometric Mean) of a_{0} and g_{0}. By symmetry in a_{0} and g_{0} we may assume a_{0}\geq g_{0} and define

19.8.2 c_{n}=\sqrt{a_{n}^{2}-g_{n}^{2}}.

Then

19.8.3 c_{n+1}=\frac{a_{n}-g_{n}}{2}=\frac{c_{n}^{2}}{4a_{n+1}},

showing that the convergence of c_{n} to 0 and of a_{n} and g_{n} to M\left(a_{0},g_{0}\right) is quadratic in each case.

The AGM appears in

and in

where a_{0}=1, g_{0}=k^{\prime}, p_{0}^{2}=1-\alpha^{2}, Q_{0}=1, and

19.8.8
p_{n+1}=\frac{p_{n}^{2}+a_{n}g_{n}}{2p_{n}},
\varepsilon_{n}=\frac{p_{n}^{2}-a_{n}g_{n}}{p_{n}^{2}+a_{n}g_{n}},
Q_{n+1}=\tfrac{1}{2}Q_{n}\varepsilon_{n},n=0,1,\dots.

Again, p_{n} and \varepsilon_{n} converge quadratically to M\left(a_{0},g_{0}\right) and 0, respectively, and Q_{n} converges to 0 faster than quadratically. If \alpha^{2}>1, then the Cauchy principal value is

where (19.8.8) still applies, but with

19.8.10 p_{0}^{2}=1-(k^{2}/\alpha^{2}).

§19.8(ii) Landen Transformations

Descending Landen Transformation

Let

19.8.11
k_{1}=\frac{1-k^{\prime}}{1+k^{\prime}},
\phi_{1}=\phi+\operatorname{arctan}\left(k^{\prime}\tan\phi\right)=%
\operatorname{arcsin}\left((1+k^{\prime})\frac{\sin\phi\cos\phi}{\sqrt{1-k^{2}%
{\sin}^{2}\phi}}\right).

(Note that 0<k<1 and 0<\phi<\pi/2 imply k_{1}<k and \phi<\phi_{1}<2\phi, and also that \phi=\pi/2 implies \phi_{1}=\pi.) Then

where

19.8.15
\omega^{2}=\frac{k^{2}-\alpha^{2}}{1-\alpha^{2}},
\alpha_{1}^{2}=\frac{\alpha^{2}\omega^{2}}{(1+k^{\prime})^{2}},
c_{1}={\csc}^{2}\phi_{1}.

Ascending Landen Transformation

Let

19.8.16
k_{2}=2\sqrt{k}/(1+k),
2\phi_{2}=\phi+\operatorname{arcsin}\left(k\sin\phi\right).

(Note that 0<k<1 and 0<\phi\leq\pi/2 imply k<k_{2}<1 and \phi_{2}<\phi.) Then

§19.8(iii) Gauss Transformation

We consider only the descending Gauss transformation because its (ascending) inverse moves F\left(\phi,k\right) closer to the singularity at k=\sin\phi=1. Let

19.8.18
k_{1}=(1-k^{\prime})/(1+k^{\prime}),
\sin\psi_{1}=\frac{(1+k^{\prime})\sin\phi}{1+\Delta},
\Delta=\sqrt{1-k^{2}{\sin}^{2}\phi}.

(Note that 0<k<1 and 0<\phi<\pi/2 imply k_{1}<k and \psi_{1}<\phi, and also that \phi=\pi/2 implies \psi_{1}=\pi/2, thus preserving completeness.) Then

where

19.8.21
\rho=\sqrt{1-(k^{2}/\alpha^{2})},
\alpha_{1}^{2}=\alpha^{2}(1+\rho)^{2}/(1+k^{\prime})^{2},
c={\csc}^{2}\phi.

If 0<\alpha^{2}<k^{2}, then \rho is pure imaginary.