■1000!は10^250より大きいか? (その26)
10!で300万を超えるから、100!、1000!はどれくらい巨大な数になるのだろうか?
===================================
1/2+7/(8logn)<(log(n!)-nlogn+n)/logn<1/2+1/logn
7/8<log{n!/n^1/2・(n/e)^n<1
exp(7/8)/sqr(2π)<n!/sqr(2nπ)・(n/e)^n<exp(1)//sqr(2π)
.957<n!/sqr(2nπ)・(n/e)^n<1.084
===================================
スターリングの公式
n→∞のとき、n!/sqr(2nπ)・(n/e)^n→1
===================================
【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/√(π)
が得られる.
===================================