■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/√(π)

が得られる.

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