ガンマ関数

とは、数学において階乗の概念を複素数全体に拡張した特殊関数。複素階乗とも。一般に \Gamma (z) と表記される。

自然数 n に対しては、ガンマ関数と n の階乗との間では次の関係式が成り立つ：
:	n!=\Gamma (n+1),\ \Gamma (n)=(n-1)!.

## 歴史
階乗を整数以外に拡張する試みは1720年代にベルヌーイとゴールドバッハにより検討され、1729年10月6日のベルヌーイからゴールドバッハへの手紙において、ベルヌーイは総乗記号を用いて記述することに成功している。
x!=\lim_{n\to\infty}\left(n+1+\frac{x}{2}\right)^{x-1}\prod_{k=1}^{n}\frac{k+1}{k+x}
これは、負の整数以外の実数に対し明確に定義を与えている。

1729年にレオンハルト・オイラーは2つの異なる記述で定義を与えた。1つは負の整数以外のすべての複素数nに対して無限積を用いることで定義される。この式は1729年10月13日付のゴールドバッハ宛の手紙に見られる。
n! = \prod_{k=1}^\infty \frac{\left(1+\frac{1}{k}\right)^n}{1+\frac{n}{k}}\,,
さらに1730年1月8日のゴールドバッハ宛の手紙には、積分での表現を発見したことを記述している。
n!=\int_0^1 (-\log s)^n\, ds\,,
オイラーはさらに反射公式を含む、この関数の重要な機能的性質のいくつかを発見した。

カール・フリードリヒ・ガウスはオイラーの総積公式を次のように記述しなおした。
\Pi(z-1) = \lim_{m\to\infty}\frac{m^z m!}{z(z+1)(z+2)\cdots(z+m)}
ガウスはこの公式を用いて乗法定理を証明し、楕円積分に関する研究につなげた。

「ガンマ関数」という名称と記号\Gamma は、1814年にルジャンドルが導入したものであり、ルジャンドルはオイラーの積分による定義を現代的な形に書き直した。ガウスによるπ（総乗記号）を用いた記述は古い文献にしばしば見られるが、現代の文献ではルジャンドルの表記が主流である（ただし、\Pi(z) = \Gamma(z+1) である）。

カール・ワイエルシュトラスはさらに別の記述を開発し、複素解析におけるガンマ関数の役割を確立した。
\Gamma(z) = \frac{e^{-\gamma z}}{z} \prod_{k=1}^\infty \left(1 + \frac{z}{k}\right)^{-1} e^\frac{z}{k},ここでγはオイラー・マスケローニ定数である。

フィリップ・J・デイヴィスは1963年のショーヴネ賞を受賞した論文でこの分野が歴史的に発展してきたこと、また今後も発展の余地が多いことを記述している。「それぞれの世代がガンマ関数について興味深いことを発見してきました。おそらく次の世代もそうすうるでしょう。」

## 定義
実部が正となる複素数 z に対して、次の広義積分で定義される複素関数：
:\Gamma(z)=\int^{\infty}_{0} t^{z-1} e^{-t} \,{\rm d}t \qquadただし複素数zの実部は正
を**ガンマ関数**と呼ぶ。この積分表示はオイラーの1730年1月8日の積分式を元にしたものであり、t = -\log xで置換積分した記述である。上述のようにガンマ関数の定義はこれ以外にも多数あり、本記事で紹介したものが全てではない。第二種オイラー積分とも呼ばれる。

一般の複素数 z に対しては解析接続もしくは次の極限で定義される。
:\Gamma(z)=\lim_{n\to\infty}\frac{n^zn!}{\prod\limits_{k=0}^{n}{(z+k)}}.
他にも互いに同値となるいくつかの定義が存在する。

## 基本的性質
0 または負の整数でない、かつ実部が正の任意の複素数 z に対して、
:\begin{align}
\Gamma(z+1) &= \int _0 ^{\infty} e^{-t} t^z \, {\rm d}t\\
&= \Bigl[ -e^{-t} t^z \Bigr] _0 ^{\infty} + z\int _0 ^{\infty} e^{-t} t^{z-1} \, {\rm d}t\\
&= z \Gamma(z) \qquad \left( \because \Bigl[ -e^{-t} t^z \Bigr] _0 ^{\infty} = 0 \right) \! ,
\end{align}
となることから、\Gamma(z+1)=z\Gamma(z) が成り立つ。またさらに、
:\begin{align}
\Gamma(1) &= \int _0 ^{\infty} e^{-t} \, {\rm d}t = \Bigl[ -e^{-t} \Bigr] _0 ^{\infty} = \lim _{t\to\infty} \! \left( -e^{-t} +1 \right)\\
&= 1,
\end{align}
である。これらの性質から、任意の正の整数 n に対して、
:\Gamma(n+1)=n\Gamma(n)=n(n-1)\Gamma(n-1)=\cdots =n!\,\Gamma(1)=n!,
より \Gamma (n+1)=n! が成り立つ。その意味でガンマ関数は階乗の定義域を複素平面に拡張したものとなっている。

