■円周率の計算(その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)
=π
となるというわけである.
===================================