自然数 n の階乗 (かいじょう、英: factorial )n ! は、1 から n までのすべての自然数の積 である。例えば、
6
!
=
6
×
5
×
4
×
3
×
2
×
1
=
720
{\displaystyle 6! = 6 \times 5 \times 4 \times 3 \times 2 \times 1 = 720 }
である、
空積 の規約のもと 0! = 1 と定義する[1] 。
階乗は数学の様々な場面に出現するが、特に組合せ論 、代数学 、解析学 などが著しい。階乗の最も基本的な出自は n -個の相異なる対象を一列に並べる方法(対象の置換 )の総数が n ! -通りであるという事実である。この事実は少なくとも12世紀にはインドの学者によって知られていた[2] 。ファビアン・ステッドマンは1677年にチェンジリンキング への応用として階乗を記述した[* 1] 。再帰的な手法による記述の後、Stedman は(独自の言葉を用いて)階乗に関しての記述を与えている:
Now the nature of these methods is such, that the changes on one number comprehends [includes] the changes on all lesser numbers, ... insomuch that a compleat Peal of changes on one number seemeth to be formed by uniting of the compleat Peals on all lesser numbers into one entire body;[4]
エクスクラメーションマーク(!)を用いた、この "n ! " という表記は1808年にChristian Krampによって発明された[5] 。
階乗数 オンライン整数列大辞典 の数列 A000142
0!
1
1!
1
2!
2
3!
6
4!
24
5!
120
6!
720
7!
5 040
8!
40 320
9!
362 880
10!
3 628 800
11!
39 916 800
12!
479 001 600
13!
6 227 020 800
14!
87 178 291 200
15!
1 307 674 368 000
16!
20 922 789 888 000
17!
355 687 428 096 000
18!
6 402 373 705 728 000
19!
121 645 100 408 832 000
20!
2 432 902 008 176 640 000
この階乗の定義は、最も重要な性質を残したまま、非整数を引数とする函数に拡張 することができる。そうすれば解析学における著しい手法などの進んだ数学を利用できるようになる。
定義 [ ]
いくつか同値な条件により定義することが可能である。
n
!
=
n
×
(
n
−
1
)
×
⋯
×
3
×
2
×
1
=
∏
k
=
1
n
k
{\displaystyle n!=n\times(n-1)\times\cdots\times3\times2\times1=\prod_{k=1}^n k}
再帰的な定義
n
!
=
{
1
if
n
=
0
,
(
n
−
1
)
!
×
n
if
n
>
0
{\displaystyle n! = \begin{cases}
1 & \text{if } n = 0, \\
(n-1)!\times n & \text{if } n > 0
\end{cases}
}
微分 に関する冪の微分法則を用いた定義[6]
n
!
=
D
n
x
n
(
n
≥
0
)
{\displaystyle n! = D^nx^n \quad(n\geq 0)}
上記の何れの定義においても、
0
!
=
1
{\displaystyle 0! = 1 }
となることが織り込み済みである(最初の定義では「 0 項の積 は 1 と定める」という規約によって)。このように定義することの理由は:
零個の対象の置換は(「何もしない」という)ちょうど一通りであること。 n > 0 のとき有効な漸化式 (n + 1)! = n ! × (n + 1) , が n = 0 の場合にも延長できること。指数函数 などの冪級数としての表示
e
x
=
∑
n
=
0
∞
x
n
n
!
{\displaystyle e^x = \sum_{n=0}^{\infty} \frac{x^n}{n!} }
など多くの公式が短く表せるようになること。組合せ論 における多くの等式が任意のサイズに適用して意味を持つこと。例えば零個の元を空集合 から選ぶ方法の総数は
(
0
0
)
=
0
!
0
!
0
!
=
1
{\displaystyle \binom{0}{0} = \frac{0!}{0!\,0!} = 1}
であり、一般に n -元集合から n -個全ての元を選び出す方法の総数は
(
n
n
)
=
n
!
n
!
0
!
=
1
{\displaystyle \binom nn = \frac{n!}{n!\,0!} = 1}
と書ける。
など様々に挙げることができる。
より進んだ数学においては、引数を非整数として階乗函数を定義することもできる(後述)。そういった一般化された定義のもとでの階乗は関数電卓 や、Maple や Mathematica などの数学ソフトウェア で利用できる。
組合せ論 [ ]
階乗を含む公式は数学の多くの分野に現れるけれども、階乗のおおもとの出自は組合せ論 にある。相異なる n -個の対象の順列 (k -順列)の総数は n ! -通りである。
階乗はしばしば「順番を無視する」という事実を反映するものとして分母 に現れる。古典的な例としては n -個の元から k -個の元を選ぶ組合せ (k -組合せ)の総数が挙げられる。このような組合せは順列から得ることができる。実際、k -順列の総数
n
k
_
=
n
(
n
−
1
)
(
n
−
2
)
⋯
(
n
−
k
+
1
)
{\displaystyle n^{\underline k}=n(n-1)(n-2)\cdots(n-k+1)}
において、順番のみが違う(k -組合せでは違いが無視される)k -順列が k ! -通りずつ存在するから、k -組合せの総数は
n
k
_
k
!
=
n
(
n
−
1
)
(
n
−
2
)
⋯
(
n
−
k
+
1
)
k
(
k
−
1
)
(
k
−
2
)
⋯
1
{\displaystyle \frac{n^{\underline k}}{k!}=\frac{n(n-1)(n-2)\cdots(n-k+1)}{k(k-1)(k-2)\cdots1}}
となる。この数は、二項冪 (1 + X )n における Xk の係数となることから、二項係数
(
n
k
)
{\displaystyle \tbinom nk}
とも呼ばれる。
代数学 に現れる階乗にはいくつも理由があるが、既述の如く二項展開 の係数として現れたり、ある種の演算の対称化 において置換 による平均化を行うなど、組合せ論的な理由で現れるものもある。
微分積分学 においても階乗は例えばテイラー級数 の分母として現れるが、これは冪函数 xn の n -階導函数が n ! であることを補正する定数である。確率論 でも階乗は用いられる。
階乗は数式操作にも有効である。例えば n の k -順列の総数を
n
k
_
=
n
!
(
n
−
k
)
!
{\displaystyle n^{\underline k}=\frac{n!}{(n-k)!}}
と書けば、(この数値を計算することを考えれば効率が悪くなるが)二項係数の対称性
(
n
k
)
=
n
k
_
k
!
=
n
!
(
n
−
k
)
!
k
!
=
n
n
−
k
_
(
n
−
k
)
!
=
(
n
n
−
k
)
{\displaystyle \binom nk=\frac{n^{\underline k}}{k!}=\frac{n!}{(n-k)!k!}=\frac{n^{\underline{n-k}}}{(n-k)!}=\binom n{n-k}}
を見るには都合がよい。
数論における階乗 [ ]
階乗は数論 にも多くの応用を持つ。特に n ! は n 以下の全ての素数 で整除されねばならない。このことの帰結として、n (> 5) が合成数 となる必要十分条件 は
(
n
−
1
)
!
≡
0
(
mod
n
)
{\displaystyle (n-1)!\ \equiv\ 0 \pmod n}
が満たされることである。より強い結果としてウィルソンの定理 は
(
p
−
1
)
!
≡
−
1
(
mod
p
)
{\displaystyle (p-1)!\ \equiv\ -1 \pmod p}
が p が素数であるための必要十分条件であることを述べる。
ルジャンドルの公式 は n ! の素因数分解に現れる p の重複度が
∑
i
=
1
∞
⌊
n
p
i
⌋
{\displaystyle \sum_{i=1}^\infty \left \lfloor \frac n {p^i} \right \rfloor}
であることを示す。これは
n
−
s
p
(
n
)
p
−
1
{\displaystyle \frac{n - s_p(n)}{p - 1}}
と書いてもよい。ただし、s p (n ) は n の p 進展開の係数の和である。
n の階乗自身が素数となる n は 2 のみである。しかし階乗素数 と呼ばれる n ! ± 1 の形の素数は無数に存在する。
1! より大きな階乗は全て偶数 である(これらは明らかに因数 2 を持ち、2 の倍数である)。同様に、5! より後の階乗は 10 の倍数(2 と 5 を因数に持つ)であり、十進展開の末尾に0が並ぶ。
ブロカールの問題 [ ]
ブロカールの問題とは、
n
!
+
1
=
m
2
{\displaystyle n!+1 = m^2}
を満たすn、mは存在するか、という問題である。現在、これを満たす(n , m )の組[* 2] は
(4,5), (5,11), (7,71)
しか見つかっていない。ABC予想 が真であれば、解は有限個しかないことが、Marius Overholtにより示されている。
階乗の解析学 [ ]
ガンマ関数 [ ]
y=Γ(x )
Γ(x +iy )の絶対値 (グラフ中「re」はx に相当、「im」はy に相当)
数学においてガンマ関数 (ガンマかんすう、英: Gamma function )とは、階乗 の概念を一般化した特殊関数 。
互いに同値となるいくつかの定義が存在するが、1729年、数学者オイラー が階乗の一般化として、最初に導入した。
定義 [ ]
実部が正となる複素数
z
{\displaystyle z}
について、次の積分で定義される関数
Γ
(
z
)
=
∫
0
∞
t
z
−
1
e
−
t
d
t
(
ℜ
z
>
0
)
{\displaystyle \Gamma(z)=\int^{\infin}_{0}t^{z-1}e^{-t}\,dt\qquad(\real{z}>0)}
をガンマ関数と呼ぶ[7] 。この積分は、ルシャンドル の定義にしたがって、第二種オイラー積分 とも呼ばれる。元は階乗の一般化としてオイラーが得たもので、
Γ
{\displaystyle \Gamma}
という記号は、ルジャンドルが用いたものである。それ以前は
Π
(
x
)
{\displaystyle \Pi(x)}
などと表記していた(ただし
Π
(
x
)
=
Γ
(
x
+
1
)
{\displaystyle \Pi(x)=\Gamma(x+1)}
)。
一般の複素数
z
{\displaystyle z}
については、解析接続もしくは次の無限乗積で定義される。
Γ
(
z
)
=
lim
n
→
∞
n
z
n
!
∏
k
=
0
n
(
z
+
k
)
{\displaystyle \Gamma(z)=\lim_{n\to\infty}\frac{n^zn!}{\prod_{k=0}^{n}{(z+k)}}}
基本的性質 [ ]
ガンマ関数は、階乗 の複素数 への拡張としてオイラーによって考案されたものであり、自然数
n
{\displaystyle n}
について
Γ
(
n
+
1
)
=
n
!
{\displaystyle \Gamma(n+1)=n! \,}
が成立する。
実際、オイラー積分による定義から
Γ
(
z
)
=
∫
t
=
0
∞
t
z
−
1
(
−
e
−
t
)
′
d
t
=
[
−
t
z
−
1
e
−
t
]
t
=
0
∞
+
(
z
−
1
)
∫
t
=
0
∞
t
z
−
2
e
−
t
d
t
=
(
z
−
1
)
Γ
(
z
−
1
)
{\displaystyle \Gamma(z)=\int_{t=0}^{\infty}{t^{z-1}(-e^{-t})'}\,dt=\left[-t^{z-1}e^{-t}\right]_{t=0}^{\infty}+(z-1)\int_{t=0}^{\infty}{t^{z-2}e^{-t}}\,dt=(z-1)\Gamma(z-1)}
Γ
(
1
)
=
∫
t
=
0
∞
e
−
t
d
t
=
[
−
e
−
t
]
t
=
0
∞
=
1
{\displaystyle \Gamma(1)=\int_{t=0}^{\infty}{e^{-t}}\,dt=\left[-e^{-t}\right]_{t=0}^{\infty}=1}
であり、自然数
n
{\displaystyle n}
について
Γ
(
n
+
1
)
=
n
!
{\displaystyle \Gamma(n+1)=n!}
が成り立つ。従って、ガンマ関数は階乗 の定義域を複素平面 に拡張したものといえる。そのような関数は無数に存在するが、正の実軸上で対数凸 である解析関数という条件を付ければ、それは一意に定まりガンマ関数に他ならない(→ボーア・モレルップの定理 )。右半平面においてオイラー積分で定義されたガンマ関数は全平面に有理型 に解析接続 する。ガンマ関数は零点 を持たず、原点と負の整数に一位の極を持つ。その留数 は、
R
e
s
(
Γ
,
−
n
)
=
(
−
1
)
n
n
!
{\displaystyle {\rm Res}(\Gamma , -n) = \frac{(-1)^n}{n!} }
である。また、非整数でのガンマ関数の値のうちでおそらく最も有名なのは、ガウス積分 になる以下の場合であろう。
Γ
(
1
2
)
=
π
{\displaystyle \Gamma\left(\frac{1}{2}\right)=\sqrt{\pi}}
これより、自然数
n
{\displaystyle n}
について
Γ
(
1
2
+
n
)
=
(
2
n
−
1
)
!
!
2
n
π
{\displaystyle
\Gamma\left(\frac{1}{2}+n\right) = \frac{(2n-1)!!}{2^{n}}\sqrt{\pi}
}
が成立することがわかる。ここで !! は二重階乗 を表す。この性質を利用して高次元 球 の体積と表面積を求めることができる。また、
Γ
(
1
2
−
n
)
=
(
−
2
)
n
(
2
n
−
1
)
!
!
π
{\displaystyle
\Gamma\left(\frac{1}{2}-n\right) = \frac{ (-2)^n}{(2n-1)!!} \sqrt{\pi}
}
定義の整合性 [ ]
定義の積分表示と乗積表示が一致することを示す。
G
n
(
z
)
=
∫
0
n
t
z
−
1
(
1
−
t
n
)
n
d
t
{\displaystyle G_n(z)=\int_{0}^{n}{t^{z-1}\left(1-\frac{t}{n}\right)^{n}}dt}
とすれば
lim
n
→
∞
(
1
−
t
/
n
)
n
=
e
−
t
{\displaystyle \lim_{n\to\infty}{(1-t/n)^n}=e^{-t}}
であるから
lim
n
→
∞
G
n
(
z
)
=
∫
0
∞
t
z
−
1
e
−
t
d
t
{\displaystyle \lim_{n\to\infty}{G_n(z)}=\int^{\infin}_{0}t^{z-1}e^{-t}dt}
である。
t
=
n
u
{\displaystyle t=nu}
の置換により
G
n
(
z
)
=
n
z
∫
0
1
u
z
−
1
(
1
−
u
)
n
d
u
{\displaystyle G_n(z)=n^{z}\int_{0}^{1}{u^{z-1}(1-u)^{n}}du}
n
z
{\displaystyle n^{z}}
を除く部分を
g
n
(
z
)
{\displaystyle g_n(z)}
として
g
0
(
z
)
=
∫
0
1
u
z
−
1
d
u
=
[
u
z
z
]
u
=
0
1
=
1
z
{\displaystyle g_0(z)=\int_{0}^{1}{u^{z-1}}du=\left[\frac{u^z}{z}\right]_{u=0}^{1}=\frac{1}{z}}
g
n
(
z
)
=
∫
0
1
(
u
z
z
)
′
(
1
−
u
)
n
d
u
=
n
z
∫
u
=
0
1
u
z
(
1
−
u
)
n
−
1
d
u
=
n
z
g
n
−
1
(
z
+
1
)
{\displaystyle g_n(z)=\int_{0}^{1}{\left(\frac{u^{z}}{z}\right)'(1-u)^{n}}du=\frac{n}{z}\int_{u=0}^{1}{u^{z}(1-u)^{n-1}}du=\frac{n}{z}g_{n-1}(z+1)}
これにより
G
n
(
z
)
=
n
z
n
!
∏
k
=
0
n
(
z
+
k
)
{\displaystyle G_n(z)=\frac{n^zn!}{\prod_{k=0}^{n}{(z+k)}}}
を得る。故に
∫
0
∞
t
z
−
1
e
−
t
d
t
=
lim
n
→
∞
G
n
(
z
)
=
lim
n
→
∞
n
z
n
!
∏
k
=
0
n
(
z
+
k
)
{\displaystyle \int^{\infin}_{0}t^{z-1}e^{-t}dt=\lim_{n\to\infty}G_n(z)=\lim_{n\to\infty}\frac{n^zn!}{\prod_{k=0}^{n}{(z+k)}}}
である。
ワイエルシュトラスの乗積表示 [ ]
オイラーの乗積表示からオイラーの定数
γ
=
lim
n
→
∞
(
∑
k
=
1
n
1
k
−
log
n
)
{\displaystyle \gamma=\lim_{n\to\infty}\left(\sum_{k=1}^{n}\frac{1}{k}-\log{n}\right)}
を括り出すと、
ワイエルシュトラスの乗積表示 が得られる。ワイエルシュトラス はガンマ関数が負の整数に極 を持つことを嫌って逆数を用いた。ガンマ関数の逆数は複素平面全体で正則 である。
1
Γ
(
z
)
=
lim
n
→
∞
∏
k
=
0
n
(
z
+
k
)
n
z
n
!
=
lim
n
→
∞
z
n
−
z
(
∏
k
=
1
n
e
z
/
k
)
(
∏
m
=
1
n
z
+
m
m
e
−
z
/
m
)
=
z
e
γ
z
∏
m
=
1
∞
(
1
+
z
m
)
e
−
z
/
m
{\displaystyle \begin{align}
\frac{1}{\Gamma(z)}&=\lim_{n\to\infty}\frac{\prod_{k=0}^{n}{(z+k)}}{n^zn!}=\lim_{n\to\infty}zn^{-z}\left(\prod_{k=1}^{n}{e^{z/k}}\right)\left(\prod_{m=1}^{n}{\frac{z+m}{m}}e^{-z/m}\right) \\
&=ze^{{\gamma}z}\prod_{m=1}^{\infty}\left(1+\frac{z}{m}\right)e^{-z/m} \\
\end{align}}
ハンケルの積分表示 [ ]
ガンマ関数は次の周回積分で表される[8] 。積分経路は正の無限大から実軸の上側に沿って原点に至り、原点を正の向きに回り、実軸の下側に沿って無限大に戻るものとする。但し、その偏角は
−
π
≤
arg
(
−
t
)
≤
π
,
0
≤
arg
(
s
)
≤
2
π
{\displaystyle -\pi\le\arg(-t)\le\pi,0\le\arg(s)\le2\pi}
とする。
Γ
(
z
)
=
i
2
sin
π
z
∫
C
(
−
t
)
z
−
1
e
−
t
d
t
(
z
∈
C
∖
Z
)
Γ
(
z
)
=
1
e
2
π
i
z
−
1
∫
C
s
z
−
1
e
−
s
d
s
(
z
∈
C
∖
Z
)
1
Γ
(
z
)
=
i
2
π
∫
C
(
−
t
)
−
z
e
−
t
d
t
(
z
∈
C
)
{\displaystyle \begin{align}
&\Gamma(z)=\frac{i}{2\sin{\pi}z}\int_C(-t)^{z-1}e^{-t}dt\qquad(z\in\mathbb{C}\setminus\mathbb{Z})\\
&\Gamma(z)=\frac{1}{e^{2{\pi}iz}-1}\int_Cs^{z-1}e^{-s}ds\qquad(z\in\mathbb{C}\setminus\mathbb{Z})\\
&\frac{1}{\Gamma(z)}=\frac{i}{2\pi}\int_C(-t)^{-z}e^{-t}dt\qquad(z\in\mathbb{C})\\
\end{align}}
これをハンケルの積分表示 と呼ぶ。このハンケルの積分表示は、積分経路を適当に変形し、数値積分でガンマ関数の値を求めるために使われることがある[9] 。
ハンケルの積分表示の導出 [ ]
極座標表示
(
−
t
)
=
r
e
i
θ
{\displaystyle (-t)=re^{i\theta}}
を用いると、実軸の上側に沿う部分は
θ
=
−
π
{\displaystyle \theta=-\pi}
で
r
=
∞
{\displaystyle r=\infty}
から
r
=
δ
{\displaystyle r=\delta}
まで、原点を回る部分は
r
=
δ
{\displaystyle r=\delta}
で
θ
=
−
π
{\displaystyle \theta=-\pi}
から
θ
=
π
{\displaystyle \theta=\pi}
まで、実軸の下側に沿う部分は
θ
=
π
{\displaystyle \theta=\pi}
で
r
=
δ
{\displaystyle r=\delta}
から
r
=
∞
{\displaystyle r=\infty}
までとなる。
∫
C
(
−
t
)
z
−
1
e
−
t
d
t
=
∫
∞
δ
(
r
e
−
π
i
)
z
−
1
e
−
r
d
r
+
∫
−
π
π
(
δ
e
i
θ
)
z
−
1
e
δ
e
i
θ
(
−
i
δ
e
i
θ
)
d
θ
+
∫
δ
∞
(
r
e
π
i
)
z
−
1
e
−
r
d
r
=
∫
∞
δ
r
z
−
1
e
−
π
i
(
z
−
1
)
e
−
r
d
r
−
∫
−
π
π
i
δ
z
e
i
θ
z
e
δ
e
i
θ
d
θ
+
∫
δ
∞
r
z
−
1
e
π
i
(
z
−
1
)
e
−
r
d
r
=
(
−
e
−
π
i
(
z
−
1
)
+
e
π
i
(
z
−
1
)
)
∫
δ
∞
r
z
−
1
e
−
r
d
r
−
∫
−
π
π
i
δ
z
e
i
θ
z
e
δ
e
i
θ
d
θ
=
−
2
i
sin
π
z
∫
δ
∞
r
z
−
1
e
−
r
d
r
−
∫
−
π
π
i
δ
z
e
i
θ
z
e
δ
e
i
θ
d
θ
{\displaystyle \begin{align}
&\int_C(-t)^{z-1}e^{-t}dt\\
&=\int_{\infty}^{\delta}(re^{-{\pi}i})^{z-1}e^{-r}dr+\int_{-\pi}^{\pi}({\delta}e^{i\theta})^{z-1}e^{{\delta}e^{i\theta}}(-i{\delta}e^{i\theta})d\theta+\int_{\delta}^{\infty}(re^{{\pi}i})^{z-1}e^{-r}dr\\
&=\int_{\infty}^{\delta}r^{z-1}e^{-{\pi}i(z-1)}e^{-r}dr-\int_{-\pi}^{\pi}i\delta^ze^{i{\theta}z}e^{{\delta}e^{i\theta}}d\theta+\int_{\delta}^{\infty}r^{z-1}e^{{\pi}i(z-1)}e^{-r}dr\\
&=\left(-e^{-{\pi}i(z-1)}+e^{{\pi}i(z-1)}\right)\int_{\delta}^{\infty}r^{z-1}e^{-r}dr-\int_{-\pi}^{\pi}i\delta^ze^{i{\theta}z}e^{{\delta}e^{i\theta}}d\theta\\
&=-2i\sin{\pi}z\int_{\delta}^{\infty}r^{z-1}e^{-r}dr-\int_{-\pi}^{\pi}i\delta^ze^{i{\theta}z}e^{{\delta}e^{i\theta}}d\theta\\
\end{align}}
ℜ
z
>
0
{\displaystyle \real{z}>0}
とすると
δ
→
0
{\displaystyle \delta \to 0}
で
δ
z
→
0
{\displaystyle \delta^z\to0}
であるから
∫
C
(
−
t
)
z
−
1
e
−
t
d
t
=
−
2
i
sin
π
z
∫
0
∞
r
z
−
1
e
−
r
d
r
=
−
2
i
sin
π
z
Γ
(
z
)
(
ℜ
z
>
0
)
{\displaystyle \begin{align}\int_C(-t)^{z-1}e^{-t}dt
&=-2i\sin{\pi}z\int_{0}^{\infty}r^{z-1}e^{-r}dr\\
&=-2i\sin{\pi}z\Gamma(z)\qquad(\real{z}>0)\\
\end{align}}
である。しかし、左辺の被積分関数は
z
{\displaystyle z}
が有界であるかぎり正則であるから、左辺は複素平面全体に解析接続する。従って、
Γ
(
z
)
=
i
2
sin
π
z
∫
C
(
−
t
)
z
−
1
e
−
t
d
t
(
z
∈
C
∖
Z
)
{\displaystyle \Gamma(z)=\frac{i}{2\sin{\pi}z}\int_C(-t)^{z-1}e^{-t}dt\qquad(z\in\mathbb{C}\setminus\mathbb{Z})}
である。
s
=
r
e
i
θ
{\displaystyle s=re^{i\theta}}
とすれば、同様にして
Γ
(
z
)
=
1
e
2
π
i
z
−
1
∫
C
s
z
−
1
e
−
t
d
s
(
z
∈
C
∖
Z
)
{\displaystyle \Gamma(z)=\frac{1}{e^{2{\pi}iz}-1}\int_Cs^{z-1}e^{-t}ds\qquad(z\in\mathbb{C}\setminus\mathbb{Z})}
を得る。また、相半公式 により、
1
Γ
(
z
)
=
sin
π
z
π
Γ
(
1
−
z
)
=
i
2
π
∫
C
(
−
t
)
−
z
e
−
t
d
t
(
z
∈
C
)
{\displaystyle \frac{1}{\Gamma(z)}=\frac{\sin{\pi}z}{\pi}\Gamma(1-z)=\frac{i}{2\pi}\int_C(-t)^{-z}e^{-t}dt\qquad(z\in\mathbb{C})}
を得る。
相半公式 [ ]
次の恒等式をオイラーの相半公式 (reflection formula)という。
Γ
(
z
)
Γ
(
1
−
z
)
=
−
z
Γ
(
z
)
Γ
(
−
z
)
=
π
sin
π
z
{\displaystyle \Gamma(z)\Gamma(1-z)=-z\Gamma(z)\Gamma(-z)=\frac{\pi}{\sin{{\pi}z}}}
この恒等式はオイラーの乗積表示から得られる。
−
z
Γ
(
z
)
Γ
(
−
z
)
=
−
z
(
lim
n
→
∞
n
z
n
!
∏
k
=
0
n
(
z
+
k
)
)
(
lim
n
→
∞
n
−
z
n
!
∏
k
=
0
n
(
−
z
+
k
)
)
=
1
z
∏
k
=
1
∞
k
2
k
2
−
z
2
=
π
π
z
∏
k
=
1
∞
k
2
−
z
2
k
2
{\displaystyle \begin{align}
-z\Gamma(z)\Gamma(-z)
&=-z\left(\lim_{n\to\infty}\frac{n^zn!}{\prod_{k=0}^{n}{(z+k)}}\right)\left(\lim_{n\to\infty}\frac{n^{-z}n!}{\prod_{k=0}^{n}{(-z+k)}}\right)\\
&=\frac{1}{z}\prod_{k=1}^{\infty}\frac{k^2}{k^2-z^2}\\
&=\frac{\pi}{{\pi}z\displaystyle\prod_{k=1}^{\infty}\displaystyle\frac{k^2-z^2}{k^2}}\\
\end{align}}
この分母は正弦関数の無限乗積展開 であるから、
Γ
(
z
)
Γ
(
1
−
z
)
=
−
z
Γ
(
z
)
Γ
(
−
z
)
=
π
sin
π
z
{\displaystyle \Gamma(z)\Gamma(1-z)=-z\Gamma(z)\Gamma(-z)=\frac{\pi}{\sin{{\pi}z}}}
である。反射公式に
z
=
1
2
{\displaystyle z=\frac{1}{2}}
を代入すれば
Γ
(
1
2
)
Γ
(
1
−
1
2
)
=
π
sin
π
2
=
π
{\displaystyle \Gamma\left(\frac{1}{2}\right)\Gamma\left(1-\frac{1}{2}\right)=\frac{\pi}{\sin{\frac{\pi}{2}}}=\pi}
となり
Γ
(
1
2
)
=
π
{\displaystyle \Gamma\left(\frac{1}{2}\right)=\sqrt{\pi}}
を得る。
乗法公式 [ ]
次の恒等式をガウスの乗法公式 (multiplication formula)という。
Γ
(
n
z
)
=
n
n
z
−
1
/
2
(
2
π
)
(
n
−
1
)
/
2
∏
k
=
0
n
−
1
Γ
(
z
+
k
n
)
{\displaystyle \Gamma(nz)=\frac{n^{nz-1/2}}{(2\pi)^{(n-1)/2}}\prod_{k=0}^{n-1}{\Gamma{\left(z+\frac{k}{n}\right)}}}
この証明を示す。両辺の比を
f
(
z
)
{\displaystyle f(z)}
とすると
f
(
z
)
=
n
n
z
−
1
/
2
∏
k
=
0
n
−
1
Γ
(
z
+
k
n
)
(
2
π
)
(
n
−
1
)
/
2
Γ
(
n
z
)
{\displaystyle \begin{align}f(z)=
&\frac{n^{nz-1/2}\prod_{k=0}^{n-1}{\Gamma{\left(z+\frac{k}{n}\right)}}}{(2\pi)^{(n-1)/2}\Gamma(nz)}\\
\end{align}}
f
(
z
+
1
)
=
n
n
z
−
1
/
2
n
n
[
∏
k
=
0
n
−
1
(
z
+
k
n
)
Γ
(
z
+
k
n
)
]
(
2
π
)
(
n
−
1
)
/
2
[
∏
k
=
0
n
−
1
(
n
z
+
k
)
]
Γ
(
n
z
)
=
n
n
z
−
1
/
2
[
∏
k
=
0
n
−
1
(
n
z
+
k
)
]
∏
k
=
0
n
−
1
Γ
(
z
+
k
n
)
(
2
π
)
(
n
−
1
)
/
2
[
∏
k
=
0
n
−
1
(
n
z
+
k
)
]
Γ
(
n
z
)
=
f
(
z
)
{\displaystyle \begin{align}f(z+1)
&=\frac{n^{nz-1/2}n^n\left[\prod_{k=0}^{n-1}\left(z+\frac{k}{n}\right)\Gamma{\left(z+\frac{k}{n}\right)}\right]}{(2\pi)^{(n-1)/2}\left[\prod_{k=0}^{n-1}(nz+k)\right]\Gamma(nz)}\\
&=\frac{n^{nz-1/2}\left[\prod_{k=0}^{n-1}\left(nz+k\right)\right]\prod_{k=0}^{n-1}\Gamma{\left(z+\frac{k}{n}\right)}}{(2\pi)^{(n-1)/2}\left[\prod_{k=0}^{n-1}(nz+k)\right]\Gamma(nz)}\\
&=f(z)\\
\end{align}}
故に、任意に大きな自然数
m
{\displaystyle m}
について
f
(
z
+
m
)
=
f
(
z
)
{\displaystyle f(z+m)=f(z)}
が成立する。スターリングの公式により
lim
ℜ
z
→
+
∞
f
(
z
)
=
lim
ℜ
z
→
+
∞
n
n
z
−
1
/
2
[
∏
k
=
0
n
−
1
2
π
z
+
k
/
n
(
z
+
k
/
n
e
)
z
+
k
/
n
]
(
2
π
)
(
n
−
1
)
/
2
2
π
n
z
(
n
z
e
)
n
z
=
lim
ℜ
z
→
+
∞
z
1
/
2
[
∏
k
=
0
n
−
1
z
k
/
n
−
1
/
2
(
1
+
k
/
n
z
)
z
+
k
/
n
−
1
/
2
e
−
k
/
n
]
=
lim
ℜ
z
→
+
∞
z
1
/
2
[
∏
k
=
0
n
−
1
z
k
/
n
−
1
/
2
e
k
/
n
e
−
k
/
n
]
=
1
{\displaystyle \begin{align}\lim_{\real{z}\to+\infty}f(z)
&=\lim_{\real{z}\to+\infty}\frac{n^{nz-1/2}\left[\prod_{k=0}^{n-1}{\sqrt{\frac{2{\pi}}{z+k/n}}\left(\frac{z+k/n}{e}\right)^{z+k/n}}\right]}{(2\pi)^{(n-1)/2}\sqrt{\frac{2{\pi}}{nz}}\left(\frac{nz}{e}\right)^{nz}}\\
&=\lim_{\real{z}\to+\infty}z^{1/2}\left[\prod_{k=0}^{n-1}z^{k/n-1/2}(1+k/nz)^{z+k/n-1/2}e^{-k/n}\right]\\
&=\lim_{\real{z}\to+\infty}z^{1/2}\left[\prod_{k=0}^{n-1}z^{k/n-1/2}e^{k/n}e^{-k/n}\right]\\
&=1
\end{align}}
途中で
lim
ℜ
z
→
+
∞
(
1
+
k
/
n
z
)
z
+
k
/
n
−
1
/
2
=
lim
ℜ
z
→
+
∞
(
1
+
k
/
n
z
)
z
=
e
n
/
k
{\displaystyle \lim_{\real{z}\to+\infty}(1+k/nz)^{z+k/n-1/2}=\lim_{\real{z}\to+\infty}(1+k/nz)^{z}=e^{n/k}}
を適用した。
f
(
z
)
=
lim
n
→
∞
f
(
z
+
n
)
=
1
{\displaystyle f(z)=\lim_{n\to\infty}f(z+n)=1}
であり、故に
Γ
(
n
z
)
=
n
n
z
−
1
/
2
(
2
π
)
(
n
−
1
)
/
2
∏
k
=
0
n
−
1
Γ
(
z
+
k
n
)
{\displaystyle \Gamma(nz)=\frac{n^{nz-1/2}}{(2\pi)^{(n-1)/2}}\prod_{k=0}^{n-1}{\Gamma{\left(z+\frac{k}{n}\right)}}}
が成立する。
微分方程式 [ ]
(
x
,
y
,
y
1
,
…
,
y
n
)
{\displaystyle (x,\ y,\ y_1,\ \ldots ,\ y_n)}
を変数とする多項式
F
(
x
,
y
,
y
1
,
…
,
y
n
)
{\displaystyle F(x,\ y,\ y_1,\ \ldots ,\ y_n)}
に対し、
F
(
x
,
y
,
y
1
,
⋯
,
y
n
)
=
0
,
y
i
=
d
i
y
d
x
i
(
i
=
1
,
⋯
,
n
)
{\displaystyle F(x, y, y_1, \cdots, y_n)=0 , \quad y_i=\frac{d^i y}{dx^i} \quad (i=1, \cdots,n)}
の形で表される微分方程式を代数的微分方程式という。ガンマ関数はいかなる代数的微分方程式も満たさないことが知られている。ヘルダーが1887年に最初に証明を与えた後
[10] 、ムーア[11] 、オストロフスキー[12] [13] 、バーンズ[14] 、ハウスドルフ[15] により、別証明や一般化がなされた。
いくつかの具体的な値 [ ]
Γ
(
−
3
2
)
=
4
π
3
≈
2.363
{\displaystyle \Gamma\left(-\frac{3}{2}\right)\,= \frac {4\sqrt{\pi}} {3} \approx 2.363\,}
Γ
(
−
1
2
)
=
−
2
π
≈
−
3.545
{\displaystyle \Gamma\left(-\frac{1}{2}\right)\,= -2\sqrt{\pi} \approx -3.545\,}
Γ
(
1
2
)
=
π
≈
1.772
{\displaystyle \Gamma\left(\frac{1}{2}\right)\,= \sqrt{\pi} \approx 1.772\,}
Γ
(
1
)
=
0
!
=
1
{\displaystyle \Gamma(1)\,=0!=1 \,}
Γ
(
3
2
)
=
π
2
≈
0.886
{\displaystyle \Gamma\left(\frac{3}{2}\right)\,= \frac {\sqrt{\pi}} {2} \approx 0.886\,}
Γ
(
2
)
=
1
!
=
1
{\displaystyle \Gamma(2)\,=1!=1 \,}
Γ
(
5
2
)
=
3
π
4
≈
1.329
{\displaystyle \Gamma\left(\frac{5}{2}\right)\,= \frac {3 \sqrt{\pi}} {4} \approx 1.329\,}
Γ
(
3
)
=
2
!
=
2
{\displaystyle \Gamma(3)\,=2!=2 \,}
Γ
(
7
2
)
=
15
π
8
≈
3.323
{\displaystyle \Gamma\left(\frac{7}{2}\right)\,= \frac {15\sqrt{\pi}} {8} \approx 3.323\,}
Γ
(
4
)
=
3
!
=
6
{\displaystyle \Gamma(4)\,=3!=6 \,}
Γ
(
n
2
)
=
π
(
n
−
2
)
!
!
2
(
n
−
1
)
/
2
,
{\displaystyle \Gamma \left (\tfrac{n}{2} \right) = \sqrt \pi \frac{(n-2)!!}{2^{(n-1)/2}}\,,}
Γ
(
1
2
+
n
)
=
(
2
n
−
1
)
!
!
2
n
π
=
(
2
n
)
!
4
n
n
!
π
Γ
(
1
2
−
n
)
=
(
−
2
)
n
(
2
n
−
1
)
!
!
π
=
(
−
4
)
n
n
!
(
2
n
)
!
π
{\displaystyle \begin{align}
\Gamma\left(\tfrac{1}{2}+n\right) &= \frac{(2n-1)!!}{2^n}\, \sqrt{\pi} = \frac{(2n)!}{4^n n!} \sqrt{\pi} \\
\Gamma\left(\tfrac{1}{2}-n\right) &= \frac{(-2)^n}{(2n-1)!!}\, \sqrt{\pi} = \frac{(-4)^n n!}{(2n)!} \sqrt{\pi}
\end{align}}
Γ
(
n
+
1
3
)
=
Γ
(
1
3
)
(
3
n
−
2
)
!
(
3
)
3
n
Γ
(
n
+
1
4
)
=
Γ
(
1
4
)
(
4
n
−
3
)
!
(
4
)
4
n
Γ
(
n
+
1
p
)
=
Γ
(
1
p
)
(
p
n
−
(
p
−
1
)
)
!
(
p
)
p
n
{\displaystyle \begin{align}
\Gamma \left(n+\tfrac13 \right) &= \Gamma \left(\tfrac13 \right) \frac{(3n-2)!^{(3)}}{3^n} \\
\Gamma \left(n+\tfrac14 \right) &= \Gamma \left(\tfrac14 \right ) \frac{(4n-3)!^{(4)}}{4^n} \\
\Gamma \left(n+\tfrac{1}{p} \right ) &= \Gamma \left(\tfrac{1}{p} \right ) \frac{(pn-(p-1))!^{(p)}}{p^n}
\end{align}}
Γ
(
1
3
)
≈
2.6789385347077476337
{\displaystyle \Gamma(\tfrac13) \approx 2.6789385347077476337}
A073005
Γ
(
1
4
)
≈
3.6256099082219083119
{\displaystyle \Gamma(\tfrac14) \approx 3.6256099082219083119}
A068466
Γ
(
1
5
)
≈
4.5908437119988030532
{\displaystyle \Gamma(\tfrac15) \approx 4.5908437119988030532}
A175380
Γ
(
1
6
)
≈
5.5663160017802352043
{\displaystyle \Gamma(\tfrac16) \approx 5.5663160017802352043}
A175379
Γ
(
1
7
)
≈
6.5480629402478244377
{\displaystyle \Gamma(\tfrac17) \approx 6.5480629402478244377}
A220086
Γ
(
1
8
)
≈
7.5339415987976119047
{\displaystyle \Gamma(\tfrac18) \approx 7.5339415987976119047}
A203142 .
Γ
(
1
4
)
=
2
π
S
,
{\displaystyle \Gamma(\tfrac14) = \sqrt{\sqrt{2 \pi} S},}
Γ
(
1
4
)
=
(
4
π
3
e
2
γ
−
δ
+
1
)
1
4
{\displaystyle \Gamma \left (\tfrac14 \right ) = \left(4 \pi^3 e^{2 \gamma -\mathrm{\delta}+1}\right)^{\frac14}}
Γ
(
1
4
)
=
(
2
π
)
3
4
∏
k
=
1
∞
tanh
(
π
k
2
)
{\displaystyle \Gamma(\tfrac14) = (2 \pi)^{\frac{3}{4}} \prod_{k=1}^\infty \tanh \left( \frac{\pi k}{2} \right)}
Γ
(
1
4
)
=
A
3
e
−
G
π
π
2
1
6
∏
k
=
1
∞
(
1
−
1
2
k
)
k
(
−
1
)
k
{\displaystyle \Gamma(\tfrac14) = A^3 e^{-\frac{G}{\pi}} \sqrt{\pi} 2^{\frac{1}{6}} \prod_{k=1}^\infty \left(1-\frac{1}{2k}\right)^{k(-1)^k}}
(
Γ
(
1
3
)
)
6
12
π
4
=
1
10
∑
k
=
0
∞
(
6
k
)
!
(
−
1
)
k
(
k
!
)
3
(
3
k
)
!
3
k
160
3
k
(
Γ
(
1
4
)
)
4
128
π
3
=
1
u
∑
k
=
0
∞
(
6
k
)
!
(
2
w
)
k
(
k
!
)
3
(
3
k
)
!
6486
3
k
{\displaystyle \begin{align}
\frac{\left(\Gamma(\tfrac13) \right )^6}{12\pi^4} &= \frac{1}{\sqrt{10}} \sum_{k = 0}^{\infty} \frac{(6k)!(-1)^k}{(k!)^{3}(3k)! 3^{k}160^{3k}} \\
\frac{\left(\Gamma(\tfrac14) \right )^4}{128\pi^3} &= \frac{1}{\sqrt{u}}\, \sum_{k = 0}^{\infty} \frac{(6k)!(2w)^k}{(k!)^{3}(3k)! 6486^{3k}}
\end{align}}
u
=
273
+
180
2
v
=
1
+
2
w
=
−
761354780
+
538359129
2
=
6486
3
2
(
u
v
2
2
)
3
{\displaystyle \begin{align}
u &= 273+180\sqrt{2}\\
v &= 1+\sqrt{2}\\
w &= -761354780+538359129\sqrt{2} = \frac{6486^3}{2\bigl(uv^2\sqrt{2}\bigr)^3}
\end{align}}
(
Γ
(
1
4
)
)
4
128
π
3
=
1
u
∑
k
=
0
∞
(
6
k
)
!
(
k
!
)
3
(
3
k
)
!
1
(
u
v
2
2
)
3
k
.
{\displaystyle \frac{\left(\Gamma(\tfrac14) \right )^4}{128\pi^3} = \frac{1}{\sqrt{u}} \sum_{k = 0}^{\infty} \frac{(6k)!}{(k!)^{3}(3k)!} \frac{1}{(uv^2\sqrt{2})^{3k}}.}
∏
r
=
1
2
Γ
(
r
3
)
=
2
π
3
≈
3.6275987284684357012
{\displaystyle \prod_{r=1}^2 \Gamma(\tfrac{r}{3}) = \frac{2\pi}{\sqrt{3}} \approx 3.6275987284684357012}
A186706
∏
r
=
1
3
Γ
(
r
4
)
=
2
π
3
≈
7.8748049728612098721
{\displaystyle \prod_{r=1}^3 \Gamma(\tfrac{r}{4}) = \sqrt{2\pi^3} \approx 7.8748049728612098721}
A220610
∏
r
=
1
4
Γ
(
r
5
)
=
4
π
2
5
≈
17.6552850814935242483
{\displaystyle \prod_{r=1}^4 \Gamma(\tfrac{r}{5}) = \frac{4\pi^2}{\sqrt{5}} \approx 17.6552850814935242483}
∏
r
=
1
5
Γ
(
r
6
)
=
4
π
5
3
≈
40.3993191220037900785
{\displaystyle \prod_{r=1}^5 \Gamma(\tfrac{r}{6}) = 4\sqrt{\frac{\pi^5}{3}} \approx 40.3993191220037900785}
∏
r
=
1
6
Γ
(
r
7
)
=
8
π
3
7
≈
93.7541682035825037970
{\displaystyle \prod_{r=1}^6 \Gamma(\tfrac{r}{7}) = \frac{8\pi^3}{\sqrt{7}} \approx 93.7541682035825037970}
∏
r
=
1
7
Γ
(
r
8
)
=
4
π
7
≈
219.8287780169572636207
{\displaystyle \prod_{r=1}^7 \Gamma(\tfrac{r}{8}) = 4\sqrt{\pi^7} \approx 219.8287780169572636207}
∏
r
=
1
n
Γ
(
r
n
+
1
)
=
(
2
π
)
n
n
+
1
{\displaystyle \prod_{r=1}^n \Gamma(\tfrac{r}{n+1}) = \sqrt{\frac{(2\pi)^n}{n+1}}}
Γ
(
1
5
)
Γ
(
4
15
)
Γ
(
1
3
)
Γ
(
2
15
)
=
2
3
20
5
6
5
−
7
5
+
6
−
6
5
4
{\displaystyle \frac{\Gamma(\tfrac{1}{5})\Gamma(\tfrac{4}{15})}{\Gamma(\tfrac{1}{3})\Gamma(\tfrac{2}{15})} = \frac{\sqrt{2}\sqrt[20]{3}}{\sqrt[6]{5} \sqrt[4]{5-\frac{7}{\sqrt{5}}+\sqrt{6-\frac{6}{\sqrt{5}}}}}}
[16]
Γ
(
1
20
)
Γ
(
9
20
)
Γ
(
3
20
)
Γ
(
7
20
)
=
5
4
(
1
+
5
)
2
{\displaystyle \frac{\Gamma(\tfrac{1}{20})\Gamma(\tfrac{9}{20})}{\Gamma(\tfrac{3}{20})\Gamma(\tfrac{7}{20})} = \frac{\sqrt[4]{5}\left(1+\sqrt{5}\right)}{2}}
[17]
Γ
(
i
)
=
G
(
1
+
i
)
G
(
i
)
=
e
−
log
G
(
i
)
+
log
G
(
1
+
i
)
.
{\displaystyle \Gamma(i) = \frac{G(1+i)}{G(i)} = e^{-\log G(i)+ \log G(1+i)}.}
Γ
(
1
+
i
)
=
i
Γ
(
i
)
≈
0.498
−
0.155
i
{\displaystyle \Gamma(1 + i) = i\Gamma(i) \approx 0.498 - 0.155i}
Γ
(
1
−
i
)
=
−
i
Γ
(
−
i
)
≈
0.498
+
0.155
i
{\displaystyle \Gamma(1 - i) = -i\Gamma(-i) \approx 0.498 + 0.155i}
Γ
(
0.5
+
0.5
i
)
≈
0.8181639995
−
0.7633138287
i
{\displaystyle \Gamma(0.5 + 0.5i) \approx 0.8181639995 - 0.7633138287 i}
Γ
(
0.5
−
0.5
i
)
≈
0.8181639995
+
0.7633138287
i
{\displaystyle \Gamma(0.5 - 0.5i) \approx 0.8181639995 + 0.7633138287 i}
Γ
(
5
+
3
i
)
≈
0.0160418827
−
9.4332932898
i
{\displaystyle \Gamma(5 + 3i) \approx 0.0160418827 - 9.4332932898 i}
Γ
(
5
−
3
i
)
≈
0.0160418827
+
9.4332932897
i
.
{\displaystyle \Gamma(5 - 3i) \approx 0.0160418827 + 9.4332932897 i.}
不完全ガンマ関数 [ ]
数学 において、不完全ガンマ関数 (ふかんぜん-かんすう incomplete gamma function)あるいは、ルジャンドルの不完全ガンマ関数 は、ガンマ関数の一般化の一つ。ガンマ関数は定積分 を用いて定義されるが、不完全ガンマ関数は不定積分 を用いて定義される。
定義 [ ]
不完全ガンマ関数には2種類あり、ガンマ関数の積分区間[0,∞]を2つに分けて以下のように定義される。
0以上の実数 x と、 実部 が正の複素数 a に対し
第1種不完全ガンマ関数
γ
(
a
,
x
)
{\displaystyle \gamma(a,x)}
γ
(
a
,
x
)
=
∫
0
x
t
a
−
1
e
−
t
d
t
{\displaystyle \gamma(a,x) = \int_0^x t^{a-1}\,e^{-t}\,dt \,\!}
第2種不完全ガンマ関数
Γ
(
a
,
x
)
{\displaystyle \Gamma(a,x)}
Γ
(
a
,
x
)
=
∫
x
∞
t
a
−
1
e
−
t
d
t
{\displaystyle \Gamma(a,x) = \int_x^{\infty} t^{a-1}\,e^{-t}\,dt \,\!}
性質 [ ]
ガンマ関数の定義は
Γ
(
a
)
=
∫
0
∞
t
a
−
1
e
−
t
d
t
{\displaystyle \Gamma(a) = \int_0^{\infty} t^{a-1}\,e^{-t}\,dt \,\!}
であるから、
γ
(
a
,
x
)
+
Γ
(
a
,
x
)
=
Γ
(
a
)
{\displaystyle \gamma(a,x) + \Gamma(a,x) = \Gamma(a) \, }
となる。
また、不完全ガンマ関数の定義式に部分積分 を用いることで
γ
(
a
+
1
,
x
)
=
a
γ
(
a
,
x
)
−
x
a
e
−
x
{\displaystyle \gamma(a+1,x) = a\gamma(a,x) - x^a e^{-x}\,}
Γ
(
a
+
1
,
x
)
=
a
Γ
(
a
,
x
)
+
x
a
e
−
x
{\displaystyle \Gamma(a+1,x) = a\Gamma(a,x) + x^a e^{-x}\,}
という関係が成り立つことも分かる。
さらに、以下のような式が成り立つ。
Γ
(
a
,
0
)
=
Γ
(
a
)
{\displaystyle \Gamma(a,0) = \Gamma(a)\, }
γ
(
a
,
x
)
→
Γ
(
a
)
a
s
x
→
∞
{\displaystyle \gamma(a,x) \rightarrow \Gamma(a)
\quad \mathrm{as\ } x \rightarrow \infty \,}
Γ
(
0
,
x
)
=
−
Ei
(
−
x
)
for
x
>
0
{\displaystyle \Gamma(0,x) = -\mbox{Ei}(-x)\mbox{ for }x>0 \,}
Γ
(
1
2
,
x
)
=
π
erfc
(
x
)
{\displaystyle \Gamma\left({1 \over 2}, x\right) = \sqrt\pi\,\mbox{erfc}\left(\sqrt x\right) \,}
γ
(
1
2
,
x
)
=
π
erf
(
x
)
{\displaystyle \gamma\left({1 \over 2}, x\right) = \sqrt\pi\,\mbox{erf}\left(\sqrt x\right) \,}
Γ
(
1
,
x
)
=
e
−
x
{\displaystyle \Gamma(1,x) = e^{-x} \,}
γ
(
1
,
x
)
=
1
−
e
−
x
{\displaystyle \gamma(1,x) = 1 - e^{-x} \,}
ここで、
Ei : 指数積分
erf : 誤差関数
erfc :相補誤差関数 で、erfc(x ) = 1 − erf(x )
であるとする。
微分法 [ ]
∂
Γ
(
a
,
x
)
∂
x
=
−
x
a
−
1
e
x
{\displaystyle \frac{\partial \Gamma (a,x) }{\partial x} = - \frac{x^{a-1}}{e^x}}
MeijerのG関数から[18] :
T
(
m
,
a
,
x
)
=
G
m
−
1
,
m
m
,
0
(
0
,
0
,
…
,
0
a
−
1
,
−
1
,
…
,
−
1
|
x
)
{\displaystyle T(m,a,x) = G_{m-1,\,m}^{\,m,\,0} \!\left( \left. \begin{matrix} 0, 0, \dots, 0 \\ a-1, -1, \dots, -1 \end{matrix} \; \right| \, x \right)}
T
(
m
,
a
,
z
)
=
−
(
−
1
)
m
−
1
(
m
−
2
)
!
d
m
−
2
d
t
m
−
2
[
Γ
(
a
−
t
)
z
t
−
1
]
|
t
=
0
+
∑
n
=
0
∞
(
−
1
)
n
z
a
−
1
+
n
n
!
(
−
a
−
n
)
m
−
1
{\displaystyle T(m,a,z) = - \frac{(-1)^{m-1} }{(m-2)! } \frac{{\rm d}^{m-2} }{{\rm d}t^{m-2} } \left[\Gamma (a-t) z^{t-1}\right]\Big|_{t=0} + \sum_{n=0}^{\infty} \frac{(-1)^n z^{a-1+n}}{n! (-a-n)^{m-1} }}
時
|
z
|
<
1
{\displaystyle |z| < 1}
∂
Γ
(
a
,
x
)
∂
a
=
ln
x
Γ
(
a
,
x
)
+
x
T
(
3
,
a
,
x
)
{\displaystyle \frac{\partial \Gamma (a,x) }{\partial a} = \ln x \Gamma (a,x) + x\,T(3,a,x) }
∂
2
Γ
(
a
,
x
)
∂
a
2
=
ln
2
x
Γ
(
a
,
x
)
+
2
x
[
ln
x
T
(
3
,
a
,
x
)
+
T
(
4
,
a
,
x
)
]
{\displaystyle \frac{\partial^2 \Gamma (a,x) }{\partial a^2} = \ln^2 x \Gamma (a,x) + 2 x[\ln x\,T(3,a,x) + T(4,a,x) ]}
∂
m
Γ
(
a
,
x
)
∂
a
m
=
ln
m
x
Γ
(
a
,
x
)
+
m
x
∑
n
=
0
m
−
1
P
n
m
−
1
ln
m
−
n
−
1
x
T
(
3
+
n
,
a
,
x
)
{\displaystyle \frac{\partial^m \Gamma (a,x) }{\partial a^m} = \ln^m x \Gamma (a,x) + m x\,\sum_{n=0}^{m-1} P_n^{m-1} \ln^{m-n-1} x\,T(3+n,a,x)}
と
P
j
n
=
(
n
j
)
j
!
=
n
!
(
n
−
j
)
!
.
{\displaystyle P_j^n = \left( \begin{array}{l} n \\ j \end{array} \right) j! = \frac{n!}{(n-j)!}.}
ポリガンマ関数 [ ]
曖昧さ回避
この項目では、ガンマ関数の対数微分で定義されるポリガンマ関数(polygamma function)について説明しています。E. Barnesによって導入された多重ガンマ関数(multiple gamma function)については「多重ガンマ関数 」をご覧ください。
実数x に対するψ(n) (x )の挙動。 オレンジがディガンマ関数、黄色がトリガンマ関数、緑がテトラガンマ関数、赤がペンタガンマ関数、青がヘキサガンマ関数に対応する。
複素平面上でのディガンマ関数ψ(z)
複素平面上でのトリガンマ関数ψ(1) (z)
複素平面上でのテトラガンマ関数ψ(2) (z)
複素平面上でのペンタガンマ関数ψ(3) (z)
数学において、ポリガンマ関数 (ぽりがんまかんすう、英: polygamma function )とはガンマ関数 の対数微分 による導関数 として定義される特殊関数 。ディガンマ関数 やトリガンマ関数 はポリガンマ関数の一種である。
定義 [ ]
ガンマ関数Γ (z )に対し、その対数微分
ψ
(
n
)
(
x
)
=
d
n
+
1
d
z
n
+
1
ln
Γ
(
z
)
=
d
n
d
z
n
ψ
(
z
)
{\displaystyle
\psi^{(n)}(x)
= \frac{d^{n+1}}{dz^{n+1}} \ln{\Gamma(z)} = \frac{d^n}{dz^n} \psi(z)
}
で、定義される関数をポリガンマ関数 と呼ぶ。
ψ(z ), ψ(1) (z ), ψ(2) (z ), ψ(3) (z )、 ψ(4) (z )は、それぞれディ -、トリ -、テトラ -、ペンタ -、ヘキサ -ガンマ関数 と呼ばれる。
ポリガンマ関数ψ(n) (z )はz =0,-1,-2,……でn+1位の極 をもち,それらの点を除く全複素平面 では解析的 になる。
漸化式 [ ]
ポリガンマ関数は次の漸化式 を満たす。
ψ
(
n
)
(
z
+
1
)
=
ψ
(
n
)
(
z
)
+
(
−
1
)
n
n
!
z
n
+
1
{\displaystyle
\psi^{(n)}(z+1) = \psi^{(n)}(z) + \frac{(-1)^n n!}{z^{n+1}}
}
級数表示 [ ]
ポリガンマ関数はz ≠0, -1, -2, -3...で次の級数 表示を持つ。
ψ
(
z
)
=
−
γ
+
∑
n
=
0
∞
(
1
z
+
n
−
1
n
+
1
)
{\displaystyle
\psi(z) = -\gamma +\sum_{n=0}^{\infty} \biggl (
\frac{1}{z+n} - \frac{1}{n+1} \biggr )
}
ψ
(
n
)
(
z
)
=
(
−
1
)
n
+
1
n
!
∑
k
=
0
∞
1
(
z
+
k
)
n
+
1
(
n
=
1
,
2
,
3
,
⋯
)
{\displaystyle
\psi^{(n)}(z) = (-1)^{n+1} n! \sum_{k=0}^{\infty} \frac{1}{(z+k)^{n+1}}
\qquad (n = 1,2,3, \cdots)
}
また、z =0でのテイラー展開 により、|z |<1の領域で次のように表される。
ψ
(
z
+
1
)
=
−
γ
+
∑
k
=
2
∞
(
−
1
)
k
ζ
(
k
)
z
k
−
1
{\displaystyle
\psi(z+1) = -\gamma +\sum_{k=2}^{\infty} (-1)^k
\zeta(k) z^{k-1}
}
ψ
(
n
)
(
z
+
1
)
=
(
−
1
)
n
+
1
∑
k
=
1
∞
(
−
1
)
k
−
1
(
n
+
k
−
1
)
!
ζ
(
n
+
k
)
z
k
−
1
(
k
−
1
)
!
(
n
=
1
,
2
,
3
,
⋯
)
{\displaystyle
\psi^{(n)}(z+1) = (-1)^{n+1} \sum_{k=1}^{\infty} \frac{(-1)^{k-1}
(n+k-1)! \zeta(n+k) z^{k-1} }{(k-1)!}
\qquad (n = 1,2,3, \cdots)
}
但し、γ =0.5772...はオイラーの定数 、ζ (n )はリーマンゼータ関数 を表す。
積分表示 [ ]
Rez >0のとき、ポリガンマ関数は次の積分表示を持つ。
ψ
(
z
)
=
−
γ
+
∫
0
∞
e
−
t
−
e
−
z
t
1
−
e
−
t
d
t
{\displaystyle
\psi(z) = -\gamma+ \int_{0}^{\infty} \frac{e^{-t}-e^{-zt}}{1-e^{-t}} dt
}
ψ
(
n
)
(
z
)
=
(
−
1
)
n
+
1
∫
0
∞
t
n
e
−
z
t
1
−
e
−
t
d
t
(
n
=
1
,
2
,
⋯
)
{\displaystyle
\psi^{(n)}(z) = (-1)^{n+1} \int_{0}^{\infty} \frac{t^ne^{-zt}}{1-e^{-t}} dt
\quad (n=1,2,\cdots)
}
相反公式 [ ]
ガンマ関数の相反公式に対し、対数微分をとることで次の関係式が導かれる。
ψ
(
n
)
(
1
−
z
)
+
(
−
1
)
n
+
1
ψ
(
n
)
(
z
)
=
π
d
n
d
z
n
cot
π
z
{\displaystyle
\psi^{(n)}(1-z) + (-1)^{n+1} \psi^{(n)}(z) = \pi \frac{d^n}{dz^n} \operatorname{cot} \pi z
}
但し、cot πz は余接関数 を表す。
漸近展開 [ ]
z →∞ (|argz | < π)のとき、ポリガンマ関数は次の漸近展開 をもつ。
ψ
(
z
)
∼
ln
z
−
1
2
z
−
∑
n
=
1
∞
B
2
n
2
n
z
2
n
{\displaystyle
\psi(z) \sim \ln{z}- \frac{1}{2z}- \sum_{n=1}^{\infty} \frac{B_{2n}}{2nz^{2n}}
}
ψ
(
n
)
(
z
)
∼
(
−
1
)
(
n
−
1
)
(
(
n
−
1
)
!
z
n
+
n
!
2
z
n
+
1
+
∑
k
=
1
∞
B
2
k
(
2
k
+
n
−
1
)
!
(
2
k
)
!
z
2
k
+
n
)
(
n
=
1
,
2
,
⋯
)
{\displaystyle
\psi^{(n)}(z) \sim (-1)^{(n-1)} \left( \frac{(n-1)!}{z^n}
+ \frac{n!}{2z^{n+1}}+ \sum_{k=1}^{\infty} \frac{B_{2k}(2k+n-1)!}{(2k)!z^{2k+n}} \right )
\quad (n=1,2,\cdots)
}
但し、B2k はベルヌーイ数 である。
特殊値 [ ]
ポリガンマ関数は、m=1において、次の値をとる。
ψ
(
1
)
=
−
γ
{\displaystyle
\psi(1) = -\gamma
}
ψ
(
n
)
(
1
)
=
(
−
1
)
n
+
1
n
!
ζ
(
n
+
1
)
(
n
=
1
,
2
,
⋯
)
{\displaystyle
\psi^{(n)}(1) = (-1)^{n+1} n! \zeta(n+1)
\quad (n=1,2,\cdots)
}
ポリガンマ関数は、m≧2の正の整数において、次の値をとる。
ψ
(
m
)
=
−
γ
+
∑
k
=
1
m
−
1
1
k
=
−
γ
+
H
m
−
1
(
m
=
2
,
3
,
4
,
⋯
)
{\displaystyle
\psi(m) = -\gamma + \sum_{k=1}^{m-1}\frac{1}{k}
= -\gamma + H_{m-1} \qquad (m = 2,3,4, \cdots)
}
ψ
(
n
)
(
m
)
=
(
−
1
)
n
n
!
{
−
ζ
(
n
+
1
)
+
∑
k
=
1
m
−
1
1
k
n
+
1
}
(
n
=
1
,
2
,
3
,
⋯
,
m
=
2
,
3
,
4
,
⋯
)
{\displaystyle \begin{align}
&\psi^{(n)}(m)
= (-1)^n n! \left \{ -\zeta(n+1)
+ \sum_{k=1}^{m-1}\frac{1}{k^{n+1}} \right \} \\
&(n = 1,2,3, \cdots,m = 2,3,4, \cdots) \\
\end{align}}
但し、γ はオイラーの定数、Hm-1 は調和数 を表す。
また、以下の等式が成り立つ。
8
G
+
π
2
+
ψ
1
(
1
4
)
=
0
{\displaystyle 8G+\pi^2+\psi_1\left(\frac{1}{4}\right)=0}
8
G
+
ψ
1
(
3
4
)
−
6
L
2
(
1
6
)
−
6
L
2
(
5
6
)
=
0
{\displaystyle 8G+\psi_1\left(\frac{3}{4}\right)-6L_2\left(\frac{1}{6}\right)-6L_2\left(\frac{5}{6}\right)=0}
ψ
1
(
5
6
)
−
5
ψ
1
(
2
3
)
+
8
L
2
(
1
6
)
−
+
8
L
2
(
5
6
)
=
0
{\displaystyle \psi_1\left(\frac{5}{6}\right)-5\psi_1\left(\frac{2}{3}\right)+8L_2\left(\frac{1}{6}\right)-+8L_2\left(\frac{5}{6}\right)=0}
ψ
1
(
1
6
)
+
ψ
1
(
5
6
)
−
4
π
2
=
0
{\displaystyle \psi_1\left(\frac{1}{6}\right)+\psi_1\left(\frac{5}{6}\right)-4\pi^2=0}
ψ
1
(
1
4
)
+
ψ
1
(
3
4
)
−
2
π
2
=
0
{\displaystyle \psi_1\left(\frac{1}{4}\right)+\psi_1\left(\frac{3}{4}\right)-2\pi^2=0}
3
ψ
1
(
1
6
)
−
15
ψ
1
(
1
3
)
+
4
π
2
=
0
{\displaystyle 3\psi_1\left(\frac{1}{6}\right)-15\psi_1\left(\frac{1}{3}\right)+4\pi^2=0}
−
ψ
1
(
5
6
)
−
ψ
1
(
1
3
)
+
4
ψ
1
(
2
3
)
=
0
{\displaystyle -\psi_1\left(\frac{5}{6}\right)-\psi_1\left(\frac{1}{3}\right)+4\psi_1\left(\frac{2}{3}\right)=0}
−
32
G
+
4
ψ
1
(
1
4
)
−
3
ψ
1
(
1
3
)
−
3
ψ
1
(
2
3
)
=
0
{\displaystyle -32G+4\psi_1\left(\frac{1}{4}\right)-3\psi_1\left(\frac{1}{3}\right)-3\psi_1\left(\frac{2}{3}\right)=0}
3
π
2
−
8
G
+
ψ
1
(
1
4
)
−
3
ψ
1
(
1
3
)
−
3
ψ
1
(
2
3
)
=
0
{\displaystyle 3\pi^2-8G+\psi_1\left(\frac{1}{4}\right)-3\psi_1\left(\frac{1}{3}\right)-3\psi_1\left(\frac{2}{3}\right)=0}
−
2
ψ
1
(
1
4
)
−
2
ψ
1
(
3
4
)
+
ψ
1
(
5
6
)
+
ψ
1
(
1
6
)
=
0
{\displaystyle -2\psi_1\left(\frac{1}{4}\right)-2\psi_1\left(\frac{3}{4}\right)+\psi_1\left(\frac{5}{6}\right)+\psi_1\left(\frac{1}{6}\right)=0}
−
4
π
3
−
ψ
2
(
1
4
)
+
ψ
2
(
3
4
)
=
0
{\displaystyle -4\pi^3-\psi_2\left(\frac{1}{4}\right)+\psi_2\left(\frac{3}{4}\right)=0}
−
64
π
3
3
−
9
ψ
2
(
1
6
)
+
63
ψ
2
(
2
3
)
=
0
{\displaystyle -64\pi^3\sqrt{3}-9\psi_2\left(\frac{1}{6}\right)+63\psi_2\left(\frac{2}{3}\right)=0}
−
52
ζ
(
3
)
−
9
ψ
2
(
2
3
)
+
ψ
2
(
5
6
)
=
0
{\displaystyle -52\zeta(3)-9\psi_2\left(\frac{2}{3}\right)+\psi_2\left(\frac{5}{6}\right)=0}
8
π
3
3
+
52
ζ
(
3
)
−
54
ψ
2
(
1
3
)
+
ψ
2
(
1
12
)
+
ψ
2
(
7
12
)
=
0
{\displaystyle 8\pi^3\sqrt{3}+52\zeta(3)-54\psi_2\left(\frac{1}{3}\right)+\psi_2\left(\frac{1}{12}\right)+\psi_2\left(\frac{7}{12}\right)=0}
−
8
ψ
2
(
1
4
)
+
ψ
2
(
1
8
)
+
ψ
2
(
5
8
)
=
0
{\displaystyle -8\psi_2\left(\frac{1}{4}\right)+\psi_2\left(\frac{1}{8}\right)+\psi_2\left(\frac{5}{8}\right)=0}
階乗の逆数和 [ ]
階乗の逆数 の総和は収斂級数
∑
n
=
0
∞
1
n
!
=
1
1
+
1
1
+
1
2
+
1
6
+
1
24
+
1
120
+
⋯
=
e
{\displaystyle \sum_{n=0}^{\infty} \frac{1}{n!} = \frac{1}{1} + \frac{1}{1} + \frac{1}{2} + \frac{1}{6} + \frac{1}{24} + \frac{1}{120} + \dotsb = e}
を与える(ネイピア数 を参照)。この和は無理数 となるけれども、階乗に適当な正整数を掛けて和が有理数となるようにすることができる:
∑
n
=
0
∞
1
(
n
+
2
)
n
!
=
1
2
+
1
3
+
1
8
+
1
30
+
1
144
+
⋯
=
1.
{\displaystyle \sum_{n=0}^{\infty} \frac{1}{(n+2)n!} = \frac{1}{2}+\frac{1}{3}+\frac{1}{8}+\frac{1}{30}+\frac{1}{144}+\dotsb=1.}
この級数の値が 1 となることを見るには、その部分和 が 1 − 1/(n +2)! であることを確認すればよい。したがって、階乗数の全体は無理列を成さない[19] 。
階乗の増大度 [ ]
n が増えるにつれて、階乗 n ! は n を変数とする任意の多項式函数 あるいは指数函数 よりも早く増加する(ただし、二重指数関数 よりは遅い)。
n ! の近似式の多くは自然対数
log
n
!
=
∑
x
=
1
n
log
x
{\displaystyle \log n! = \sum_{x=1}^n \log x}
に基づく。もっとも単純に得られる log(n !) の近似値を評価する式は、上記の式と以下の積分:
∫
1
n
log
x
d
x
≤
∑
x
=
1
n
log
x
≤
∫
0
n
log
(
x
+
1
)
d
x
{\displaystyle \int_1^n \log x \, dx \leq \sum_{x=1}^n \log x \leq \int_0^n \log (x+1) \, dx}
によって与えられる。積分を評価すれば
n
log
(
n
e
)
+
1
≤
log
n
!
≤
(
n
+
1
)
log
(
n
+
1
e
)
+
1
{\displaystyle n\log\left(\frac{n}{e}\right)+1 \leq \log n! \leq (n+1)\log\left( \frac{n+1}{e} \right) + 1}
を得る。これは、ランダウの記号 を用いれば log(n !) のオーダーは Θ(n log n ) であることを言っているのであり、この結果はソートアルゴリズム の計算量 を測るのに重要な役割を果たす。さて上記の log(n !) の評価から
e
(
n
e
)
n
≤
n
!
≤
e
(
n
+
1
e
)
n
+
1
{\displaystyle e\left(\frac ne\right)^n \leq n! \leq e\left(\frac{n+1}e\right)^{n+1}}
がわかる。実用上はより弱い結果だがより評価のしやすいものを用いることもある。上記の式から簡単な評価をしてみると、任意の n に対して (n /3)n < n ! であり、また n ≥ 6 のとき n ! < (n /2)n であることなどが分かる。
大きな n に対して n ! をもっとましに評価するにはスターリングの公式
n
!
≈
2
π
n
(
n
e
)
n
{\displaystyle n!\approx\sqrt{2\pi n}\left(\frac{n}{e}\right)^n}
を利用する。実は任意の n に対し
n
!
>
2
π
n
(
n
e
)
n
{\displaystyle n! > \sqrt{2\pi n}\left(\frac{n}{e}\right)^n}
であることが証明できる[要出典 ] 。
log(n !) の別な近似はシュリニヴァーサ・ラマヌジャンにより
log
n
!
≈
n
log
n
−
n
+
log
(
n
(
1
+
4
n
(
1
+
2
n
)
)
)
6
+
log
(
π
)
2
{\displaystyle \log n! \approx n\log n - n + \frac {\log(n(1+4n(1+2n)))} {6} + \frac {\log(\pi)} {2}}
=
n
log
n
−
n
+
log
(
1
+
1
/
(
2
n
)
+
1
/
(
8
n
2
)
)
6
+
log
(
2
n
)
2
+
log
(
π
)
2
.
{\displaystyle = n\log n - n + \frac {\log(1 +1/(2n) +1/(8n^2))} {6} + \frac {\log (2n)} {2} + \frac {\log(\pi)} {2}.}
と与えられている[20] 。この近似の誤差は、スターリングの公式の次の誤差項 1 ⁄12n よりも小さい。
スターリングの近似 [ ]
(ln n !) と (n ln n − n ) は n が増大するにつれて近づいていく
スターリングの近似 (英: Stirling's approximation )またはスターリングの公式 (英: Stirling's formula )は、階乗 、或いはガンマ関数の漸近近似 である。数学者ジェイムズ・スターリングに因む。
応用上よく使われる形の公式は、ランダウの記号 を用いて、
log
(
n
!
)
=
n
log
(
n
)
−
n
+
O
(
log
(
n
)
)
{\displaystyle \log(n!) = n\log(n) - n +O(\log(n))}
である。O (log(n )) における次の項は (1/2)log(2πn ) である。よって公式のより正確な形は
n
!
∼
2
π
n
(
n
e
)
n
{\displaystyle n! \sim \sqrt{2\pi n} \left(\frac{n}{e}\right)^n}
である。これは漸近公式 であり、スターリングの近似は
lim
n
→
∞
n
!
2
π
n
(
n
/
e
)
n
=
1
{\displaystyle \lim_{n\to\infty} \frac{n!}{\sqrt{2\pi n} (n/e)^n} = 1}
という性質を持つ。
n ! の漸近近似よりもむしろ上下からの評価が必要なことがある。任意の正の整数 n に対して、
2
π
n
n
+
1
/
2
e
−
n
≤
n
!
≤
e
n
n
+
1
/
2
e
−
n
{\displaystyle \sqrt{2\pi}\ n^{n+1/2}e^{-n} \le n! \le e\ n^{n+1/2}e^{-n} }
が成り立ち、従って任意の n ≥ 1 に対して比
n
!
n
n
+
1
/
2
e
−
n
{\displaystyle \frac{n!}{n^{n+1/2}e^{-n}}}
は常に
2
π
=
2.5066...
{\displaystyle \sqrt{2\pi} = 2.5066...}
A019727 と
e
=
2.71828...
{\displaystyle e = 2.71828...}
A001113 の間にある。
ガンマ関数
Γ
(
z
)
{\displaystyle \Gamma(z)}
は正の整数 n に対して
Γ
(
n
)
=
(
n
−
1
)
!
{\displaystyle \Gamma(n)=(n-1)!}
を満たす。スターリングの近似はガンマ関数に対しても考えることができ
Γ
(
z
+
1
)
≈
2
π
z
(
z
e
)
z
(
|
arg
z
|
<
π
,
|
z
|
≫
0
)
{\displaystyle \Gamma(z+1)\approx\sqrt{2\pi z}\left(\frac{z}{e}\right)^z\qquad(|\arg z|<\pi,|z|\gg0)}
が成り立つ。これは
lim
z
→
∞
Γ
(
z
+
1
)
2
π
z
(
z
/
e
)
z
=
1
(
|
arg
z
|
<
π
)
{\displaystyle \lim_{z\to\infty}\frac{\Gamma(z+1)}{\sqrt{2\pi z}(z/e)^z}=1\qquad(|\arg z|<\pi)}
を意味する。但し、
|
arg
z
|
=
π
{\displaystyle |\arg z|=\pi}
に近づくにつれ近似の誤差が大きくなる(極限の収束が遅くなる)ため、応用上は反射公式などを用いて
|
arg
z
|
≤
π
/
2
{\displaystyle |\arg z|\le\pi/2}
程度に制限することが多い。
導出 [ ]
初等的な方法による導出 [ ]
スターリングの公式の導出を厳密にするためにはオイラーの和公式 、或いは鞍点法 など複素解析の技法を用いる必要があるが、多少の無理を許せば初等的に導くことも可能である。先ず、階乗の対数 を積分で近似する。
log
(
n
!
)
=
∑
k
=
1
n
log
k
≈
∑
k
=
1
n
∫
k
−
1
/
2
k
+
1
/
2
log
x
d
x
=
∫
1
/
2
n
+
1
/
2
log
x
d
x
{\displaystyle \log(n!)=\sum_{k=1}^{n}{\log k}
\approx\sum_{k=1}^{n}\int_{k-1/2}^{k+1/2}\log x\, dx
=\int_{1/2}^{n+1/2}\log x\, dx}
∫
1
/
2
n
+
1
/
2
log
x
d
x
=
[
x
log
x
−
x
]
1
/
2
n
+
1
/
2
=
(
n
+
1
2
)
log
(
n
+
1
2
)
−
(
n
+
1
2
)
−
(
1
2
)
log
(
1
2
)
+
1
2
{\displaystyle \begin{align}
&\int_{1/2}^{n+1/2}\log x\, dx=\left[x\log{x}-x\right]_{1/2}^{n+1/2} \\
=&\left(n+\frac{1}{2}\right) \log\left(n+\frac{1}{2}\right)-\left(n+\frac{1}{2}\right)
-\left(\frac{1}{2}\right) \log\left(\frac{1}{2}\right)+\frac{1}{2} \\
\end{align}}
然し、
log
1
−
∫
1
/
2
3
/
2
log
x
d
x
≃
0.045
log
2
−
∫
3
/
2
5
/
2
log
x
d
x
≃
0.011
log
3
−
∫
5
/
2
7
/
2
log
x
d
x
≃
0.005
{\displaystyle \begin{align}
\log1-\int_{1/2}^{3/2}\log x\, dx &\simeq0.045 \\
\log2-\int_{3/2}^{5/2}\log x\, dx &\simeq0.011 \\
\log3-\int_{5/2}^{7/2}\log x\, dx &\simeq0.005 \\
\end{align} }
などという誤差があるから、定数項と誤差をまとめて log c とし、c の値は後に求める。
log
(
n
!
)
≈
log
c
+
(
n
+
1
2
)
log
(
n
+
1
2
)
−
(
n
+
1
2
)
{\displaystyle \log(n!)\approx\log{c}+\left(n+\frac{1}{2}\right)\log\left(n+\frac{1}{2}\right)-\left(n+\frac{1}{2}\right)}
両辺の真数を取ると
n
!
≈
c
(
n
+
1
2
)
n
+
1
/
2
(
1
e
)
n
+
1
/
2
=
c
e
−
1
/
2
(
n
+
1
2
)
1
/
2
(
1
+
1
2
n
)
n
(
n
e
)
n
{\displaystyle n!\approx c\left(n+\frac{1}{2}\right)^{n+1/2}\left(\frac{1}{e}\right)^{n+1/2}
=ce^{-1/2}\left(n+\frac{1}{2}\right)^{1/2}\left(1+\frac{1}{2n}\right)^n \left(\frac{n}{e}\right)^n}
となり
lim
n
→
∞
(
n
+
1
2
)
1
/
2
=
n
1
/
2
=
n
{\displaystyle \lim_{n\to\infty}\left(n+\frac{1}{2}\right)^{1/2}=n^{1/2}=\sqrt{n}}
lim
n
→
∞
(
1
+
1
2
n
)
n
=
e
1
/
2
{\displaystyle \lim_{n\to\infty}\left(1+\frac{1}{2n}\right)^n=e^{1/2}}
であるから
n
!
≈
c
n
(
n
e
)
n
{\displaystyle n!\approx c\sqrt{n}\left(\frac{n}{e}\right)^n}
を得る。ウォリス積 を変形して
π
2
=
∏
k
=
1
∞
(
2
k
)
2
(
2
k
+
1
)
(
2
k
−
1
)
≈
1
2
n
+
1
(
(
2
n
)
!
!
(
2
n
−
1
)
!
!
)
2
=
1
2
n
+
1
(
(
2
n
)
!
!
(
2
n
)
!
!
(
2
n
)
!
)
2
=
2
4
n
(
n
!
)
4
(
2
n
+
1
)
(
(
2
n
)
!
)
2
{\displaystyle \begin{align}
\frac{\pi}{2}=&\prod_{k=1}^{\infty}\frac{(2k)^2}{(2k+1)(2k-1)}\approx\frac{1}{2n+1}\left(\frac{(2n)!!}{(2n-1)!!}\right)^2=\frac{1}{2n+1}\left(\frac{(2n)!!(2n)!!}{(2n)!}\right)^2 \\
=&\frac{2^{4n}(n!)^4}{(2n+1)((2n)!)^2} \\
\end{align}}
これを代入すると
π
2
=
2
4
n
c
4
n
2
(
n
/
e
)
4
n
(
2
n
+
1
)
c
2
(
2
n
)
(
2
n
/
e
)
4
n
=
c
2
n
2
(
2
n
+
1
)
(
2
n
)
≈
c
2
4
{\displaystyle \frac{\pi}{2}=\frac{2^{4n}c^4n^2(n/e)^{4n}}{(2n+1)c^2(2n)(2n/e)^{4n}}=\frac{c^2n^2}{(2n+1)(2n)}\approx\frac{c^2}{4}}
n
!
≈
2
π
n
(
n
e
)
n
{\displaystyle n!\approx\sqrt{2\pi n}\left(\frac{n}{e}\right)^n}
を得る。
精度を改善するために
n
!
≈
2
π
n
(
n
e
)
n
(
1
+
α
n
+
β
n
2
+
γ
n
3
)
{\displaystyle n!\approx\sqrt{2\pi n}\left(\frac{n}{e}\right)^n\left(1+\frac{\alpha}{n}+\frac{\beta}{n^2}+\frac{\gamma}{n^3}\right)}
とすると
(
n
+
1
)
!
≈
2
π
(
n
+
1
)
(
n
+
1
e
)
n
+
1
(
1
+
α
n
+
1
+
β
(
n
+
1
)
2
+
γ
(
n
+
1
)
3
)
≈
2
π
(
n
+
1
)
(
n
+
1
e
)
n
+
1
(
1
+
α
n
−
α
n
2
+
α
n
3
−
α
n
4
+
β
n
2
−
2
β
n
3
+
3
β
n
4
+
γ
n
3
−
3
γ
n
3
+
O
(
1
n
5
)
)
{\displaystyle \begin{align}(n+1)!
&\approx\sqrt{2\pi{(n+1)}}\left(\frac{n+1}{e}\right)^{n+1}\left(
1+\frac{\alpha}{n+1}+\frac{\beta}{(n+1)^2}+\frac{\gamma}{(n+1)^3}\right)\\
&\approx\sqrt{2\pi{(n+1)}}\left(\frac{n+1}{e}\right)^{n+1}\left(
1+\frac{\alpha}{n}-\frac{\alpha}{n^2}+\frac{\alpha}{n^3}-\frac{\alpha}{n^4}
+\frac{\beta}{n^2}-\frac{2\beta}{n^3}+\frac{3\beta}{n^4}
+\frac{\gamma}{n^3}-\frac{3\gamma}{n^3}+O\left(\frac{1}{n^5}\right)\right)
\end{align}}
となる。然し
(
n
+
1
)
!
=
(
n
+
1
)
n
!
≈
(
n
+
1
)
2
π
n
(
n
e
)
n
(
1
+
α
n
+
β
n
2
+
γ
n
3
)
{\displaystyle (n+1)!=(n+1)n!\approx(n+1)\sqrt{2\pi n}\left(\frac{n}{e}\right)^n
\left(1+\frac{\alpha}{n}+\frac{\beta}{n^2}+\frac{\gamma}{n^3}\right)}
であるから
(
n
+
1
n
)
n
+
1
/
2
1
e
(
1
+
α
n
−
α
n
2
+
α
n
3
−
α
n
4
+
β
n
2
−
2
β
n
3
+
3
β
n
4
+
γ
n
3
−
3
γ
n
3
+
O
(
1
n
5
)
)
=
1
+
α
n
+
β
n
2
+
γ
n
3
{\displaystyle \left(\frac{n+1}{n}\right)^{n+1/2}\frac{1}{e}
\left(1+\frac{\alpha}{n}-\frac{\alpha}{n^2}+\frac{\alpha}{n^3}-\frac{\alpha}{n^4}
+\frac{\beta}{n^2}-\frac{2\beta}{n^3}+\frac{3\beta}{n^4}
+\frac{\gamma}{n^3}-\frac{3\gamma}{n^3}+O\left(\frac{1}{n^5}\right)\right)
=1+\frac{\alpha}{n}+\frac{\beta}{n^2}+\frac{\gamma}{n^3}}
となるはずである。再び対数を用いる。
(
n
+
1
2
)
log
(
n
+
1
n
)
−
1
≈
(
n
+
1
2
)
(
1
n
−
1
2
n
2
+
1
3
n
3
−
1
4
n
4
+
1
5
n
5
−
O
(
1
n
6
)
)
−
1
≈
1
12
n
2
−
1
12
n
3
+
3
40
n
4
+
O
(
1
n
5
)
{\displaystyle \begin{align}\left(n+\frac{1}{2}\right)\log\left(\frac{n+1}{n}\right)-1
&\approx\left(n+\frac{1}{2}\right)
\left(\frac{1}{n}-\frac{1}{2n^2}+\frac{1}{3n^3}-\frac{1}{4n^4}+\frac{1}{5n^5}-O\left(\frac{1}{n^6}\right)\right)-1\\
&\approx\frac{1}{12n^2}-\frac{1}{12n^3}+\frac{3}{40n^4}+O\left(\frac{1}{n^5}\right)
\end{align}}
(
n
+
1
n
)
n
+
1
/
2
1
e
≈
1
+
1
12
n
2
−
1
12
n
3
+
3
40
n
4
+
1
2
(
12
n
2
)
2
+
O
(
1
n
5
)
{\displaystyle \left(\frac{n+1}{n}\right)^{n+1/2}\frac{1}{e}
\approx1+\frac{1}{12n^2}-\frac{1}{12n^3}+\frac{3}{40n^4}+\frac{1}{2(12n^2)^2}+O\left(\frac{1}{n^5}\right)}
これを左辺に代入すると
1
+
α
n
+
1
n
2
(
−
α
+
β
+
1
12
)
+
1
n
3
(
13
12
α
−
2
β
+
γ
−
1
12
)
+
1
n
4
(
−
7
6
α
+
37
12
β
−
3
γ
+
113
1440
)
+
O
(
1
n
5
)
=
1
+
α
n
+
β
n
2
+
γ
n
3
{\displaystyle 1+\frac{\alpha}{n}+\frac{1}{n^2}\left( -\alpha+\beta+\frac{1}{12} \right)
+\frac{1}{n^3}\left( \frac{13}{12}\alpha-2\beta+\gamma-\frac{1}{12} \right)
+\frac{1}{n^4} \left( -\frac{7}{6}\alpha+\frac{37}{12}\beta-3\gamma+\frac{113}{1440} \right)
+O\left(\frac{1}{n^5}\right)=1+\frac{\alpha}{n}+\frac{\beta}{n^2}+\frac{\gamma}{n^3}}
となり、両辺の係数を比較して α=1/12,β=1/288,γ=-130/51840 を得る。結局、
n
!
≈
2
π
n
(
n
e
)
n
(
1
+
1
12
n
+
1
288
n
2
−
139
51840
n
3
)
{\displaystyle n!\approx\sqrt{2\pi n}\left(\frac{n}{e}\right)^n\left(1+\frac{1}{12n}+\frac{1}{288n^2}-\frac{139}{51840n^3}\right)}
を得る。
オイラーの和公式による導出 [ ]
オイラーの乗積表示によるガンマ関数の定義の対数をとり
log
Γ
(
z
)
=
log
(
(
z
−
1
)
Γ
(
z
−
1
)
)
=
lim
N
→
∞
(
z
−
1
)
log
N
+
∑
n
=
1
N
(
log
n
−
log
(
n
+
z
−
1
)
)
{\displaystyle \begin{align}\log\Gamma(z)&=\log\big((z-1)\Gamma(z-1)\big)\\
&=\lim_{N\to\infty}(z-1)\log{N}+\sum_{n=1}^{N}\big(\log{n}-\log{(n+z-1)}\big)\\
\end{align}}
f
(
n
)
=
log
n
−
log
(
n
+
z
−
1
)
{\displaystyle f(n)=\log{n}-\log{(n+z-1)}}
にオイラーの和公式 を適用すれば
log
Γ
(
z
)
=
lim
N
→
∞
(
z
−
1
)
log
N
+
∫
n
=
1
N
f
(
n
)
d
n
+
1
2
(
f
(
N
)
+
f
(
1
)
)
+
∑
k
=
1
m
B
2
k
(
2
k
)
!
(
f
(
2
k
−
1
)
(
N
)
−
f
(
2
k
−
1
)
(
1
)
)
+
∫
n
=
1
N
B
2
m
+
1
(
n
−
⌊
n
⌋
)
(
2
m
+
1
)
!
f
(
2
m
+
1
)
(
n
)
d
n
=
lim
N
→
∞
(
z
−
1
)
log
N
+
[
n
log
n
−
(
n
+
z
−
1
)
log
(
n
+
z
−
1
)
]
n
=
1
N
+
1
2
(
log
N
−
log
(
N
+
z
−
1
)
−
log
z
)
+
∑
k
=
1
m
B
2
k
(
2
k
)
(
2
k
−
1
)
(
1
N
2
k
−
1
−
1
(
N
+
z
−
1
)
2
k
−
1
−
1
+
1
z
2
k
−
1
)
+
∫
n
=
1
N
B
2
m
+
1
(
n
−
⌊
n
⌋
)
2
m
+
1
(
1
n
2
m
+
1
−
1
(
n
+
z
−
1
)
2
m
+
1
)
d
n
=
lim
N
→
∞
(
N
+
z
−
1
2
)
(
log
N
−
log
(
N
+
z
−
1
)
)
+
(
z
+
1
2
)
log
z
+
∑
k
=
1
m
B
2
k
(
2
k
)
(
2
k
−
1
)
(
1
N
2
k
−
1
−
1
(
N
+
z
−
1
)
2
k
−
1
−
1
+
1
z
2
k
−
1
)
+
∫
n
=
1
N
B
2
m
+
1
(
n
−
⌊
n
⌋
)
2
m
+
1
(
1
n
2
m
+
1
−
1
(
n
+
z
−
1
)
2
m
+
1
)
d
n
=
−
z
+
1
+
(
z
−
1
2
)
log
z
−
∑
k
=
1
m
B
2
k
(
2
k
)
(
2
k
−
1
)
(
1
−
1
z
2
k
−
1
)
+
∫
n
=
1
N
B
2
m
+
1
(
n
−
⌊
n
⌋
)
2
m
+
1
(
1
n
2
m
+
1
−
1
(
n
+
z
−
1
)
2
m
+
1
)
d
n
{\displaystyle \begin{align}\log\Gamma(z)
&=\lim_{N\to\infty}(z-1)\log{N}+\int_{n=1}^{N}f(n)dn+\frac{1}{2}\big(f(N)+f(1)\big)\\
&\qquad+\sum_{k=1}^{m}\frac{B_{2k}}{(2k)!}\left(f^{(2k-1)}(N)-f^{(2k-1)}(1)\right)+\int_{n=1}^{N}\frac{B_{2m+1}(n-\lfloor{n}\rfloor)}{(2m+1)!}f^{(2m+1)}(n)dn\\
&=\lim_{N\to\infty}(z-1)\log{N}+\bigg[{n}\log{n}-{(n+z-1)}\log(n+z-1)\bigg]_{n=1}^{N}+\frac{1}{2}\big(\log{N}-\log{(N+z-1)}-\log{z}\big)\\
&\qquad+\sum_{k=1}^{m}\frac{B_{2k}}{(2k)(2k-1)}\left(\frac{1}{N^{2k-1}}-\frac{1}{(N+z-1)^{2k-1}}-1+\frac{1}{z^{2k-1}}\right)+\int_{n=1}^{N}\frac{B_{2m+1}(n-\lfloor{n}\rfloor)}{2m+1}\left(\frac{1}{n^{2m+1}}-\frac{1}{(n+z-1)^{2m+1}}\right)dn\\
&=\lim_{N\to\infty}\left(N+z-\frac{1}{2}\right)\big(\log{N}-\log{(N+z-1)}\big)+\left(z+\frac{1}{2}\right)\log{z}\\
&\qquad+\sum_{k=1}^{m}\frac{B_{2k}}{(2k)(2k-1)}\left(\frac{1}{N^{2k-1}}-\frac{1}{(N+z-1)^{2k-1}}-1+\frac{1}{z^{2k-1}}\right)+\int_{n=1}^{N}\frac{B_{2m+1}(n-\lfloor{n}\rfloor)}{2m+1}\left(\frac{1}{n^{2m+1}}-\frac{1}{(n+z-1)^{2m+1}}\right)dn\\
&=-z+1+\left(z-\frac{1}{2}\right)\log{z}-\sum_{k=1}^{m}\frac{B_{2k}}{(2k)(2k-1)}\left(1-\frac{1}{z^{2k-1}}\right)+\int_{n=1}^{N}\frac{B_{2m+1}(n-\lfloor{n}\rfloor)}{2m+1}\left(\frac{1}{n^{2m+1}}-\frac{1}{(n+z-1)^{2m+1}}\right)dn\\
\end{align}}
となる。右辺の定数を集めて
C
=
1
−
∑
k
=
1
m
B
2
k
(
2
k
)
(
2
k
−
1
)
+
∫
n
=
1
N
B
2
m
+
1
(
n
−
⌊
n
⌋
)
d
n
(
2
m
+
1
)
n
2
m
+
1
{\displaystyle C=1-\sum_{k=1}^{m}\frac{B_{2k}}{(2k)(2k-1)}+\int_{n=1}^{N}\frac{B_{2m+1}(n-\lfloor{n}\rfloor)dn}{(2m+1)n^{2m+1}}}
とすれば
log
Γ
(
z
)
=
C
−
z
+
(
z
−
1
2
)
log
z
+
∑
k
=
1
m
B
2
k
(
2
k
)
(
2
k
−
1
)
z
2
k
−
1
−
∫
n
=
1
N
B
2
m
+
1
(
n
−
⌊
n
⌋
)
d
n
(
2
m
+
1
)
(
n
+
z
−
1
)
2
m
+
1
{\displaystyle \log\Gamma(z)=C-z+\left(z-\frac{1}{2}\right)\log{z}+\sum_{k=1}^{m}\frac{B_{2k}}{(2k)(2k-1)z^{2k-1}}-\int_{n=1}^{N}\frac{B_{2m+1}(n-\lfloor{n}\rfloor)dn}{(2m+1)(n+z-1)^{2m+1}}}
となり、この主要部をガンマ関数の反射公式に代入して
z
→
i
∞
{\displaystyle z\to{i\infty}}
とすれば
となるが
log
sin
π
z
=
log
(
e
π
i
z
−
e
−
π
i
z
)
−
π
i
2
−
log
2
≈
π
i
z
−
π
i
2
−
log
2
{\displaystyle \log\sin{{\pi}z}=\log\left(e^{{\pi}iz}-e^{-{\pi}iz}\right)-\frac{{\pi}i}{2}-\log{2}\approx{\pi}iz-\frac{{\pi}i}{2}-\log{2}}
であるから
C
=
log
2
π
2
{\displaystyle C=\frac{\log{2\pi}}{2}}
を得る。剰余項については
α
=
2
1
+
cos
(
arg
z
)
{\displaystyle \alpha=\frac{2}{1+\cos(\arg{z})}}
として
|
∫
n
=
1
N
B
2
m
+
1
(
n
−
⌊
n
⌋
)
d
n
(
2
m
+
1
)
(
n
+
z
−
1
)
2
m
+
1
|
≤
|
B
2
m
|
2
m
∫
n
=
1
N
d
n
|
n
+
z
−
1
|
2
m
+
1
≤
|
B
2
m
|
2
m
α
2
m
+
1
∫
n
=
1
N
d
n
(
n
+
|
z
|
−
1
)
2
m
+
1
=
|
B
2
m
|
α
2
m
+
1
|
z
|
2
m
=
O
(
z
−
2
m
)
{\displaystyle \begin{align}
\left|\int_{n=1}^{N}\frac{B_{2m+1}(n-\lfloor{n}\rfloor)dn}{(2m+1)(n+z-1)^{2m+1}}\right|
&\le\frac{\left|B_{2m}\right|}{2m}\int_{n=1}^{N}\frac{dn}{\left|n+z-1\right|^{2m+1}}\\
&\le\frac{\left|B_{2m}\right|}{2m\alpha^{2m+1}}\int_{n=1}^{N}\frac{dn}{(n+\left|z\right|-1)^{2m+1}}
=\frac{\left|B_{2m}\right|}{\alpha^{2m+1}\left|z\right|^{2m}}=O\left(z^{-2m}\right)\\
\end{align}}
である。故に
log
Γ
(
z
)
=
log
2
π
2
−
z
+
(
z
−
1
2
)
log
z
+
∑
k
=
1
m
B
2
k
(
2
k
)
(
2
k
−
1
)
z
2
k
−
1
+
O
(
z
−
2
m
)
{\displaystyle \log\Gamma(z)=\frac{\log{2\pi}}{2}-z+\left(z-\frac{1}{2}\right)\log{z}+\sum_{k=1}^{m}\frac{B_{2k}}{(2k)(2k-1)z^{2k-1}}+O\left(z^{-2m}\right)}
を得る。最初の数項を書き下せば
log
Γ
(
z
)
≈
log
2
π
−
z
+
(
z
−
1
2
)
log
z
+
1
12
z
−
1
360
z
3
+
1
1260
z
5
−
1
1680
z
7
+
1
1188
z
9
{\displaystyle \log\Gamma(z)\approx\log\sqrt{2\pi}-z+\left(z-\frac{1}{2}\right)\log{z}
+\frac{1}{12z}-\frac{1}{360z^3}+\frac{1}{1260z^5}-\frac{1}{1680z^7}+\frac{1}{1188z^9}}
Γ
(
z
)
≈
2
π
z
(
z
e
)
z
exp
(
1
12
z
−
1
360
z
3
+
1
1260
z
5
−
1
1680
z
7
+
1
1188
z
9
)
{\displaystyle \Gamma(z)\approx\sqrt{\frac{2\pi}{z}}\left(\frac{z}{e}\right)^z
\exp\left(\frac{1}{12z}-\frac{1}{360z^3}+\frac{1}{1260z^5}-\frac{1}{1680z^7}+\frac{1}{1188z^9}\right)}
とやり、指数関数のテイラー展開 により
Γ
(
z
)
≈
2
π
z
(
z
e
)
z
(
1
+
1
12
z
+
1
288
z
2
−
139
51840
z
3
)
{\displaystyle \Gamma(z)\approx\sqrt{\frac{2\pi}{z}}\left(\frac{z}{e}\right)^z
\left(1+\frac{1}{12z}+\frac{1}{288z^2}-\frac{139}{51840z^3}\right)}
となる。
鞍点法による導出 [ ]
スターリングの公式は鞍点法 の好適例とされることが多いが、実際に複素平面全体(負の実数を除く)で漸近近似が成立することを鞍点法によって示すのは困難であるから、ここでは
z
{\displaystyle z}
を正の実数 に限定する。ガンマ関数 は
t
=
z
(
1
+
u
)
{\displaystyle t=z(1+u)}
の置換により
Γ
(
z
+
1
)
=
∫
0
∞
t
z
e
−
t
d
t
=
∫
−
1
∞
z
z
(
1
+
u
)
z
e
−
z
−
z
u
z
d
u
=
z
z
+
1
e
−
z
∫
−
1
∞
e
−
z
(
u
−
log
(
1
+
u
)
)
d
u
=
z
z
+
1
e
−
z
(
∫
−
1
−
ϵ
e
−
z
(
u
−
log
(
1
+
u
)
)
d
u
+
∫
−
ϵ
ϵ
e
−
z
(
u
−
log
(
1
+
u
)
)
d
u
+
∫
ϵ
∞
e
−
z
(
u
−
log
(
1
+
u
)
)
d
u
)
(
ϵ
≪
1
)
{\displaystyle \begin{align}\Gamma(z+1)
&=\int_{0}^{\infty}{t^ze^{-t}dt}\\
&=\int_{-1}^{\infty}{z^z(1+u)^ze^{-z-zu}zdu}\\
&=z^{z+1}e^{-z}\int_{-1}^{\infty}{e^{-z\left(u-\log(1+u)\right)}du}\\
&=z^{z+1}e^{-z}\left(\int_{-1}^{-\epsilon}{e^{-z\left(u-\log(1+u)\right)}du}+\int_{-\epsilon}^{\epsilon}{e^{-z\left(u-\log(1+u)\right)}du}+\int_{\epsilon}^{\infty}{e^{-z\left(u-\log(1+u)\right)}du}\right)\qquad(\epsilon\ll1)\\
\end{align}}
となるが、
z
{\displaystyle z}
が十分に大きければ
u
=
0
{\displaystyle u= 0}
の附近が支配的であるから
Γ
(
z
+
1
)
≈
z
z
+
1
e
−
z
∫
−
ϵ
ϵ
e
−
z
(
u
−
log
(
1
+
u
)
)
d
u
≈
z
z
+
1
e
−
z
∫
−
ϵ
ϵ
e
−
z
u
2
/
2
d
u
≈
z
z
+
1
e
−
z
∫
−
∞
∞
e
−
z
u
2
/
2
d
u
{\displaystyle \Gamma(z+1)\approx{z^{z+1}e^{-z}\int_{-\epsilon}^{\epsilon}{e^{-z\left(u-\log(1+u)\right)}du}}\approx{z^{z+1}e^{-z}\int_{-\epsilon}^{\epsilon}{e^{-zu^2/2}du}}\approx{z^{z+1}e^{-z}\int_{-\infty}^{\infty}{e^{-zu^2/2}du}}}
という近似が許され、ガウス積分 により
Γ
(
z
+
1
)
≈
z
z
+
1
e
−
z
2
π
z
=
2
π
z
(
z
e
)
z
{\displaystyle \Gamma(z+1)\approx{z^{z+1}e^{-z}\sqrt\frac{2\pi}{z}=\sqrt{2{\pi}z}\left(\frac{z}{e}\right)^z}}
を得る。
ϵ
=
z
−
1
/
3
{\displaystyle \epsilon=z^{-1/3}}
として、近似の誤差は
であり
z
2
/
3
e
−
z
1
/
3
/
2
≪
z
−
1
/
2
(
z
→
∞
)
{\displaystyle {z^{2/3}}e^{-z^{1/3}/2}\ll{z^{-1/2}}\qquad(z\to\infty)}
であるから
Γ
(
z
+
1
)
=
2
π
z
(
z
e
)
z
(
1
+
O
(
z
−
1
/
3
)
)
{\displaystyle \Gamma(z+1)=\sqrt{2{\pi}z}\left(\frac{z}{e}\right)^z\left(1+O(z^{-1/3})\right)}
を得る。これは
lim
z
→
∞
Γ
(
z
+
1
)
2
π
z
(
z
e
)
z
=
1
(
|
arg
z
|
<
π
)
{\displaystyle \lim_{z\to\infty}\frac{\Gamma(z+1)}{\sqrt{2{\pi}z}\left(\frac{z}{e}\right)^z}=1\qquad(|\arg{z}|<{\pi})}
を示すに十分である。但し、実際の誤差は
O
(
z
−
1
)
{\displaystyle O(z^{-1})}
であるが、それを鞍点法で示すのは困難である。
収束の速度と誤差見積もり [ ]
より正確に記すと、次のようになる。
n
!
=
2
π
n
(
n
e
)
n
e
λ
n
{\displaystyle n! = \sqrt{2 \pi n} \; \left(\frac{n}{e}\right)^{n}e^{\lambda_n}}
ここで
1
12
n
+
1
<
λ
n
<
1
12
n
.
{\displaystyle \frac{1}{12n+1} < \lambda_n < \frac{1}{12n}.}
スターリングの公式は以下の級数(スターリング級数 )の近似(初項で打ち切ったもの)である。
n
!
=
2
π
n
(
n
e
)
n
(
1
+
1
12
n
+
1
288
n
2
−
139
51840
n
3
−
571
2488320
n
4
+
⋯
)
.
{\displaystyle
n!=\sqrt{2\pi n}\left({n\over e}\right)^n
\left(
1
+{1\over12n}
+{1\over288n^2}
-{139\over51840n^3}
-{571\over2488320n^4}
+ \cdots
\right).
}
n
→
∞
{\displaystyle n \to \infty}
としたとき、省かれた級数はその最初の項とそれ以降が相殺するように漸近していく。これは漸近展開 の一例である。
以下のような階乗の対数の漸近展開も「スターリング級数」と呼ぶ。
ln
n
!
=
n
ln
n
−
n
+
1
2
ln
(
2
π
n
)
+
1
12
n
−
1
360
n
3
+
1
1260
n
5
−
1
1680
n
7
+
⋯
.
{\displaystyle
\ln n!=n\ln n - n + {1\over 2}\ln(2\pi n)
+{1\over12n}
-{1\over360n^3}
+{1\over1260n^5}
-{1\over 1680n^7}
+\cdots.
}
この場合、誤差は打ち切った級数の初項と同じ符号で同程度の大きさであることが知られている。
ガンマ関数に対するスターリングの公式 [ ]
すべての正の整数に対して、
n
!
=
Π
(
n
)
=
Γ
(
n
+
1
)
{\displaystyle n! = \Pi(n) = \Gamma(n+1)}
が成り立つ。ここで Γ はガンマ関数を表す。
しかしながら、π関数
Π
(
z
)
=
Γ
(
z
+
1
)
{\displaystyle \Pi(z) = \Gamma(z+1)}
は、階乗とは異なり、より広く、正でない整数を除いてすべての複素数に対して定義される。それにもかかわらず、スターリングの公式をなお適用することができる。Re(z ) > 0 であれば
log
(
Γ
(
z
)
)
=
(
z
−
1
2
)
log
(
z
)
−
z
+
1
2
log
(
2
π
)
+
2
∫
0
∞
arctan
(
t
/
z
)
exp
(
2
π
t
)
−
1
d
t
{\displaystyle \log(\Gamma (z)) =\left(z-\tfrac{1}{2}\right)\log(z) -z + \tfrac{1}{2}\log(2 \pi)
+2\int_0^\infty \frac{\arctan(t/z)}{\exp(2\pi t)-1}\,{\rm d}t}
部分積分を繰り返すことで次が得られる
log
(
Γ
(
z
)
)
∼
(
z
−
1
2
)
log
(
z
)
−
z
+
1
2
log
(
2
π
)
+
∑
n
=
1
∞
B
2
n
2
n
(
2
n
−
1
)
z
2
n
−
1
{\displaystyle \log(\Gamma (z)) \sim \left(z-\tfrac{1}{2}\right)\log(z) -z + \tfrac{1}{2}\log(2 \pi)
+\sum_{n=1}^\infty \frac{B_{2n}}{2n(2n-1)z^{2n-1}}}
ここで Bn は n 番目のベルヌーイ数 である。(無限和は収束しないので、この公式は漸近展開 にすぎないことに注意する。)公式はεを正数として |arg(z )| < π−ε であるときに絶対値の十分大きい z に対して成り立つ。最初の m 項が使われるとき誤差項は
O
(
z
−
2
m
−
1
)
{\displaystyle O(z^{-2m-1})}
である。対応する近似は
Γ
(
z
)
=
2
π
z
(
z
e
)
z
(
1
+
O
(
1
z
)
)
{\displaystyle \Gamma(z) = \sqrt{\frac{2 \pi}{z}}~{\left( \frac{z}{e} \right)}^z \left( 1 + O \left( \frac{1}{z} \right) \right)}
のように書ける。この漸近展開のより進んだ応用は Re(z ) が定数の複素変数 z に対してである。例えば直線 1/4 + it 上でリーマン・ジーゲルテータ関数 の Im(z ) において適用されたスターリングの公式を見よ。
ビネーの公式 [ ]
スターリングの公式は収束しない級数を伴うので解析的に扱いづらいが、収束しない級数を収束する積分に換えたものとしてビネーの(第二)公式がある。
Γ
(
z
+
1
)
=
2
π
z
(
z
e
)
z
e
μ
(
z
)
,
μ
(
z
)
=
2
∫
0
∞
arctan
(
t
/
z
)
e
2
π
t
−
1
d
t
,
(
ℜ
z
>
0
)
{\displaystyle \Gamma(z+1)=\sqrt{2\pi z}\left(\frac{z}{e}\right)^ze^{\mu(z)},\quad
\mu(z)=2\int_0^\infty \frac{\arctan(t/z)}{e^{2\pi t}-1}dt,\quad(\real{z}>0)}
ビネーの公式は、スターリングの級数を形式的に(収束条件を無視して)操作することによっても導かれるが、厳密には対数ガンマ関数の導関数にアベル・プラナの和公式 を適用して得られる。
log
Γ
(
z
)
=
lim
N
→
∞
z
log
N
+
∑
n
=
1
N
log
n
−
∑
n
=
0
N
log
(
n
+
z
)
{\displaystyle \log\Gamma(z)=\lim_{N\to\infty}z\log{N}+\sum_{n=1}^{N}\log{n}-\sum_{n=0}^{N}\log(n+z)}
d
d
z
log
Γ
(
z
)
=
lim
N
→
∞
log
N
−
∑
n
=
0
N
1
n
+
z
{\displaystyle \frac{d}{dz}\log\Gamma(z)=\lim_{N\to\infty}\log{N}-\sum_{n=0}^{N}\frac{1}{n+z}}
d
2
d
z
2
log
Γ
(
z
)
=
∑
n
=
0
∞
1
(
n
+
z
)
2
{\displaystyle \frac{d^2}{dz^2}\log\Gamma(z)=\sum_{n=0}^{\infty}\frac{1}{(n+z)^2}}
ℜ
z
>
0
{\displaystyle \real{z}>0}
なら
f
=
(
n
+
z
)
−
2
{\displaystyle f=(n+z)^{-2}}
は右半平面において正則であるからプラナの和公式により
d
2
d
z
2
log
Γ
(
z
)
=
∫
0
∞
1
(
n
+
z
)
2
d
t
+
1
2
z
2
+
i
∫
0
∞
(
z
+
i
t
)
−
2
−
(
z
−
i
t
)
−
2
e
2
π
t
−
1
d
t
=
1
z
+
1
2
z
2
+
i
∫
0
∞
(
z
+
i
t
)
−
2
−
(
z
−
i
t
)
−
2
e
2
π
t
−
1
d
t
{\displaystyle \begin{align}\frac{d^2}{dz^2}\log\Gamma(z)
&=\int_0^\infty \frac{1}{(n+z)^2}dt+\frac{1}{2z^2}+i\int_0^\infty \frac{(z+it)^{-2}-(z-it)^{-2}}{e^{2\pi t}-1}dt\\
&=\frac{1}{z}+\frac{1}{2z^2}+i\int_0^\infty \frac{(z+it)^{-2}-(z-it)^{-2}}{e^{2\pi t}-1}dt\\
\end{align}}
積分して
d
d
z
log
Γ
(
z
)
=
C
1
+
log
z
−
1
2
z
+
i
∫
0
∞
−
(
z
+
i
t
)
−
1
+
(
z
−
i
t
)
−
1
e
2
π
t
−
1
d
t
{\displaystyle \frac{d}{dz}\log\Gamma(z)=C_1+\log{z}-\frac{1}{2z}
+i\int_0^\infty \frac{-(z+it)^{-1}+(z-it)^{-1}}{e^{2\pi t}-1}dt}
log
Γ
(
z
)
=
C
2
+
C
1
z
+
z
log
z
−
log
z
2
+
i
∫
0
∞
−
log
(
z
+
i
t
)
+
log
(
z
−
i
t
)
e
2
π
t
−
1
d
t
=
C
2
+
C
1
z
+
(
z
−
1
2
)
log
z
+
2
∫
0
∞
arctan
(
t
/
z
)
e
2
π
t
−
1
d
t
{\displaystyle \begin{align}\log\Gamma(z)
&=C_2+C_1z+z\log z-\frac{\log z}{2}+i\int_0^\infty \frac{-\log{(z+it)}+\log{(z-it)}}{e^{2\pi t}-1}dt\\
&=C_2+C_1z+\left(z-\frac{1}{2}\right)\log z+2\int_0^\infty \frac{\arctan(t/z)}{e^{2\pi t}-1}dt\\
\end{align}}
ℜ
z
>
0
{\displaystyle \real{z}>0}
なら
|
tan
(
t
/
z
)
|
<
M
{\displaystyle |\tan(t/z)|<M}
は有界であるから
|
∫
0
∞
arctan
(
t
/
z
)
e
2
π
t
−
1
d
t
|
≤
|
∫
0
|
z
|
1
/
2
arctan
(
t
/
z
)
e
2
π
t
−
1
d
t
|
+
|
∫
|
z
|
1
/
2
∞
arctan
(
t
/
z
)
e
2
π
t
−
1
d
t
|
≤
∫
0
|
z
|
1
/
2
∑
k
=
1
∞
|
t
/
z
|
k
d
t
2
π
t
+
M
∫
|
z
|
1
/
2
∞
e
2
π
d
t
(
e
2
π
−
1
)
e
t
e
2
π
(
|
z
|
>
1
)
≤
∫
0
|
z
|
1
/
2
d
t
2
π
(
|
z
|
−
t
)
+
M
e
2
π
e
2
π
−
1
∫
t
=
|
z
|
1
/
2
∞
e
−
2
π
t
d
t
{\displaystyle \begin{align}
\left|\int_0^\infty \frac{\arctan(t/z)}{e^{2\pi t}-1}dt\right|
&\le\left|\int_0^{|z|^{1/2}}\frac{\arctan(t/z)}{e^{2\pi t}-1}dt\right|
+\left|\int_{|z|^{1/2}}^\infty \frac{\arctan(t/z)}{e^{2\pi t}-1}dt\right|\\
&\le\int_0^{|z|^{1/2}}\frac{\sum_{k=1}^\infty |t/z|^kdt}{2\pi{t}}
+M\int_{|z|^{1/2}}^\infty \frac{e^{2\pi}dt}{(e^{2\pi}-1)e^{t}e^{2\pi}}\qquad(|z|>1)\\
&\le\int_0^{|z|^{1/2}}\frac{dt}{2\pi(|z|-t)}+\frac{Me^{2\pi}}{e^{2\pi}-1}\int_{t=|z|^{1/2}}^\infty e^{-2\pi t}dt\\
\end{align}}
lim
z
→
∞
|
∫
0
∞
arctan
(
t
/
z
)
e
2
π
t
−
1
d
t
|
≤
lim
z
→
∞
1
2
π
[
log
(
|
z
|
−
t
)
]
0
|
z
|
1
/
2
+
lim
z
→
∞
M
e
2
π
2
π
(
e
2
π
−
1
)
[
−
e
−
2
π
t
]
|
z
|
1
/
2
∞
=
0
{\displaystyle \lim_{z\to\infty}\left|\int_0^\infty \frac{\arctan(t/z)}{e^{2\pi t}-1}dt\right|
\le\lim_{z\to\infty}\frac{1}{2\pi}\Big[\log\left(|z|-t\right)\Big]_0^{|z|^{1/2}}
+\lim_{z\to\infty}\frac{Me^{2\pi}}{2\pi\left(e^{2\pi}-1\right)}\Big[-e^{-2\pi t}\Big]_{|z|^{1/2}}^\infty
=0}
である。
スターリングの公式と比較して積分定数を求め
log
Γ
(
z
)
=
1
2
log
2
π
−
z
+
(
z
−
1
2
)
log
z
+
2
∫
0
∞
arctan
(
t
/
z
)
e
2
π
t
−
1
d
t
{\displaystyle \log\Gamma(z)=\frac{1}{2}\log2\pi-z+\left(z-\frac{1}{2}\right)\log{z}
+2\int_0^\infty \frac{\arctan(t/z)}{e^{2\pi t}-1}dt}
真数に直して
Γ
(
z
+
1
)
=
z
Γ
(
z
)
=
2
π
z
(
z
e
)
z
e
μ
(
z
)
,
μ
(
z
)
=
2
∫
0
∞
arctan
(
t
/
z
)
e
2
π
t
−
1
d
t
,
(
ℜ
z
>
0
)
{\displaystyle \Gamma(z+1)=z\Gamma(z)=\sqrt{2\pi z}\left(\frac{z}{e}\right)^ze^{\mu(z)},\quad
\mu(z)=2\int_0^\infty \frac{\arctan(t/z)}{e^{2\pi t}-1}dt,\quad(\real{z}>0)}
を得る。なお、ビネーの公式を元にして部分積分を繰り返すとスターリングの級数が得られる。
収束級数形式のスターリングの公式 [ ]
トーマス・ベイズの John Canton への書簡が1763年に王立協会により公表されている。それによると、スターリングの公式は収束級数 ではないとされていた[2] 。
スターリングの公式の収束級数形式を得るには以下を評価する。
∫
0
∞
2
arctan
(
t
/
z
)
exp
(
2
π
t
)
−
1
d
t
=
ln
Γ
(
z
)
−
(
z
−
1
2
)
ln
z
+
z
−
1
2
ln
(
2
π
)
{\displaystyle \int_0^\infty \frac{2\arctan(t/z)}{\exp(2\pi t)-1}\, dt
=\ln\Gamma (z) - \left( z-\frac12 \right) \ln z +z - \frac12\ln(2\pi)}
1つの方法として、階乗冪 の逆数の収束級数を使う方法がある。
z
n
¯
=
z
(
z
+
1
)
⋯
(
z
+
n
−
1
)
{\displaystyle z^{\overline n} = z(z+1) \cdots (z+n-1)}
としたとき、次のようになる。
∫
0
∞
2
arctan
(
t
/
z
)
exp
(
2
π
t
)
−
1
d
t
=
∑
n
=
1
∞
c
n
(
z
+
1
)
n
¯
{\displaystyle \int_0^\infty \frac{2\arctan(t/z)}{\exp(2\pi t)-1} \, dt
=\sum_{n=1}^\infty \frac{c_n}{(z+1)^{\overline n}}}
ここで
c
n
=
1
n
∫
0
1
x
n
¯
(
x
−
1
2
)
d
x
{\displaystyle c_n =\frac{1}{n} \int_0^1 x^{\overline n} \left( x-\frac{1}{2} \right)\, dx}
である。以上から次のようなスターリング級数が得られる。
ln
Γ
(
z
)
=
(
z
−
1
2
)
ln
z
−
z
+
1
2
ln
(
2
π
)
+
1
12
(
z
+
1
)
+
1
12
(
z
+
1
)
(
z
+
2
)
+
59
360
(
z
+
1
)
(
z
+
2
)
(
z
+
3
)
+
29
60
(
z
+
1
)
(
z
+
2
)
(
z
+
3
)
(
z
+
4
)
+
⋯
{\displaystyle \begin{align}
\ln \Gamma (z) &= \left( z-\frac{1}{2} \right) \ln z -z + \frac{1}{2}\ln(2\pi) \\
&\quad +\frac{1}{12(z+1)} + \frac{1}{12(z+1)(z+2)} + \frac{59}{360(z+1)(z+2)(z+3)} + \frac{29}{60(z+1)(z+2)(z+3)(z+4)} + \cdots
\end{align} }
これは、
ℜ
(
z
)
>
0
{\displaystyle \Re(z)>0}
のとき収束する。
計算機向けの変形 [ ]
ガンマ関数の(関数電卓などの)計算機向けの近似として次の式がある。
Γ
(
z
)
≈
2
π
z
(
z
e
z
sinh
1
z
+
1
810
z
6
)
z
,
{\displaystyle \Gamma(z) \approx \sqrt{\frac{2 \pi}{z} } \left( \frac{z}{e} \sqrt{ z \sinh \frac{1}{z} + \frac{1}{810z^6} } \right)^{z},}
これは、次と同等である。
2
ln
Γ
(
z
)
≈
ln
(
2
π
)
−
ln
z
+
z
(
2
ln
z
+
ln
(
z
sinh
1
z
+
1
810
z
6
)
−
2
)
,
{\displaystyle 2 \ln \Gamma(z) \approx \ln(2 \pi) - \ln z + z \left(2 \ln z + \ln \left( z \sinh \frac{1}{z} + \frac{1}{810z^6} \right) - 2 \right),}
これらはスターリングの公式を組み替えて、その結果生じる冪級数と双曲線正弦 関数のテイラー展開 の間の合致を観察することで得られる。この近似は z の実数部が 8 以上のとき、小数点以下 8 桁を超える精度を持つ。2002年、Robert H. Windschitl がリソースの制限された計算機(電卓など)でのそれなりの正確性を持った近似としてこれを示した(参考文献の 'Toth' 参照)。
Gergő Nemes は 2007年にほぼ同程度の結果を与える近似式を提案した。こちらはより単純である。
Γ
(
z
)
≈
2
π
z
(
1
e
(
z
+
1
12
z
−
1
10
z
)
)
z
,
{\displaystyle \Gamma(z) \approx \sqrt{\frac{2 \pi}{z} } \left( \frac{1}{e} \left( z + \frac{1}{12z- \frac{1}{10z}} \right) \right)^{z},}
これは、次と同等である。
ln
Γ
(
z
)
≈
1
2
(
ln
(
2
π
)
−
ln
z
)
+
z
(
ln
(
z
+
1
12
z
−
1
10
z
)
−
1
)
.
{\displaystyle \ln \Gamma(z) \approx \frac{1}{2} \left(\ln(2 \pi) - \ln z \right)
+ z \left(\ln \left( z + \frac{1}{12z- \frac{1}{10z}} \right)-1\right). }
歴史 [ ]
この公式を最初に発見したのはアブラーム・ド・モアブルで、次のような形式だった。
n
!
∼
[
c
o
n
s
t
a
n
t
]
⋅
n
n
+
1
/
2
e
−
n
{\displaystyle n!\sim [{\rm constant}]\cdot n^{n+1/2} e^{-n}}
スターリングの貢献は定数が
2
π
{\displaystyle \sqrt{2\pi}}
であることを示したことである。より正確な形式は Jacques Binet が見出した。
スターリングの近似の「一次」バージョン
n
!
=
n
n
{\displaystyle n!=n^n }
は、マックス・プランクが1901年の黒体 放射の論文で使用した。これは多量の光子 や振動子についての黒体放射エネルギーの方程式にリンクしている。この近似は量子論 でよく使われ、例えばピーター・デバイとルイ・ド・ブロイも使っている。アルベルト・アインシュタインとサティエンドラ・ボースは違う方式を採用した。非常に大きな n について確率分布をグラフに描画してみると、両者はほぼ平行になる。
連続変数への補間 [ ]
ガンマ函数とパイ函数 [ ]
階乗函数は負の整数を除く任意の実数に対するものに一般化することができる。例えば * 0! = 1! = 1 , * (−1/2)! = √ π , * (1/2)! = √ π / 2 .
負の整数を除けば、階乗函数は非整数の値に対しても定義することができるが、そのためには解析学 の道具立てが必要である。そのように階乗の値を「補間」して得られるものの一つがガンマ函数 Γ(z ) である(ただし引数が 1 だけずれる)。これは負の整数を除く任意の複素数 z に対して定義される。z の実部が正である場合には
Γ
(
z
)
=
∫
0
∞
t
z
−
1
e
−
t
d
t
{\displaystyle \Gamma(z)=\int_0^\infty t^{z-1} e^{-t}\, dt}
で与えられる。ガンマ函数と階乗との関係は、任意の自然数 n に対して
n
!
=
Γ
(
n
+
1
)
{\displaystyle n!=\Gamma(n+1)}
が成り立つことである。オイラーのもともとの定義式は
Γ
(
z
)
=
lim
n
→
∞
n
z
n
!
∏
k
=
0
n
(
z
+
k
)
{\displaystyle \Gamma(z)=\lim_{n\to\infty}\frac{n^zn!}{\displaystyle\prod_{k=0}^n (z+k)}}
である。ガウス の導入した別表記として、負でない実数 z に対するパイ函数 Π(z ) は
Π
(
z
)
=
∫
0
∞
t
z
e
−
t
d
t
{\displaystyle \Pi(z)=\int_0^\infty t^{z} e^{-t}\,dt}
を満たす。ガンマ函数との関係は
Π
(
z
)
=
Γ
(
z
+
1
)
{\displaystyle \Pi(z) = \Gamma(z+1)}
である。
Π
(
n
)
=
n
!
for
n
∈
N
{\displaystyle \Pi(n) = n!\text{ for }n \in \mathbf{N}}
が成り立つことを思えば、こちらのほうが階乗を補完した函数としては適していると言えるかもしれない。さてパイ函数は階乗が満たすのと同じ漸化式
Π
(
z
)
=
z
Π
(
z
−
1
)
{\displaystyle \Pi(z) = z\Pi(z-1)}
を、しかし定義される限り任意の複素数 z に対して満たす。事実としてはこれはもう漸化式ではなくて函数等式 と見るべきものであるが。この函数等式をガンマ函数に関するものに書き換えれば
Γ
(
n
+
1
)
=
n
Γ
(
n
)
{\displaystyle \Gamma(n+1)=n\Gamma(n)}
となる。階乗を延長したものがパイ函数なのだから、定義可能な任意の複素数 z に対して
z
!
:=
Π
(
z
)
{\displaystyle z! := \Pi(z)}
と定めることは可能である。これらの補間函数を用いて半整数 における階乗の値をさだめるならば、例えば
Γ
(
1
2
)
=
(
−
1
2
)
!
=
Π
(
−
1
2
)
=
π
{\displaystyle \Gamma\left (\frac{1}{2}\right )=\left (-\frac{1}{2}\right )!=\Pi\left (-\frac{1}{2}\right ) = \sqrt{\pi}}
が成り立ち、さらに自然数 n ∈ N に対して
Γ
(
1
2
+
n
)
=
(
−
1
2
+
n
)
!
=
Π
(
−
1
2
+
n
)
=
π
∏
k
=
1
n
2
k
−
1
2
=
(
2
n
)
!
4
n
n
!
π
=
(
2
n
−
1
)
!
2
2
n
−
1
(
n
−
1
)
!
π
{\displaystyle \Gamma\left (\frac{1}{2}+n\right ) = \left (-\frac{1}{2}+n\right )! = \Pi\left (-\frac{1}{2}+n\right ) = \sqrt{\pi} \prod_{k=1}^n {2k - 1 \over 2} = {(2n)! \over 4^n n!} \sqrt{\pi} = {(2n-1)! \over 2^{2n-1}(n-1)!} \sqrt{\pi}}
が得られる。例えば
Γ
(
4.5
)
=
3.5
!
=
Π
(
3.5
)
=
1
2
⋅
3
2
⋅
5
2
⋅
7
2
π
=
8
!
4
4
4
!
π
=
7
!
2
7
3
!
π
=
105
16
π
≈
11.63.
{\displaystyle \Gamma\left (4.5 \right ) = 3.5! = \Pi\left (3.5\right ) = {1\over 2}\cdot{3\over 2}\cdot{5\over 2}\cdot{7\over 2} \sqrt{\pi} = {8! \over 4^4 4!} \sqrt{\pi} = {7! \over 2^7 3!} \sqrt{\pi} = {105 \over 16} \sqrt{\pi} \approx 11.63.}
同様に n ∈ N に対して
Γ
(
1
2
−
n
)
=
(
−
1
2
−
n
)
!
=
Π
(
−
1
2
−
n
)
=
π
∏
k
=
1
n
2
1
−
2
k
=
(
−
4
)
n
n
!
(
2
n
)
!
π
{\displaystyle \Gamma\left (\frac{1}{2}-n\right ) = \left (-\frac{1}{2}-n\right )! = \Pi\left (-\frac{1}{2}-n\right ) = \sqrt{\pi} \prod_{k=1}^n {2 \over 1 - 2k} = {(-4)^n n! \over (2n)!} \sqrt{\pi}}
が成り立ち、例えば
Γ
(
−
2.5
)
=
(
−
3.5
)
!
=
Π
(
−
3.5
)
=
2
−
1
⋅
2
−
3
⋅
2
−
5
π
=
(
−
4
)
3
3
!
6
!
π
=
−
8
15
π
≈
−
0.9453.
{\displaystyle \Gamma\left (-2.5 \right ) = (-3.5)! = \Pi\left (-3.5\right ) = {2\over -1}\cdot{2\over -3}\cdot{2\over -5} \sqrt{\pi} = {(-4)^3 3! \over 6!} \sqrt{\pi} = -{8 \over 15} \sqrt{\pi} \approx -0.9453.}
パイ函数が殆ど全ての複素数値に対して定義される階乗の延長として唯一のものであるかは確かなことではない、それは定義される限り解析的 としても同じことである。しかし、ふつうはこれが階乗の複素函数への最も自然な延長であるものと考える。例えば、ボーア・モレルップの定理 はガンマ函数が Γ(1) = 1 かつ函数等式 Γ(n + 1) = n Γ(n ) を満足する、ガウス平面の全域で有理型 かつ実軸の正の部分で対数凸 となるような唯一の函数であることを述べる。同様の主張はパイ函数に関しても、函数等式 Π(n ) = n Π(n − 1) に関して述べられる。
そうは言うものの、解析的函数論の意味で恐らくより簡明な、階乗の値を補間する複素函数は存在する。例えばアダマール の「ガンマ」函数[21] はガンマ函数とは異なり整函数 になる[22] 。
オイラーはまた非整数の階乗に対する近似無限乗積
n
!
=
Π
(
n
)
=
∏
k
=
1
∞
(
k
+
1
k
)
n
k
n
+
k
=
[
(
2
1
)
n
1
n
+
1
]
[
(
3
2
)
n
2
n
+
2
]
[
(
4
3
)
n
3
n
+
3
]
⋯
.
{\displaystyle \begin{align}n! = \Pi(n) &= \prod_{k = 1}^\infty \left(\frac{k+1}{k}\right)^n\!\!\frac{k}{n+k} \\ &= \left[ \left(\frac{2}{1}\right)^n\frac{1}{n+1}\right]\left[ \left(\frac{3}{2}\right)^n\frac{2}{n+2}\right]\left[ \left(\frac{4}{3}\right)^n\frac{3}{n+3}\right]\cdots. \end{align}}
についても考察している。これは上記のガンマ函数に関する公式と同じものと見做すことができる。しかしこの公式は収束が遅く、実用的な意味でパイ函数やガンマ函数の値を計算することに利用することはできない。
ガウス平面上での挙動 [ ]
複素変数に対する階乗の絶対値と偏角を、単位長さ間隔で −3 ≤ x ≤ 3 , −2 ≤ y ≤ 2 の範囲で描いた等高線。太くなぞった等高線は φ = ±π である。
複素変数の階乗の値をガンマ函数による表現を通して評価することができる。絶対値 ρ と偏角 φ を用いて
f
=
ρ
exp
(
i
φ
)
=
(
x
+
i
y
)
!
=
Γ
(
x
+
i
y
+
1
)
{\displaystyle f=\rho \exp(i\varphi)=(x+iy)!=\Gamma(x+iy+1)}
と書けば、絶対値一定曲線 ρ = (定数) と偏角一定曲線 φ = (定数) を等値線として格子を描くことができる。一定間隔で引いた等値線の間にさらに細かく等値線を引けば、それが補間で得られる値である。極である負の整数においては絶対値と偏角が定義できず、またその周辺で等値線は密になる。
展開の係数の最初の方
n
gn
近似値
0
1
1
1
−γ
−0.5772156649
2
π
2
12
+
γ
2
2
{\displaystyle \frac{\pi^2}{12}+\frac{\gamma^2}{2}}
0.9890559955
3
−
ζ
(
3
)
3
−
π
2
γ
12
−
γ
3
6
{\displaystyle -\frac{\zeta(3)}{3}-\frac{\pi^2\gamma}{12}-\frac{\gamma^3}{6}}
−0.9074790760
γ はオイラー・マスケローニ定数 、ζ はリーマンゼータ函数 である。
|z | < 1 に対してはテイラー展開
z
!
=
∑
n
=
0
∞
g
n
z
n
{\displaystyle z!=\sum_{n=0}^{\infty} g_n z^n}
が利用できる。この展開のより多くの項は、Sage のような計算機代数システム で計算できる。
階乗の近似 [ ]
展開の係数 a n [23]
n
an
0
1 ⁄12
1
1 ⁄30
2
53 ⁄210
3
195 ⁄371
4
22999 ⁄22737
5
29944523 ⁄19733142
6
109535241009 ⁄48264275462
大きな値に対する階乗の値の近似をディガンマ函数 の積分を通じて連分数 表示を用いて記述できる。この方法はスティルチェス による[24] もので、z ! = exp(P (z )) と書けば P (z ) は
P
(
z
)
=
p
(
z
)
+
log
(
2
π
)
/
2
−
z
+
(
z
+
1
2
)
log
(
z
)
{\displaystyle P(z) = p(z) + \log(2\pi)/2 - z + \left(z+\frac{1}{2}\right)\log(z)}
で、スティルチェスはこの第一項 p (z ) の連分数展開
p
(
z
)
=
a
0
z
+
a
1
z
+
a
2
z
+
a
3
z
+
⋱
{\displaystyle
p(z)=\cfrac{a_0}{z+
\cfrac{a_1}{z+
\cfrac{a_2}{z+
\cfrac{a_3}{z+\ddots}}}}
}
を与えた。
さて、任意の複素数 z ≠ 0 に対して log(z !) = P (z ) あるいは log(Γ(z + 1)) = P (z ) とするのは誤りであり[要出典 ] 、実際には実軸の近くの特定の範囲の z でしか成り立たない(一方 |ℑ(Γ(z + 1))| < π である。引数の実部は大きいほど、虚部はより小さくなければならない。しかし逆の関係式 z ! = exp(P (z )) は原点を除くガウス平面の全域で有効である。ただし実軸の負の部分では収束性は弱くなる[要出典 ] (特異点の周辺ではどのような近似もよい収束性を得ることが難しい)。一方、|ℑ(z )| > 2 または ℜ(z ) > 2 の範囲では上記の六つの係数は double
精度の複素数に対してその階乗の近似値を得るのに十分である。より高い精度でより多くの係数を計算するには rational QD-scheme (H. Rutishauser 's QD algorithm )[25] を用いる。
負の整数に対する拡張不能性 [ ]
関係式 n ! = n × (n − 1)! を使えばある整数に対する階乗をそれより「小さい」整数の階乗から計算できる。この関係式を逆に使えば、「大きい」整数に対して与えられた階乗から
(
n
−
1
)
!
=
n
!
n
{\displaystyle (n-1)! = \frac{n!}{n}}
と計算することも可能である。しかし注意すべきは、これでは負の整数に関する階乗を計算することはできないということである(この式に従って (−1)! を計算するには零除算が必要となりこれ以下の負の整数における階乗の値の計算は不可能となる)。このことはガンマ函数においても同じことで、ガンマ函数は負の整数を除くガウス平面の全域において定義できるにも拘らず、負の整数における値だけは定義することができない。
一般化 [ ]
多重指数記法 [ ]
多重指数
α
=
(
α
1
,
α
2
,
…
,
α
n
)
{\displaystyle \alpha = (\alpha_1, \alpha_2,\ldots,\alpha_n)}
に対し階乗は、
α
!
=
α
1
!
⋅
α
2
!
⋯
α
n
!
{\displaystyle \alpha ! = \alpha_1! \cdot \alpha_2! \cdots \alpha_n!}
と定義できる。これは例えば、多変数関数の展開に使われる。
デデキント環への拡張 [ ]
マンジュル・バルガヴァ は階乗を一般のデデキント環 上で定義し、いくつかの古典的な問題を解決するために用いた[26] 。それらの階乗は整数ではなく、イデアル となる。
階乗に類似する概念 [ ]
二重階乗の例
(-9)!!
= 1 ⁄105
(-7)!!
= −1 ⁄15
(-5)!!
= 1 ⁄3
(-3)!!
= −1
(-1)!!
= 1
0!!
= 1
1!!
= 1
2!!
= 2
3!!
= 3
4!!
= 8
5!!
= 15
6!!
= 48
7!!
= 105
8!!
= 384
9!!
= 945
10!!
= 3840
11!!
= 10395
12!!
= 46080
13!!
= 135135
14!!
= 645120
15!!
= 2027025
16!!
= 10321920
17!!
= 34459425
18!!
= 185794560
19!!
= 654729075
20!!
= 3715891200
二重階乗 [ ]
階乗の類似として、二重階乗 n !! は自然数 n に対し一つ飛ばしに積を取る。二重階乗 n !! は階乗 n ! の二回反復合成 (n !)! とは異なる。
(
2
n
)
!
!
=
(
2
n
)
(
2
n
−
2
)
⋯
(
2
)
=
2
n
n
!
{\displaystyle (2n)!!=(2n)(2n-2)\cdots(2)=2^n n!}
(
2
n
+
1
)
!
!
=
(
2
n
+
1
)
(
2
n
−
1
)
⋯
(
1
)
=
(
2
n
+
1
)
!
(
2
n
)
!
!
{\displaystyle (2n+1)!!=(2n+1)(2n-1)\cdots(1)=\frac{(2n+1)!}{(2n)!!}}
奇数 n = 1, 3, 5, 7, … に対する二重階乗の最初の方の値は
1, 3, 15, 105, 945, 10395, 135135, …, オンライン整数列大辞典 の数列 A001147
偶数 n = 0, 2, 4, 6, 8, … に対する二重階乗の値の最初の方は
1, 2, 8, 48, 384, 3840, 46080, 645120, … オンライン整数列大辞典 の数列 A000165
で与えられる。
負の奇数にも拡張される(
(
−
(
2
n
+
1
)
)
!
!
=
(
−
1
)
n
(
2
n
−
1
)
!
!
{\displaystyle (-(2n+1))!! = \frac{(-1)^n}{(2n-1)!!}}
)。また、複素数値への拡張として、以下が知られている[27] 。
z
!
!
=
2
[
1
+
2
z
−
cos
(
π
z
)
]
/
4
π
[
cos
(
π
z
)
−
1
]
/
4
Γ
(
1
+
1
2
z
)
{\displaystyle z!!={2}^{[1+2z-\cos(\pi z)]/4}{\pi}^{ [\cos (\pi z)-1 ]/4}\Gamma (1+\frac{1}{2} z)}
多重階乗 [ ]
より一般に多重階乗 (multifactorial ) は、連続した整数の積である通常の階乗 n ! 、一つ飛ばしの積である二重階乗 n !! 、二つ飛ばしの積である三重階乗 n !!! または n !3 、三つ飛ばしの四重階乗 n !!!! または n !4 などを総称して言う。一般の k -重階乗 n !k は正整数 n に関して帰納的に
n
!
k
=
{
n
if
0
≤
n
<
k
;
n
(
(
n
−
k
)
!
k
)
if
n
≥
k
{\displaystyle n!_k= \begin{cases}
n &\text{if }0\le n<k;\\
n\,((n-k)!_k) &\text{if }n\ge k
\end{cases}}
と定義できる。
これとは異なる多重階乗の定義として
z
!
(
k
)
=
z
(
z
−
k
)
⋯
(
k
+
1
)
=
k
(
z
−
1
)
/
k
(
z
k
)
(
z
−
k
k
)
⋯
(
k
+
1
k
)
=
k
(
z
−
1
)
/
k
Γ
(
z
k
+
1
)
Γ
(
1
k
+
1
)
.
{\displaystyle z!^{(k)} = z(z-k)\cdots (k+1)
= k^{(z-1)/k}\left(\frac{z}{k}\right)\left(\frac{z-k}{k}\right)\cdots \left(\frac{k+1}{k}\right)
= k^{(z-1)/k} \frac{\Gamma\left(\frac{z}{k}+1\right)}{\Gamma\left(\frac{1}{k}+1\right)}\,.}
とするものもある。
階乗冪 [ ]
n の k -順列の総数 n k は n から始めて上から k 個の連続する整数の積を取るある意味で不完全な階乗と呼べる階乗の類似物であった。これを下降階乗冪 と呼ぶ。その反対に n から始めて下から k 個の連続する整数の積をとったもの n k を上昇階乗冪 といい、これら二つを総称して階乗冪と呼ぶ。ただし一般に整数 n ではなく実数 x を変数として考えることが多い。
x
n
_
=
∏
k
=
0
n
−
1
(
x
−
k
)
x
n
¯
=
∏
k
=
0
n
−
1
(
x
+
k
)
{\displaystyle \begin{align}
x^{\underline{n}} & =\prod_{k=0}^{n-1}(x-k)\\
x^{\overline{n}} &=\prod_{k=0}^{n-1}(x+k)
\end{align}}
明らかに自然数 n に対して
n
k
_
=
n
!
(
n
−
k
)
!
,
n
k
¯
=
(
n
+
k
−
1
)
!
(
k
−
1
)
!
,
{\displaystyle n^{\underline{k}} = \frac{n!}{(n-k)!},\quad n^{\overline{k}}=\frac{(n+k-1)!}{(k-1)!},}
n
!
=
n
n
_
=
1
n
¯
.
{\displaystyle n!= n^{\underline{n}} = 1^{\overline{n}}.}
また実数 x ≠ 0 に対して
x
0
_
=
x
0
¯
=
1
{\displaystyle x^{\underline{0}}=x^{\overline{0}}=1}
だが x = 0 のときは定義されない。
素数階乗 [ ]
素数階乗 (Primorial ) n # は最初の n -個の素数の総乗
n
#
=
∏
i
=
1
n
p
i
{\displaystyle n\# = \prod_{i=1}^{n} p_i}
である[28] 。オンライン整数列大辞典 の数列 A002110 。
これは、素数が無限に存在するという命題 の証明に用いられることがある。
superfactorial (1) [ ]
Clifford Pickover 版の超階乗 (superfactorial )は、階乗を入れ子に拡張したものである。ドル記号 $を用いて書かれる。
n
$
=
n
!
n
!
=
n
!
n
!
n
!
⋅
⋅
⋅
n
!
⏟
n
!
{\displaystyle n\$={}^{n!}n! = \underbrace{n!^{n!^{\scriptstyle n!^{{\textstyle\,\cdot}^{{\textstyle\,\cdot}^{{\textstyle\,\cdot\,}^{\scriptstyle n!}}}}}}}_{n!}}
nが3以上になると、非常に大きい値になる。
これとは異なる種類の超階乗 (superfactorial) の定義がある(次節参照)。
superfactorial (2) [ ]
ニール・スローンとサイモン・プラウフは、スーパー階乗 (superfactorial )をThe Encyclopedia of Integer Sequences (N.J.A. Sloane 1995 ) の中で定義した。例として、4のスーパー階乗は次のようになる。
s
f
(
4
)
=
1
!
×
2
!
×
3
!
×
4
!
=
288.
{\displaystyle \mathrm{sf}(4)=1! \times 2! \times 3! \times 4!=288. \,}
一般的にスーパー階乗は下の式で定義される。
s
f
(
n
)
=
∏
k
=
1
n
k
!
=
∏
k
=
1
n
k
n
−
k
+
1
=
1
n
⋅
2
n
−
1
⋅
3
n
−
2
⋯
(
n
−
1
)
2
⋅
n
1
.
{\displaystyle
\mathrm{sf}(n)
=\prod_{k=1}^n k! =\prod_{k=1}^n k^{n-k+1}
=1^n\cdot2^{n-1}\cdot3^{n-2}\cdots(n-1)^2\cdot n^1.
}
次のような定義もある。
s
f
(
n
)
=
∏
0
≤
i
<
j
≤
n
(
j
−
i
)
{\displaystyle
\mathrm{sf}(n)
=\prod_{0 \le i < j \le n} (j-i)
}
最初のいくつかの値は、次のようになる[29] 。
1, 1, 2, 12, 288, 34560, 24883200, 125411328000, ……
スーパー階乗は、複素数値にも拡張できる。その結果はバーンズのG関数と呼ばれる。定義は次のようになる。
G
(
z
+
1
)
=
(
2
π
)
z
/
2
exp
(
−
z
+
z
2
(
1
+
γ
)
2
)
∏
k
=
1
∞
{
(
1
+
z
k
)
k
exp
(
z
2
2
k
−
z
)
}
{\displaystyle G(z+1)=(2\pi)^{z/2} \text{exp}\left(- \frac{z+z^2(1+\gamma)}{2} \right) \, \prod_{k=1}^{\infty} \left\{ \left(1+\frac{z}{k}\right)^k \text{exp}\left(\frac{z^2}{2k}-z\right) \right\}}
自然数に対しては、以下が成り立っている。
G
(
n
+
2
)
=
s
f
(
n
)
=
{
0
if
n
=
−
1
,
−
2
,
…
∏
i
=
0
n
i
!
if
n
=
0
,
1
,
2
,
…
{\displaystyle G(n+2)= \mathrm{sf}(n)=\begin{cases} 0&\text{if }n=-1,-2,\dots\\ \prod_{i=0}^{n} i!&\text{if }n=0,1,2,\dots\end{cases}}
hyperfactorial [ ]
ハイパー階乗 (hyperfactorial )は、以下で定義される。
H
(
n
)
=
∏
k
=
1
n
k
k
=
1
1
⋅
2
2
⋅
3
3
⋯
(
n
−
1
)
n
−
1
⋅
n
n
{\displaystyle H(n)= \prod_{k=1}^n k^k= 1^1\cdot2^2\cdot3^3\cdots(n-1)^{n-1}\cdot n^n}
これはとても大きくなっていく。最初のいくつかの値はつぎの通りである[30] 。
1 , 4 , 108 , 27648, 86400000, ……
ハイパー階乗は定義域を複素数にまで拡張できる。それはK函数と呼ばれ、以下で定義される。
K
(
z
)
=
(
2
π
)
(
−
z
+
1
)
/
2
exp
[
(
z
2
)
+
∫
0
z
−
1
ln
(
t
!
)
d
t
]
.
{\displaystyle K(z)=(2\pi)^{(-z+1)/2} \exp\left[\begin{pmatrix} z\\ 2\end{pmatrix}+\int_0^{z-1} \ln(t!)\,dt\right].}
自然数nに対し、次が成り立つ。
K
(
n
+
1
)
=
1
1
2
2
3
3
⋯
n
n
.
{\displaystyle K(n+1)=1^1\, 2^2\, 3^3 \cdots n^n.}
exponential factorial [ ]
exponential factorial [* 3] n ! (感嘆符は右肩に添字として書く)は
n
!
=
{
1
(
n
=
0
)
n
(
n
−
1
)
!
(
n
>
0
)
{\displaystyle n^! = \begin{cases}
1 & (n = 0)\\
n^{{(n-1)}^!} & (n > 0)
\end{cases}}
で与えられる。つまり、自然数 n に対して
n
!
=
n
(
n
−
1
)
(
n
−
2
)
.
.
.
3
2
1
{\displaystyle n^! = n_{}^{(n-1)^{{(n-2)}^{.\,^{.\,^{.\,^{{3}^{{2}^{1}}}}}}}}}
であり、最初の5つの値は次のようになる。
0! = 1, 1! = 1, 2! = 2, 3! = 9, 4! = 262144, …… (オンライン整数列大辞典 の数列 A049384 )
5! の値は十進展開で183231桁にも及ぶきわめて大きな自然数であり、これ以降も非常に大きくなる。
全ての自然数のexponential factorialの逆数の総和は、
∑
n
=
1
∞
1
n
!
=
1.611114925808376736
11
⋯
11
⏟
183213
272243
⋯
{\displaystyle \sum_{n=1}^{\infty} \frac{1}{n^!} =1.611114925808376736\underbrace{11\cdots11}_{183213}272243\cdots}
となる。この数は、超越数 であり、リウヴィル数 である[32] 。
また、高次 exponential factorialが定義される。例として、二次 exponential factorialは、
n
!
!
=
n
!
2
=
n
!
(
n
−
1
)
!
(
n
−
2
)
!
.
.
.
3
!
2
!
1
!
{\displaystyle n^{!!} = n^{{!}^{2}} = {n^{{!}{(n-1)^{{!}{{(n-2)}^{{!}{{.}^{{.}^{{.}^{3^{{!}{2^{{!}{1^!}}}}}}}}}}}}}}}
となる。一般の m -次 exponential factorial は、
n
!
m
=
n
!
(
m
−
1
)
(
n
−
1
)
!
m
=
n
!
(
m
−
1
)
(
n
−
1
)
!
(
m
−
1
)
(
n
−
2
)
!
(
m
1
)
.
.
.
3
!
(
m
−
1
)
2
!
(
m
−
1
)
1
!
(
m
−
1
)
{\displaystyle n^{{!}^{m}} = n^{{{!}^{(m-1)}{(n-1)}^{!^{m}}}} = {n^{{!^{(m-1)}}{(n-1)^{{!^{(m-1)}}{{(n-2)}^{{!^{(m1)}}{{.}^{{.}^{{.}^{3^{{!^{(m-1)}}{2^{{!^{(m-1)}}{1^{!^{(m-1)}}}}}}}}}}}}}}}}}
で与えられる。ただし、n, m は自然数である。
脚注 [ ]
注釈 [ ]
↑ Graham & Knuth Patashnik , p. 111
↑ Biggs , pp. 109−136
↑ Stedman 1677 , pp. 6–9
↑ Stedman 1677 , p. 8
↑ Higgins , p. 12
↑ http://ocw.mit.edu/courses/mathematics/18-01-single-variable-calculus-fall-2006/lecture-notes/lec4.pdf
↑ Wolfram mathworld: Gamma Function
↑ Springer Online Reference Works: Gamma-function
↑ Schmelzer & Trefethen (2007), Computing the Gamma function using contour integrals and rational approximations
↑ Otto Ludwig Hölder, "Über die Eigenschaft der Gammafunction keiner algebraischen Differentialgleichung zu genügen," Math. Ann. , 28 , (1887) pp. 1–13.
↑ Eliakim Hastings Moore, "Concerning transcendentally transcendental functions," Math. Ann. , 48 (1897), pp. 49–74.
↑ A. Ostrowski, "Neuer Beweis der Hölderschen Satzes, dass die Gammafunktion keiner algebraischen Differntialgleichung genügt." Math. Ann. 79 (1919), pp. 286–288.
↑ A. Ostrowski, "Zum Hölderschen Satz über Γ(x). Math. Ann. 94 (1925), pp. 248–251.
↑ E. W. Barnes, "The theory of the Gamma function," Messenger of Math. 29 (1900), pp. 64–128.
↑ F. Hausdorff, "Zum Hölderschen Satz über Γ(x)," Math. Ann. 94 (1925), pp. 244–247.
↑ Raimundas Vidūnas, Expessions for Values of the Gamma Function
↑ Weisstein, Eric W. , "Gamma Function " - MathWorld.(英語)
↑ K.O. Geddes, M.L. Glasser, R.A. Moore and T.C. Scott, Evaluation of Classes of Definite Integrals Involving Elementary Functions via Differentiation of Special Functions , AAECC (Applicable Algebra in Engineering, Communication and Computing), vol. 1, (1990), pp. 149-165, [1]
↑ Guy 2004 , p. 346
↑ Ramanujan 1988 , p. 339
↑ Hadamard 1894
↑ Peter Luschny, Hadamard versus Euler - Who found the better Gamma function? .
↑ Digital Library of Mathematical Functions, http://dlmf.nist.gov/5.10
↑ Hadamard 1894
↑ Peter Luschny, On Stieltjes' Continued Fraction for the Gamma Function. .
↑ The Factorial Function and Generalizations
↑ Weisstein, Eric W. , "Double Factorial " - MathWorld.(英語)
↑ Weisstein, Eric W. , "Primorial " - MathWorld .(英語)
↑ オンライン整数列大辞典 の数列 A000178
↑ オンライン整数列大辞典 の数列 A002109
↑ 巨大数研究 Wiki 指数階乗
↑ Sondow, Jonathan, "Exponential Factorial " - MathWorld. (英語)
↑ The publisher is given as "W.S." who may have been William Smith, possibly acting as agent for the Society of College Youths , to which society the "Dedicatory" is addressed.[3]
↑ このような(n , m )を、ブラウン数 (英: Brown numbers ) と呼ぶ。
↑ 中国語 : 阶幂 (階冪)、指数階乗[31]
出典 [ ]
参考文献 [ ]
ドナルド・E・クヌース、ロナルド・L・グレアム・オーレン・パタシュニク 『コンピュータの数学』 有澤誠・ほか訳、共立出版、1993年8月。ISBN 4-320-02668-3 : 原著 Graham, Ronald L.; Knuth, Donald E.; Patashnik, Oren (1988), Concrete Mathematics]] , Addison-Wesley, Reading MA, ISBN 0-201-14236-8
Keith B. Oldham他 『関数事典(CD-ROM付)』 河村哲也監訳、朝倉書店、2013年12月、ISBN 978-4-254-11136-1 。
Biggs, N. L. (1979), The roots of combinatorics , Historia Math. 6
Stedman, Fabian (1677), Campanalogia , London
Higgins, Peter (2008), Number Story: From Counting to Cryptography , New York: Copernicus, ISBN 978-1-84800-000-1
Guy, Richard K. (2004), “E24 Irrationality sequences” , Unsolved problems in number theory (3rd ed.), Springer-Verlag , ISBN 0-387-20860-7 , Zbl 1058.11001 , http://books.google.com/books?id=1AP2CEGxTkgC&pg=PA346
Ramanujan, Srinivasa (1988), The lost notebook and other unpublished papers , Springer Berlin, ISBN 3-540-18726-X
Hadamard, M. J. (1894) (French), Sur L’Expression Du Produit 1·2·3· · · · ·(n−1) Par Une Fonction Entière , OEuvres de Jacques Hadamard , Centre National de la Recherche Scientifiques, Paris, 1968, http://www.luschny.de/math/factorial/hadamard/HadamardFactorial.pdf
Abramowitz, M. and Stegun, I., Handbook of Mathematical Functions , http://www.math.hkbu.edu.hk/support/aands/toc.htm
Paris, R. B., and Kaminsky, D., Asymptotics and the Mellin-Barnes Integrals , Cambridge University Press, 2001
Whittaker, E. T., and Watson, G. N., A Course in Modern Analysis , fourth edition, Cambridge University Press, 1963. ISBN 0-521-58807-3
Toth, V. T. Programmable Calculators: Calculators and the Gamma Function . http://www.rskey.org/gamma.htm , modified 2006
Elic W. Weisstein, Stirling's Approximation at MathWorld
Stirling's approximation at PlanetMath
テンプレート:仮リンク 『ガンマ関数入門 』 上野健爾 訳・解説、日本評論社〈はじめよう数学6〉、2002年10月25日 。ISBN 4-535-60846-6 。
小松勇作 『特殊函数 』 朝倉書店〈近代数学講座5〉、2004年3月15日 、復刊。ISBN 978-4-254-11655-7 。
寺沢寛一 『自然科学者のための数学概論 』 岩波書店、1983年5月18日 、増訂版。ISBN 978-4-00-005480-5 。
Milton Abramowitz; Irene A. Stegun, ed (1965-06-01). Handbook of Mathematical Functions: with Formulas, Graphs, and Mathematical Tables . Dover Books on Mathematics. Dover Publications. ISBN 0-486-61272-4 .
E. T. Whittaker ; G. N. Watson (1996-09-13) [1927]. A Course of Modern Analysis . Cambridge Mathematical Library (4th ed.). Cambridge University Press. ISBN 0-521-58807-3 .
M. Abramowitz and I. A. Stegun, eds. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables . New York: Dover, 1972. (See Chapter 6.)
G. Arfken and H. Weber. Mathematical Methods for Physicists . Harcourt/Academic Press, 2000. (See Chapter 10.)
W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling. Numerical Recipes in C . Cambridge, UK: Cambridge University Press, 1988. (See Section 6.2.)
Milton Abramowitz and Irene A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables . Dover (1965) ISBN 978-0486612720
E. T. Whittaker and G. N. Watson, A Course of Modern Analysis . Cambridge University Press (1927; reprinted 1996) ISBN 978-0521588072
George B. Arfken and Hans J. Weber, Mathematical Methods for Physicists , Academic Press; ジョージ.ブラウン・アルフケン、ハンス.J・ウェーバー (著)、 権平健一郎 、神原武志 、小山直人 (翻訳) 『基礎物理数学第4版Vol.3 特殊関数』 講談社 (2001) ISBN 978-4061539792
関連項目 [ ]
巨大数
総乗
三角数 : 階乗の加法的な対応物
関数一覧
特殊関数
複素解析
ベータ関数
リーマンゼータ関数
分数階微積分学
外部リンク [ ]
Hazewinkel, Michiel, ed. (2001), “Factorial” , Encyclopaedia of Mathematics , Springer , ISBN 978-1556080104 , http://eom.springer.de/p/f038080.htm
Weisstein, Eric W. , "Factorial " - MathWorld.(英語)
この記事は、クリエイティブ・コモンズ・ライセンス 表示-継承 3.0 非移植 のもと提供されているオンライン数学辞典『PlanetMath』の項目Factorial の本文を含む
Weisstein, Eric W. , "Gamma Function " - MathWorld.(英語)
ガンマ関数とベータ関数 (PDF ) (日本語)
http://factorielle.free.fr
階乗の計算
A collected derivation of Stirling's equation by Timothy Jones