■円周率・レムニスケート周率・算術幾何平均(その3)

(1)2数a0,b0をとり,a1=(a0+b0)/2,b1=√a0b0を計算する.次に,a2=(a1+b1)/2,b2=√a1b1とする.

 0<a<bのとき

 a=a0<a1<a2<・・・>b2>b1>b0>b

M(a,b)=liman=limbn

M(ra,rb)=rM(a,b)

M(a,b)=bM(a/b,1)

M(a,b)=aM(1,b/a)

k’=b/a,k^2=1−k’^2

  M(a,b)=1/(2/π∫(0,1)dx/√(1−x^2)(1−k^2x^2))

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

(2)a0=1,b0=√2から始めると、その極限値は、π/ω=1.19・・・

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

驚くべきことに算術幾何平均は楕円積分にも現れます。

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

 f=G,g=Aならば極限M(a,b)は楕円積分

  M(a,b)=1/(2/π∫(0,π/2)dφ/√{(acosφ)^2+(bsinφ)^2})

により表すことができる(ガウス).

∫(0,π/2)dφ/√{(acosφ)^2+(bsinφ)^2}=∫(0,π/2)dφ/√{((a+b)/2・cosφ)^2+(√ab・/sinφ)^2}

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

算術幾何平均の一般化

(a,b)→(f(a,b),g(a,b))

f(a,b)=(a+b)/2,g(a,b)=(ab)^1/2ですが、

f(a,b)=(a+2b)/3,g(a,b)={(a^2+ab+b^2)b/3}^1/3

f(a,b)={(a^2+2b^2)/4}^1/2,g(a,b)={(a^2+b^2)b^2/2}^1/4

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