■円周率の計算(その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))・・・

が得られる.

===================================