■円周率の計算(その12)
Stirlingの公式の大体の形は(n!/n^n+1/2・e(−n)がn→∞のときに一定の極限値をもつ)はいろいろな方法で導かれますが,その極限値を正しく√(2π)とするにはWallisの公式が不可欠のようです.実際,これまでの証明はすべてWallisの公式(ないしそれと同等の結果)を活用しております.
===================================
【1】ウォリスの公式
1/2B(1/2,(n+1)/2)=integral(0-π/2)(sinθ)^ndθ
この値をSnとおくと,部分積分により漸化式
Sn=(n-1)/nSn-2
が得られますから,
n=2k(偶数)なら1・3・・・(2k-1)/2・4・・・(2k)*π/2
n=2k+1(奇数)なら2・4・・・(2k)/1・3・・・(2k+1)
これより,
lim1・3・・・(2k-1)/2・4・・・(2k)*√(k)=1/√(π)
変形するとウォリスの公式
(2n)!/(2^nn!)^2√(n)=1/√(π)
が得られる.
===================================
【2】ウォリスの公式(その2)
ウォリスの公式は
π/2=(2・2/1・3)(4・4)/(3・5)(6・6)/(5・7)・・・(2n・2n)/((2n−1)・(2n+1))・・・
と記すとわかりやすい.この公式の不思議なところは有理数の無限積→πになっている点である.
===================================
【3】ヴィエトの公式
まず,任意のxに対して,無限積公式
sinx/x=cosx/2cosx/4cosx/8・・・
も示しておこう.
(証明)
sinx=2sinx/2cosx/2
=4sinx/4cosx/4cosx/2
=8sinx/8cosx/8cosx/4cosx/2
・・・・・
=2^nsinx/2^ncosx/2^n・・・cosx/2
書き直すと
sinx=x[sin(x/2^n)/(x/2^n)]cosx/2・・・cosx/2^n
ここで,n→∞のとき,
sin(x/2^n)/(x/2^n)→1
であるから,sinxの無限積表示
sinx=xΠcosx/2^n
=x(1−x^2/π^2)(1−x^2/4π^2)(1−x^2/9π^2)・・・
が得られる.この結果は,sinxがx=0,±π,±2π,±3π,・・・のとき,0になることに一致している.
sinx/x=cosx/2cosx/4cosx/8・・・
において,x=π/2とおくと,
2/π=cosπ/2^2cosπ/2^3cosπ/2^4・・・
===================================
【4】ヴィエトの公式(その2)
半角の公式
cosx/2=√(1+cosx)/2
また,
cosπ/2^2=1/√2
cosπ/2^3=√(1+cosπ/2^2)/2=√(1/2+1/2√1/2)
cosπ/2^4=√(1+cosπ/2^3)/2=√(1/2+1/2(1/2+1/2√1/2))
以下同様の計算を続けることによって,
2/π=√(1/2)√(1/2+1/2√1/2)√(1/2+1/2(1/2+1/2√1/2))・・・
が得られる.
===================================