■円周率の計算(その34)
[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))
の前に
[3]π^2=9/8Σ1/64^n(16/(6n+1)+8/(6n+2)−2/(6n+4)−1/(6n+5))
の導出を試みたい.定積分がπ^2になる原始関数は何だろうか?
===================================
0<t<1のとき,
1/(1−t^12)=Σt^12n=1+t^12+t^24+・・・
t^k-1/(1−t^12)=t^k-1Σt^12n=t^k-1(1+t^12+t^24+・・・)
また, 上式を0から1/√2まで積分する
Ik=∫(0,1/√2)t^k-1/(1−t^12)dt
=∫(0,1/√2)Σt^k+12n-1dt
=Σ∫(0,1/√2)t^k+12n-1dt
=Σt^k+12n/(k+12n)
=1/2^k/2Σ(1/64^n・1/(12n+k))
S(a,b)=Σ(1/64^n・1/(an+b))
とおくと,この公式は
Σ1/64^n{16/(6n+1)+8/(6n+2)−2/(6n+4)−1/(6n+5)}
=32S(12,2)+16S(12,4)−4S(12,8)−2S(12,10)
=∫(0,1/√2)(64t+64t^3−64t^7−64t^9)/(1−t^12)dt
=64∫(0,1/√2)(t+t^3−t^7−t^9)/(1−t^12)dt
=64∫(0,1/√2)(t+t^3)/(1+t^6)dt
の定数倍である.
この不定積分は公式集を参照することにしたのであるが,結構面倒な形になる.次回の宿題としたい.
===================================