■円周率の計算(その39)
1976年,arctan(x)の展開公式よりも格段に優れた新しい公式が発表されました.サラミンとブレントは独立に楕円積分の計算と関係したガウス・ルジャンドルの算術幾何平均法という強力な武器を提唱しました.
===================================
【1】算術幾何平均
2数a0 ,b0 をとり,それらの算術平均a1 =(a0 +b0 )/2,幾何平均b1 =√a0 b0 を計算する.次に,a1 ,b1 の平均を計算し,a2 =(a1 +b1 )/2,b2 =√a1 b1 とする.すると,an とbn は急速に同じ極限に到達する.
b0≦b1≦・・・≦bn≦an≦・・・≦a1≦a0
これを算術幾何平均とよぶ.
cn=√(an^2−bn^2)
と定義すると
c1=(a0−b0)/2≦√(a^2−b^2)=c0/2
・・・・・・・・・・・・・・・・・・・・・・・・
0≦cn≦c0/2^n
===================================
【2】ガウスの公式
π=2{AGM(1,1/√2)}^2/(1−Σ2^ncn^2)
AGM(1,1/√2)=0.8472130847
より,ガウスは
AGM(1,1/√2)=π/2・1/K(1/√2)
一般に
AGM(1,m)=π/2・1/K(√1−m^2)
であることを見抜き,その証明を与えた.
1/K(√k)=2/π・AGM(1,√1−k^2)
E(√k)/K(√k)=1−Σ2^n-1cn^2
===================================