■シンク関数の2分割(その19)
【2】スターリングの近似公式別証
(Q)任意の自然数nに対して
n!=√(2πn)n^nexp(−n)exp(θ/12n)
を満たす0≦θ≦1が存在することを証明せよ.
(A)f(x)=(1/2)log(1+x)/(1−x)−x
とおく(0≦x≦1).
f(0)=0
f’(x)=x^2/(1−x^2)>0
さらに,
g(x)=f(x)−x^3/3(1−x^2)
とおけば,
g(0)=0
g’(x)=−2x^4/3(1−x^2)^2<0
よって,
0≦(1/2)log(1+x)/(1−x)−x≦x^3/3(1−x^2)
が成り立つ.x=1/(2n+1)を代入すると
0≦(n+1/2)log(n+1)/n−1≦1/12(1/n−1/(n+1))
an=n^(n+1/2)exp(-n)/n!,bn=anexp(1/12n),an<bnとおくと
logan+1/an≧0よりan+1≧an
logbn+1/bn≦0よりbn+1≦bn
よって
cn^(n+1/2)exp(−n)≦n!≦cn^(n+1/2)exp(−n)exp(1/12n)
===================================
【3】ウォリスの公式
さらに,ウォリスの公式より,c=√2πが示される.
(証) 1/2B(1/2,(n+1)/2)=∫(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/√(π)
が得られる.
===================================