歴史的には、ガンマ関数は「階乗の複素数への拡張となるもの」（複素階乗）の実例として、オイラーにより考案された。階乗の複素数への拡張となる関数は無数に存在するが「正の実軸上で対数凸である解析関数」という条件を付ければ、それは一意に定まりガンマ関数に他ならない（ボーア・モレルップの定理）。

右半平面においてオイラー積分で定義されたガンマ関数は全平面に有理型に解析接続する。

ガンマ関数は零点を持たず、原点と負の整数に一位の極を持つ。その留数は、
:\operatorname{Res}(\Gamma ,\, -n) = \frac{(-1)^n}{n!}. 
である。

また、1/2 に対するガンマ関数の値は、ガウス積分の結果に一致する。
:\Gamma\!\left(\frac{1}{2}\right) = \sqrt{\pi}.
これより、自然数 n に対して、
:
\Gamma\!\left(\frac{1}{2}+n\right) = \frac{(2n-1)!!}{2^{n}}\sqrt{\pi},

が成立することがわかる。ここで !! は二重階乗を表す。この性質を利用して高次元の球の体積と表面積を求めることができる。また、
:
\Gamma\!\left(\frac{1}{2}-n\right) = \frac{ (-2)^n}{(2n-1)!!} \sqrt{\pi}.


## 定義の整合性
定義の積分表示と極限表示が一致することを示す。
{\rm d}t}}
とすれば
:\lim_{n\to\infty}{\left(1-\dfrac{t}{n}\right)^n}=e^{-t}
であるから直感的には
:\lim_{n\to\infty}{G_n(z)}=\int^{\infin}_{0}t^{z-1}e^{-t}{\rm d}t
である。（厳密にははさみうちの原理によって証明される） の置換により
{\rm d}u}}
となる． を除く部分を  として
{\rm d}u=\left[\frac{u^z}{z}\right]_{u=0}^{1}=\frac{1}{z}
g_n(z)=\int_{0}^{1}{\left(\frac{u^{z}}{z}\right)'(1-u)^{n}}{\rm d}u=\frac{n}{z}\int_{u=0}^{1}{u^{z}(1-u)^{n-1}}{\rm d}u=\frac{n}{z}g_{n-1}(z+1)}}
これにより
}}
を得る。故に
}}
である。

