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

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

【１】スターリングの公式の誘導

　スターリングの公式を誘導してみましょう．

　　ｌｏｇｎ！＝ｌｏｇ１＋ｌｏｇ２＋・・・＋ｌｏｇｘ＝Σｌｏｇｘ

ここで，ｙ＝ｌｏｇｘのグラフを幅が１の長方形に分割していくと，ｘが十分大きければ相対的に和の間隔が小さくなるので，和は積分に置き換えられます．

　　Σｌｏｇｘ≒∫(1,x)ｌｏｇｔｄｔ

ｌｏｇｘの原始関数は置換積分よりｘｌｏｇｘ−ｘ＋Ｃと計算されますから，右辺はｘｌｏｇｘ−ｘ＋１となります．したがって，

　　ｎ！≒ｅｎ^nｅｘｐ（−ｎ）

が得られます．

　　ｌｏｇｎ！＝ｎｌｏｇｎ−ｎ＋ｏ（ｎ）

　　　ただし，ｌｉｍｏ（ｎ）／ｎ＝０

としても大体了解されますが，もっと正確に近似すると

　　∫(0,n)ｌｏｇｔｄｔ<ｌｏｇｎ！<∫(1,n+1)ｌｏｇｔｄｔ

より

　　ｎｌｏｇｎ−ｎ<ｌｏｇｎ！<（ｎ＋１）ｌｏｇ（ｎ＋１）−ｎ

したがって，両辺の相加平均に近い（ｎ＋１／２）ｌｏｇｎ−ｎでｌｏｇｎ！を近似できることになり，

　　∫(1,x)ｌｏｇｔｄｔ

　＝ｌｏｇ１＋ｌｏｇ２＋・・・＋ｌｏｇｘ−１／２ｌｏｇｘ＋δ

であること，また，ウォリスの公式：

　　√π〜（ｎ！）^2２^2n／（２ｎ）！√ｎ

より，結局，

　　ｎ！〜√（２πｎ）ｎ^nｅｘｐ（−ｎ）

にたどりつきます．

−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−

　スターリングの近似公式は階乗の一般化であるガンマ関数の近似値としても使われています．

　　Γ（ｘ＋１）＝∫ｅｘｐ（−ｔ）ｔ^xｄｔ〜√（２πｘ）ｘ^xｅｘｐ（−ｘ）

近似の程度を進めると

　　Γ（ｘ＋１）〜√（２πｘ）ｘ^xｅｘｐ{-x[1+1/(12x)+1/(288x^2)-139/(51840x^3)-.....}

が得られます．これらの公式ではｘが大きくなるほど相対誤差は小さくなります．

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