■円周率の計算(その35)
[2]パーシバル,ベラード(1997年)
π=1/2^6Σ(−1)^n/2^10n(2^5/(4n+1)−1/(4n+3)+2^8/(10n+1)−2^6/(10n+3)−2^2/(10n+5)−2^2/(10n+7)+1/(10n+9))
の導出を試みたい.
===================================
0<t<1のとき,
1/(1+t^20)=Σt^20n=1−t^20+t^40−・・・
t^k-1/(1+t^20)=t^k-1Σt^20n=t^k-1(1−t^20+t^40−・・・)
また, 上式を0から1/√2まで積分する
Ik=∫(0,1/√2)t^k-1/(1+t^20)dt
=∫(0,1/√2)Σ(−1)^nt^k+20n-1dt
=Σ∫(0,1/√2)(−1)^nt^k+20n-1dt
=Σ(−1)^nt^k+20n/(k+20n)
=1/2^k/2Σ(−1)^n/2^10n・1/(20n+k)
S(a,b)=Σ(−1)^n/2^10n・1/(an+b)
とおくと,この公式は
1/2^6Σ(−1)^n/2^10n(2^5/(4n+1)−1/(4n+3)+2^8/(10n+1)−2^6/(10n+3)−2^2/(10n+5)−2^2/(10n+7)+1/(10n+9))
=2^-15S(20,5)−2^-65S(20,15)+2^3S(20,2)−2S(20,6)−2^-3S(20,10)−2^-3S(20,14)+2^-5S(20,18)
=∫(0,1/√2)(10√2t^4−10√2t^14+16t−16t^5−4t^9−16t^13+16t^17)/(1+t^20)dt
この不定積分も結構面倒な形になる.次回の宿題としたい.
===================================