## ワイエルシュトラスの乗積表示
オイラーの乗積表示からオイラーの定数
:\gamma=\lim_{n\to\infty}\left(\sum_{k=1}^{n}\frac{1}{k}-\log{n}\right)
を括り出すと**ワイエルシュトラスの乗積表示**が得られる。ワイエルシュトラスはガンマ関数が負の整数に極を持つことを嫌って逆数を用いた。ガンマ関数の逆数は複素平面全体で正則である。
{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^

## ハンケルの積分表示

ガンマ関数は次の周回積分で表される。積分経路は正の無限大から実軸の上側に沿って原点に至り、原点を正の向きに回り、実軸の下側に沿って無限大に戻るものとする。但し、その偏角は-\pi\le\arg(-t)\le\pi,0\le\arg(s)\le2\piとする。

これを**ハンケルの積分表示**と呼ぶ。このハンケルの積分表示は、積分経路を適当に変形し、数値積分でガンマ関数の値を求めるために使われることがある。

### ハンケルの積分表示の導出
極座標表示 (-t)=re^{i\theta} を用いると、実軸の上側に沿う部分は \theta=-\pi で r=\infty から r=\delta まで、原点を回る部分は r=\delta で \theta=-\pi から \theta=\pi まで、実軸の下側に沿う部分は \theta=\pi で r=\delta から r=\infty までとなる。

となり

を得る。

## ルジャンドルの倍数公式
次の恒等式を**ルジャンドルの倍数公式**と呼ぶ。これは**ガウスの乗法公式**の特別な場合である。

### 証明
ベータ関数は以下のように表される。
: \Beta (z_1,z_2)=\frac{\Gamma (z_1)\Gamma (z_2)}{\Gamma (z_1+z_2)}=\int_0^1 t^{z_1-1}(1-t)^{z_2-1} \, dt

ここで z_1=z_2=z とおくと、
: \frac{\Gamma^2(z)}{\Gamma (2z)}=\int_0^1 t^{z-1}(1-t)^{z-1} \, dt

t=\frac{1+x}{2} とおくと
: \frac{\Gamma^2(z)}{\Gamma (2z)}=\frac{1}{2^{2z-1}}\int_{-1}^1 \left(1-x^{2}\right)^{z-1} \, dx
(1-x^2)^{z-1} は偶関数なので
: 2^{2z-1}\Gamma^2(z)=2\Gamma (2z)\int_0^1 (1-x^2)^{z-1} \, dx
ここで
: \Beta \left(\frac{1}{2},z\right)=\int_0^1 t^{\frac{1}{2}-1}(1-t)^{z-1} \, dt, \quad t=s^2
とすると
: \Beta \left(\frac{1}{2},z\right)=2\int_0^1 (1-s^2)^{z-1} \, ds = 2\int_0^1 (1-x^2)^{z-1} \, dx
よって
: 2^{2z-1}\Gamma^2(z)=\Gamma (2z)\Beta \left(\frac{1}{2},z\right)
よって
: \Beta \left(\frac{1}{2},z\right)=\frac{\Gamma \left(\frac{1}{2}\right)\Gamma (z)}{\Gamma \left(z+\frac{1}{2}\right)}, \quad \Gamma \left(\frac{1}{2}\right)=\sqrt{\pi}
よって以下の式が成り立つ。
: \Gamma (z)\Gamma \left(z+\frac{1}{2}\right)=2^{1-2z}\sqrt{\pi} \; \Gamma (2z).

## 乗法公式
次の恒等式を**ガウスの乗法公式**(multiplication formula)という。
{(2\pi)^{(n-1)/2}}\prod_{k=0}^{n-1}{\Gamma{\left(z+\frac{k}{n}\right)}}}}
### 証明
両辺の比をf(z)とすると
}{(2\pi)^{(n-1)/2}\Gamma(nz)}\\
\end{align}
\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についてf(z+m)=f(z)が成立する。スターリングの公式により
{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}}}
途中で

を適用した。

であり、故に
{(2\pi)^{(n-1)/2}}\prod_{k=0}^{n-1}{\Gamma{\left(z+\frac{k}{n}\right)}}}}
が成立する。
## 微分方程式
(x,\ y,\ y_1,\ \ldots ,\ y_n)を変数とする多項式F(x,\ y,\ y_1,\ \ldots ,\ y_n)に対し、

の形で表される微分方程式を代数的微分方程式という。ガンマ関数はいかなる代数的微分方程式も満たさないことが知られている。このことはヘルダーが1887年に最初に証明を与えた後
、、、、ハウスドルフにより、別証明や一般化がなされた。

## いくつかの特殊値

 {3} \approx 2.363\,
\Gamma\left(-\frac{1}{2}\right)\,= -2\sqrt{\pi} \approx -3.545\,
\Gamma\left(\frac{1}{2}\right)\,= \sqrt{\pi} \approx 1.772\,
\Gamma(1)\,=0!=1 \,
\Gamma\left(\frac{3}{2}\right)\,= \frac {\sqrt{\pi}} {2} \approx 0.886\,
\Gamma(2)\,=1!=1 \,
\Gamma\left(\frac{5}{2}\right)\,= \frac {3 \sqrt{\pi}} {4} \approx 1.329\,
\Gamma(3)\,=2!=2 \,
\Gamma\left(\frac{7}{2}\right)\,= \frac {15\sqrt{\pi}} {8} \approx 3.323\,
\Gamma(4)\,=3!=6 \,
}}

## ポリガンマ関数

ガンマ関数の対数微分

:\psi(z)=\frac{d}{dz}\log \Gamma(z)

をディガンマ関数（Digamma function）と呼ぶ。同様の対数微分を繰り返した関数

:\psi^{(n)}(z)=\frac{d^{n+1}}{dz^{n+1}}\log \Gamma(z)

を、ポリガンマ関数（Polygamma function）と呼ぶ。

## 脚注

## 参考文献
*
*
* 柴垣和三雄：「ガンマ函数の理論と応用 : 並びに小数6位まで有效な複素変数のガンマ函数表」、岩波書店 (解析数学叢書)、（1952年）．
*
*
*
*
* 

## 関連項目
*特殊関数
*階乗
*関数一覧
*不完全ガンマ関数
*複素解析
*オイラー積分
*ベータ関数
*リーマンゼータ関数
*分数階微積分学
*Abramowitz and Stegun
*

## 外部リンク
*
*
*