Editing Gauss–Legendre algorithm (section)

== Mathematical background ==

=== Limits of the arithmetic–geometric mean ===

The [[arithmetic–geometric mean]] of two numbers, a<sub>0</sub> and b<sub>0</sub>, is found by calculating the limit of the sequences

:<math>\begin{align} a_{n+1} & = \frac{a_n+b_n}{2}, \\[6pt]
                     b_{n+1} & = \sqrt{a_n b_n},
       \end{align}
</math>

which both converge to the same limit.<br />
If <math>a_0=1</math> and <math>b_0=\cos\varphi</math> then the limit is <math display="inline">{\pi \over 2K(\sin\varphi)}</math> where <math>K(k)</math> is the [[Elliptic integral#Complete elliptic integral of the first kind|complete elliptic integral of the first kind]]

:<math>K(k) = \int_0^{\pi/2} \frac{d\theta}{\sqrt{1-k^2 \sin^2\theta}}.</math>

If <math>c_0 = \sin\varphi</math>, <math>c_{i+1} = a_i - a_{i+1}</math>, then

:<math>\sum_{i=0}^\infty 2^{i-1} c_i^2 = 1 - {E(\sin\varphi)\over K(\sin\varphi)}</math>

where <math>E(k)</math> is the [[Elliptic integral#Complete elliptic integral of the second kind|complete elliptic integral of the second kind]]:

:<math>E(k) = \int_0^{\pi/2}\sqrt {1-k^2 \sin^2\theta}\; d\theta</math> 

Gauss knew of these two results.<ref name="brent">{{Citation
 | last=Brent
 | first=Richard
 | author-link=Richard Brent (scientist)
 | publication-date=
 | year=1975
 | title=Multiple-precision zero-finding methods and the complexity of elementary function evaluation
 | periodical=Analytic Computational Complexity
 | series=
 | publication-place=New York
 | place=
 | publisher=Academic Press
 | editor-last=Traub
 | editor-first=J F
 | volume=
 | issue=
 | pages=151–176
 | url=http://wwwmaths.anu.edu.au/~brent/pub/pub028.html
 | issn=
 | doi=
 | oclc=
 | accessdate=8 September 2007
 | archive-date=23 July 2008
 | archive-url=https://web.archive.org/web/20080723170157/http://wwwmaths.anu.edu.au/~brent/pub/pub028.html
 | url-status=dead
 }}</ref>
<ref name="salamin1">[[Eugene Salamin (mathematician)|Salamin, Eugene]], ''Computation of pi'', Charles Stark Draper Laboratory ISS memo 74–19, 30 January 1974, Cambridge, Massachusetts</ref>
<ref name="salamin2">{{Citation
 | last=Salamin
 | first=Eugene
 | author-link=Eugene Salamin (mathematician)
 | publication-date=
 | year=1976
 | title=Computation of pi Using Arithmetic–Geometric Mean
 | periodical=Mathematics of Computation
 | series=
 | publication-place=
 | place=
 | publisher=
 | editor-last=
 | editor-first=
 | volume=30
 | issue=135
 | pages=565–570
 | url=
 | issn=0025-5718
 | doi=10.2307/2005327
 | jstor=2005327
 | oclc=
 | accessdate=
}}</ref>

=== Legendre’s identity ===
Legendre proved the following identity: 
:<math display="block">K(\cos \theta) E(\sin \theta ) + K(\sin \theta ) E(\cos \theta) - K(\cos \theta) K(\sin \theta) = {\pi \over 2},</math>
for all <math>\theta</math>.<ref name="brent" />

=== Elementary proof with integral calculus ===

The Gauss-Legendre algorithm can be proven to give results converging to π using only integral calculus. This is done here<ref>{{citation|title=Recent Calculations of π: The Gauss-Salamin Algorithm|last1=Lord|first1=Nick|doi=10.2307/3619132|year=1992|journal=The Mathematical Gazette|volume=76|issue=476|pages=231–242|jstor=3619132|s2cid=125865215 }}</ref> and here.<ref>{{citation|title=Easy Proof of Three Recursive π-Algorithms|last1=Milla|first1=Lorenz|arxiv=1907.04110|year=2019}}</ref>