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)
=π
となるというわけである.
===================================