■円周率の計算(その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

===================================