■円周率の計算(その33)

 BBP公式

[1]π=Σ1/16^n(4/(8n+1)-2/(8n+4)-1/(8n+5)-1/(8n+6))

の導出にならって,AW公式

[2]アダムチック,ワゴン(1997年)

  π=Σ(-1)^n/4^n(2/(4n+1)+2/(4n+2)+1/(4n+3))

の導出を試みたい.

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

 0<t<1のとき,

  1/(1+t^4)=Σt^4=1-t^4+t^8+・・・

  t^k-1/(1+t^4)=t^k-1Σt^4=t^k-1(1-t^4+t^8+・・・)

 また, 上式を0から1/√2まで積分する

  Ik=∫(0,1/√2)(-1)^nt^k-1/(1+t^4)dt

=∫(0,1/√2)Σ(-1)^nt^k+4n-1dt

=Σ∫(0,1/√2)(-1)^nt^k+4n-1dt

=Σ(-1)^nt^k+4n/(k+4n)

=1/2^k/2Σ((-1)^n/4^n・1/(4n+k))

  S(a,b)=Σ((-1)^n/4^n・1/(an+b))

とおいて,S(a,b)=πとなる組み合わせを考える.

 AW公式は

  2S(4,1)+2S(4,2)+S(4,3)

=Σ1/4^n{2/(4n+1)+2/(4n+2)+1/(4n+3)}

=∫(0,1/√2)(2√2+4t+2√2t^2)/(1+t^4)dt

=2√2∫(0,1/√2)(1+√2t+t^2)/(1+t^4)dt

=2√2∫(0,1/√2)1/(1-√2t+t^2)dt

=2√2∫(0,1/√2)1/(t-1/√2)^2+1/2)dt

=4√2∫(0,1/√2)1/(√2t-1)^2+1)dt

=4∫(0,1)1/(x-1)^2+1)dx

=4arctan(x-1)|(0,1)

=π

となるというわけである.

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