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

=π

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

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