Γ関数

同義語:ガンマ関数第二種Euler積分gamma function

概要

Γ関数(gamma function)とは、$\operatorname{Re}z>0$ で $\Gamma(z)=\int_0^\infty t^{z-1}e^{-t}\,dt$ で定まり、関数等式 $\Gamma(z+1)=z\Gamma(z)$ により階乗を $\Gamma(n+1)=n!$ と延長する関数である。Gauss の極限公式と Weierstrass の乗積表示を持ち、$z=-n$ を留数 $(-1)^n/n!$ の1位の極とする零点のない有理型関数に解析接続される。$(0,\infty)$ 上で正値かつ対数凸で、Bohr–Mollerup の定理により、この性質と関数等式および $\Gamma(1)=1$ で特徴づけられる。$\Gamma(1/2)=\sqrt{\pi}$ や相反公式 $\Gamma(z)\Gamma(1-z)=\pi/\sin\pi z$ も成り立つ。

$$\newcommand{AA}[0]{\mathscr{A}} \newcommand{abs}[1]{\left\lvert#1\right\rvert} \newcommand{Arg}[0]{\operatorname{Arg}} \newcommand{BB}[0]{\mathscr{B}} \newcommand{C}[0]{\mathbb{C}} \newcommand{CC}[0]{\mathscr{C}} \newcommand{floor}[1]{\left\lfloor#1\right\rfloor} \newcommand{ind}[0]{\mathrm{ind}} \newcommand{mmod}[1]{\ \left(\mathrm{mod}\ #1\right)} \newcommand{Mod}[1]{\ \left(\mathrm{mod}\ #1\right)} \newcommand{ord}[0]{\mathrm{ord}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{rank}[0]{\mathrm{rank}} \newcommand{SS}[0]{\mathscr{S}} \newcommand{TT}[0]{\mathscr{T}} \newcommand{UU}[0]{\mathscr{U}} \newcommand{wenvert}[1]{\left\lvert\left\lvert#1\right\rvert\right\rvert} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: 階乗, 広義積分, 部分積分, 正則関数, 凸関数

$\Gamma$ 関数は、階乗を複素数へ拡張する関数である。$\operatorname{Re}z>0$ では積分 $\int_0^\infty t^{z-1}e^{-t}\,dt$ で定まり、非負整数 $n$ に対して $\Gamma(n+1)=n!$ を満たす。さらに解析接続により $0,-1,-2,\ldots$ を除く複素平面上の正則関数になる。Beta関数、Stirlingの公式、ガンマ分布など、解析学や確率論の多くの場面に現れる。たとえば $\Gamma(1)=\int_0^\infty e^{-t}\,dt=1$ である。

定義

以下、複素数 $z$ に対し $\sigma:=\operatorname{Re}z$ とし、$t>0$ に対し $t^{z-1}:=e^{(z-1)\log t}$ とする。$|e^{w}|=e^{\operatorname{Re}w}$ なので $|t^{z-1}|=t^{\sigma-1}$ である。複素数値関数の広義積分は、実部・虚部それぞれの広義積分の和として定める。

$\Gamma$ 関数

$\operatorname{Re}z>0$ を満たす複素数 $z$ に対し、
$$ \Gamma(z):=\int_0^\infty t^{z-1}e^{-t}\,dt $$
と定める(この広義積分は thm-gamma-convergence により絶対収束する)。$\Gamma$ を $\Gamma$関数(gamma function、ガンマ関数)といい、右辺の積分を第二種 Euler 積分という。$\Gamma$ は thm-gamma-continuation により $\mathbb{C}\setminus\{0,-1,-2,\ldots\}$ 上の正則関数に一意に延長され、延長後の関数も同じ記号 $\Gamma$ で表す。

指数減衰による優関数

実数 $a$ に対し、正の定数 $C_a$ が存在して $t\ge1$ で $t^{a}e^{-t/2}\le C_a$ が成り立つ。特に $t\ge1$ で $t^{a}e^{-t}\le C_ae^{-t/2}$ である。

$a\le0$ なら $t\ge1$ で $t^{a}\le1$, $e^{-t/2}\le1$ なので $C_a=1$ でよい。$a>0$ とし、$m\ge a$ なる正の整数 $m$ をとる。$u>0$ に対し $e^{u}=\sum_{k\ge0}u^k/k!\ge u^m/m!$ なので、$u=t/2$ とおいて $e^{-t/2}\le2^mm!\,t^{-m}$ である。よって $t\ge1$ で $t^ae^{-t/2}\le2^mm!\,t^{a-m}\le2^mm!$($a-m\le0$)となり、$C_a=2^mm!$ でよい。後半は $t^ae^{-t}=(t^ae^{-t/2})e^{-t/2}$ による。$\square$

第二種 Euler 積分の収束域

複素数 $z$ について次が成り立つ。

  1. $\sigma=\operatorname{Re}z>0$ ならば $\int_0^\infty|t^{z-1}e^{-t}|\,dt<\infty$ である。したがって $\Gamma(z)$ は絶対収束する広義積分として定まる。
  2. $\sigma\le0$ ならば $\int_0^1|t^{z-1}e^{-t}|\,dt=+\infty$ である。すなわち $\operatorname{Re}z>0$ は絶対収束のための必要十分条件である。

$|t^{z-1}e^{-t}|=t^{\sigma-1}e^{-t}$ は $(0,\infty)$ 上連続なので、任意の有界閉区間上で積分可能である。積分区間を $(0,1]$ と $[1,\infty)$ に分ける。
(1) $\sigma>0$ とする。$0< t\le1$ で $0< t^{\sigma-1}e^{-t}\le t^{\sigma-1}$ であり、$\int_0^1t^{\sigma-1}\,dt=1/\sigma<\infty$ なので、広義積分 の定理「広義積分の比較判定法」(左の端で同じ議論をした形)により $\int_0^1t^{\sigma-1}e^{-t}\,dt$ は収束する。$t\ge1$ では lem-gamma-exponential-bound($a=\sigma-1$)により $t^{\sigma-1}e^{-t}\le C_{\sigma-1}e^{-t/2}$ であり、$\int_1^\infty e^{-t/2}\,dt=2e^{-1/2}<\infty$ なので $\int_1^\infty t^{\sigma-1}e^{-t}\,dt$ も収束する。被積分関数の実部・虚部の絶対値はいずれも $t^{\sigma-1}e^{-t}$ 以下なので、$\Gamma(z)$ は絶対収束する。
(2) $\sigma\le0$ とする。$0< t\le1$ で $e^{-t}\ge e^{-1}$ なので $t^{\sigma-1}e^{-t}\ge e^{-1}t^{\sigma-1}$ であり、$1-\sigma\ge1$ なので $\int_0^1t^{\sigma-1}\,dt=+\infty$ である。同じ比較判定法の発散側により主張が従う。$\square$

Euler の定数の存在

$H_N:=1+\frac12+\cdots+\frac1N$ とおくと、数列 $H_N-\log N$($N\ge1$)は単調減少で $0$ 以上であり、したがって収束する。その極限
$$ \gamma:=\lim_{N\to\infty}\bigl(H_N-\log N\bigr) $$
を Eulerの定数(Euler–Mascheroni 定数)という。$\gamma=0.5772\ldots$ である。

$n\le t\le n+1$ で $\frac1{n+1}\le\frac1t\le\frac1n$ なので、$\frac1{n+1}\le\int_n^{n+1}\frac{dt}{t}=\log(n+1)-\log n\le\frac1n$ である。前者から $(H_{N+1}-\log(N+1))-(H_N-\log N)=\frac1{N+1}-\bigl(\log(N+1)-\log N\bigr)\le0$ で単調減少である。後者を $n=1,\ldots,N-1$ で加えると $\log N\le H_{N-1}\le H_N$ となり、$H_N-\log N\ge0$ である。下に有界な単調減少列は収束する。数値は Apo76 による。$\square$

公式と成立条件

関数等式と階乗

$\operatorname{Re}z>0$ に対し、関数等式 $\Gamma(z+1)=z\,\Gamma(z)$ が成り立つ。また $\Gamma(1)=1$ であり、非負整数 $n$ に対し $\Gamma(n+1)=n!$ である。したがって $\Gamma$ は階乗を $\operatorname{Re}z>0$ へ延長する関数である。

$0<\varepsilon< R<\infty$ に対し、部分積分により
$$ \int_\varepsilon^R t^{z}e^{-t}\,dt=\bigl[-t^{z}e^{-t}\bigr]_\varepsilon^R+z\int_\varepsilon^R t^{z-1}e^{-t}\,dt =\varepsilon^{z}e^{-\varepsilon}-R^{z}e^{-R}+z\int_\varepsilon^R t^{z-1}e^{-t}\,dt $$
である。$|\varepsilon^{z}e^{-\varepsilon}|\le\varepsilon^{\sigma}\to0$($\varepsilon\to0$、$\sigma>0$)、$|R^{z}e^{-R}|=R^{\sigma}e^{-R}\to0$($R\to\infty$、lem-gamma-exponential-bound)であり、両辺の積分は thm-gamma-convergence により収束するので、極限をとって $\Gamma(z+1)=z\Gamma(z)$ を得る。$\Gamma(1)=\int_0^\infty e^{-t}\,dt=1$ であり、$\Gamma(n+1)=n!$ は $n$ に関する帰納法で $\Gamma(n+2)=(n+1)\Gamma(n+1)=(n+1)\,n!$ による。$\square$

指数関数の近似

正の整数 $n$ と $0\le u\le n$ に対し
$$ 0\le e^{-u}-\Bigl(1-\frac un\Bigr)^n\le\frac{u^2}{n}e^{-u} $$
が成り立つ。

$0\le x\le1$ に対し $e^{x}=1+\int_0^xe^{s}\,ds\ge1+x$ かつ $e^{-x}=1-\int_0^xe^{-s}\,ds\ge1-x$ である。$x=u/n$ を代入して $n$ 乗すると $(1-u/n)^n\le e^{-u}$ および $e^{u}\ge(1+u/n)^n$ を得る。前者が左の不等式である。後者から
$$ e^{-u}-\Bigl(1-\frac un\Bigr)^n=e^{-u}\Bigl(1-e^{u}\Bigl(1-\frac un\Bigr)^n\Bigr)\le e^{-u}\Bigl(1-\Bigl(1+\frac un\Bigr)^n\Bigl(1-\frac un\Bigr)^n\Bigr)=e^{-u}\Bigl(1-\Bigl(1-\frac{u^2}{n^2}\Bigr)^n\Bigr) $$
であり、$0\le u^2/n^2\le1$ に対する Bernoulliの不等式(不等式) $(1-u^2/n^2)^n\ge1-u^2/n$ から右の不等式を得る。$\square$

Gauss の極限公式

$\operatorname{Re}z>0$ に対し
$$ \Gamma(z)=\lim_{N\to\infty}\frac{N^{z}\,N!}{z(z+1)\cdots(z+N)} $$
が成り立つ。

$\operatorname{Re}z>0$ を固定し、$g_n(z):=\int_0^1(1-t)^nt^{z-1}\,dt$($n\ge0$)とおく。$g_0(z)=\int_0^1t^{z-1}\,dt=1/z$ であり、$n\ge1$ のとき部分積分により
$$ g_n(z)=\Bigl[\frac{(1-t)^nt^{z}}{z}\Bigr]_0^1+\frac nz\int_0^1(1-t)^{n-1}t^{z}\,dt=\frac nz\,g_{n-1}(z+1) $$
である(境界項は $t=1$ で $0$、$t\to0$ で $|t^{z}|=t^{\sigma}\to0$)。帰納法により
$$ g_N(z)=\frac{N!}{z(z+1)\cdots(z+N)} $$
を得る。よって主張の右辺の第 $N$ 項は $\Pi_N(z):=N^{z}g_N(z)$ である。置換 $u=Nt$ により
$$ \Pi_N(z)=\int_0^1(1-t)^NN^{z}t^{z-1}\,dt=\int_0^N\Bigl(1-\frac uN\Bigr)^Nu^{z-1}\,du $$
となる。$h_N(z):=\int_0^N\bigl(e^{-u}-(1-u/N)^N\bigr)u^{z-1}\,du$ とおくと、lem-gamma-exp-approximation により
$$ |h_N(z)|\le\int_0^N\frac{u^2}{N}e^{-u}u^{\sigma-1}\,du\le\frac1N\int_0^\infty u^{\sigma+1}e^{-u}\,du=\frac{\Gamma(\sigma+2)}{N}\longrightarrow0 $$
である。一方 $\int_0^Ne^{-u}u^{z-1}\,du\to\Gamma(z)$($N\to\infty$)なので、$\Pi_N(z)=\int_0^Ne^{-u}u^{z-1}\,du-h_N(z)\to\Gamma(z)$ である。$\square$

Weierstrass の乗積表示

$\operatorname{Re}z>0$ に対し、無限積
$$ \frac1{\Gamma(z)}=z\,e^{\gamma z}\prod_{n=1}^\infty\Bigl[\Bigl(1+\frac zn\Bigr)e^{-z/n}\Bigr] $$
が成り立つ($\gamma$ は Eulerの定数)。特に $\operatorname{Re}z>0$ で $\Gamma(z)\ne0$ である。

$H_N=\sum_{n=1}^N1/n$ とおく。$N^{z}=e^{z\log N}$ と $z(z+1)\cdots(z+N)=z\,N!\prod_{n=1}^N(1+z/n)$ から、有限の恒等式
$$ \frac{z(z+1)\cdots(z+N)}{N^{z}\,N!}=z\,e^{z(H_N-\log N)}\prod_{n=1}^N\Bigl[\Bigl(1+\frac zn\Bigr)e^{-z/n}\Bigr] $$
が成り立つ。左辺の逆数 $\Pi_N(z):=N^{z}N!/\bigl(z(z+1)\cdots(z+N)\bigr)$ は thm-gamma-gauss-limit により $\Gamma(z)$ に収束し、lem-gamma-euler-constant により $e^{z(H_N-\log N)}\to e^{\gamma z}\ne0$ である。
有限積 $P_N(z):=\prod_{n=1}^N(1+z/n)e^{-z/n}$ が $N\to\infty$ で $0$ でない極限に収束することを示す。$R>0$ をとり $|z|\le R$ とする。$n\ge2R$ なら $w=z/n$ は $|w|\le1/2$ を満たし、主値(対数関数の主値)の対数について $\log(1+w)-w=\sum_{k\ge2}(-1)^{k-1}w^k/k$ なので $|\log(1+w)-w|\le\frac{|w|^2}{2}\sum_{j\ge0}|w|^j\le|w|^2$ である。よって $\sum_{n\ge2R}\bigl(\log(1+z/n)-z/n\bigr)$ は $|z|\le R$ 上一様に絶対収束し($\sum_n R^2/n^2<\infty$)、$N_0:=\lceil2R\rceil$ とおくと $P_N(z)=P_{N_0}(z)\exp\bigl(\sum_{n=N_0+1}^N(\log(1+z/n)-z/n)\bigr)$ は、$P_{N_0}(z)\ne0$($z\notin\{-1,\ldots,-N_0\}$ のとき)である限り $0$ でない極限 $P(z)$ に収束する。$\operatorname{Re}z>0$ ならこの条件は満たされる。
以上により、恒等式の右辺は $0$ でない極限 $z\,e^{\gamma z}P(z)$ を持つ。左辺は $1/\Pi_N(z)$ であり、もし $\Gamma(z)=0$ なら $|1/\Pi_N(z)|\to\infty$ となって右辺の収束に反する。よって $\Gamma(z)\ne0$ であり、$1/\Pi_N(z)\to1/\Gamma(z)$ となるので、恒等式で $N\to\infty$ として主張を得る。$\square$

乗積表示による解析接続

prf-gamma-weierstrass の議論は、無限積 $z e^{\gamma z}\prod_{n\ge1}(1+z/n)e^{-z/n}$ がすべての $z\in\mathbb{C}$ で収束し(有界閉集合上一様)、$\mathbb{C}$ 全体で正則な関数を定めることも示している。この関数の零点は $z=0,-1,-2,\ldots$(いずれも1位)だけである。したがってその逆数は $\mathbb{C}\setminus\{0,-1,-2,\ldots\}$ 上の正則関数であり、$\operatorname{Re}z>0$ で $\Gamma$ と一致するので、thm-gamma-continuation の解析接続と一致する。文献によっては、この乗積表示を $\Gamma$ の定義として採用する(WW96 §12.1)。無限積の収束と正則性に関する一般論は 無限積の定理「正則関数の無限積」にある(Ahl79 Chapter 5 §2.2)。

Euler の乗積表示

$\operatorname{Re}z>0$ に対し
$$ \Gamma(z)=\frac1z\prod_{n=1}^\infty\Bigl[\Bigl(1+\frac1n\Bigr)^{z}\Bigl(1+\frac zn\Bigr)^{-1}\Bigr] $$
が成り立つ。実際、$\prod_{n=1}^{N-1}(1+1/n)=N$ なので
$$ \frac{N^{z}\,N!}{z(z+1)\cdots(z+N)}=\frac1z\prod_{n=1}^{N-1}\Bigl(1+\frac1n\Bigr)^{z}\prod_{n=1}^{N}\Bigl(1+\frac zn\Bigr)^{-1} =\frac1z\Bigl(1+\frac1N\Bigr)^{-z}\prod_{n=1}^{N}\Bigl[\Bigl(1+\frac1n\Bigr)^{z}\Bigl(1+\frac zn\Bigr)^{-1}\Bigr] $$
であり、$(1+1/N)^{-z}\to1$ なので thm-gamma-gauss-limit から従う。歴史的には、Euler がこの形の乗積で $\Gamma$ を導入し、Weierstrass に先立って Newman が Euler の乗積表示から Weierstrass の乗積表示を導いていた(WW96 §12.1)。

解析接続と極・留数

$\Gamma$ は $\mathbb{C}\setminus\{0,-1,-2,\ldots\}$ 上の正則関数へ一意に解析接続される。延長された $\Gamma$ について次が成り立つ。

  1. 正の整数 $m$ と、$\operatorname{Re}z>-m$ かつ $z\notin\{0,-1,\ldots,-(m-1)\}$ なる $z$ に対し $\Gamma(z)=\dfrac{\Gamma(z+m)}{z(z+1)\cdots(z+m-1)}$。
  2. $\Gamma$ は $\mathbb{C}$ 上の有理型関数であり、極(複素解析)は $z=0,-1,-2,\ldots$ に限る。
  3. 各非負整数 $n$ に対し $z=-n$ は1位の極であり、留数は $\operatorname{Res}_{z=-n}\Gamma(z)=\dfrac{(-1)^n}{n!}$ である。
  4. 関数等式 $\Gamma(z+1)=z\Gamma(z)$ は $z\notin\{0,-1,-2,\ldots\}$ で成り立つ。

prop-gamma-holomorphic により $\Gamma$ は $\operatorname{Re}z>0$ で正則である。正の整数 $m$ に対し $\Omega_m:=\{z\mid\operatorname{Re}z>-m\}\setminus\{0,-1,\ldots,-(m-1)\}$、$G_m(z):=\Gamma(z+m)/\bigl(z(z+1)\cdots(z+m-1)\bigr)$ とおく。$z\in\Omega_m$ なら $\operatorname{Re}(z+m)>0$ なので分子は正則、分母は $\Omega_m$ 上で $0$ にならない多項式であり、$G_m$ は $\Omega_m$ 上正則である。
延長であること:$\operatorname{Re}z>0$ で prop-gamma-functional-equation を $z,z+1,\ldots,z+m-1$ に順に適用すると $\Gamma(z+m)=(z+m-1)\cdots z\,\Gamma(z)$ なので $G_m(z)=\Gamma(z)$ である。整合性:$m< m'$, $z\in\Omega_m$ なら、$\operatorname{Re}(z+m)>0$ で同様に $\Gamma(z+m')=(z+m'-1)\cdots(z+m)\Gamma(z+m)$ なので $G_{m'}(z)=G_m(z)$ である。よって $G_m$ たちは $\bigcup_m\Omega_m=\mathbb{C}\setminus\{0,-1,-2,\ldots\}$ 上の一つの正則関数を定め、(1) が従う。(4) は (1) を $z$ と $z+1$ に適用して比較すれば従う。
一意性:$\mathbb{C}\setminus\{0,-1,-2,\ldots\}$ は連結開集合であり、$\{\operatorname{Re}z>0\}$ はその空でない開部分集合なので、一致の定理により延長は一意である。
極と留数:非負整数 $n$ を固定し $m=n+1$ とすると、$\operatorname{Re}z>-(n+1)$ で $\Gamma(z)=\Gamma(z+n+1)/\bigl(z(z+1)\cdots(z+n)\bigr)$ である。$q(z):=z(z+1)\cdots(z+n-1)$($n=0$ のときは $q\equiv1$)とおくと $q(-n)=(-n)(-n+1)\cdots(-1)=(-1)^nn!\ne0$ であり、$z=-n$ の近傍で $(z+n)\Gamma(z)=\Gamma(z+n+1)/q(z)$ は正則で、$z=-n$ での値は $\Gamma(1)/q(-n)=(-1)^n/n!\ne0$ である。よって $z=-n$ は1位の極で留数は $(-1)^n/n!$ である。他の点では十分大きい $m$ をとれば $\Omega_m$ 上正則なので、極は $\{0,-1,-2,\ldots\}$ に限る。$\square$

性質

右半平面での正則性と導関数

$\Gamma$ は $\{z\mid\operatorname{Re}z>0\}$ 上で正則であり、その導関数は
$$ \Gamma'(z)=\int_0^\infty(\log t)\,t^{z-1}e^{-t}\,dt\qquad(\operatorname{Re}z>0) $$
で与えられる(右辺は絶対収束する)。

正の整数 $n$ に対し $\Gamma_n(z):=\int_{1/n}^{n}t^{z-1}e^{-t}\,dt$ とおく。被積分関数は $(t,z)$ について連続で $z$ について正則なので、$\Gamma_n$ は $\mathbb{C}$ 全体で正則であり、$\Gamma_n'(z)=\int_{1/n}^n(\log t)t^{z-1}e^{-t}\,dt$ である(有界閉区間上のパラメータつき積分の正則性は、積分を Riemann 和の一様極限として Ahl79 Chapter 5 §1.1 の Weierstrass の定理に帰着させるか、Fubiniの定理 と Moreraの定理 を用いて示される標準的な事実であり、本記事では証明を割愛する)。
$0<\alpha<\beta$ とし、帯領域 $S:=\{z\mid\alpha\le\operatorname{Re}z\le\beta\}$ 上で $\Gamma_n\to\Gamma$ が一様収束することを示す。$z\in S$ なら
$$ |\Gamma(z)-\Gamma_n(z)|\le\int_0^{1/n}t^{\sigma-1}e^{-t}\,dt+\int_n^\infty t^{\sigma-1}e^{-t}\,dt \le\int_0^{1/n}t^{\alpha-1}\,dt+C_{\beta-1}\int_n^\infty e^{-t/2}\,dt =\frac{n^{-\alpha}}{\alpha}+2C_{\beta-1}e^{-n/2} $$
である(第1項では $t\le1$, $\sigma\ge\alpha$ より $t^{\sigma-1}\le t^{\alpha-1}$、第2項では $t\ge1$, $\sigma\le\beta$ より $t^{\sigma-1}\le t^{\beta-1}$ と lem-gamma-exponential-bound)。右辺は $z$ によらず $0$ に収束する。正則関数列の局所一様収束極限は正則であり、導関数も局所一様収束する(Weierstrass の定理。Ahl79 Chapter 5 §1.1)ので、$\Gamma$ は $\operatorname{Re}z>0$ で正則で $\Gamma'=\lim\Gamma_n'$ である。最後に $|(\log t)t^{z-1}e^{-t}|=|\log t|\,t^{\sigma-1}e^{-t}$ は、$t\le1$ で $|\log t|\le t^{-\alpha/2}\cdot(2/(e\alpha))$($0< s\le1$, $a>0$ に対する $s^{a}\log(1/s)\le1/(ea)$ を $a=\alpha/2$ で用いる。左辺は $s=e^{-1/a}$ で最大値 $1/(ea)$ をとる)、$t\ge1$ で $\log t\le t$ なので、thm-gamma-convergence と同様に $(0,\infty)$ 上可積分であり、$\Gamma_n'(z)\to\int_0^\infty(\log t)t^{z-1}e^{-t}\,dt$ である。$\square$

対数微分の級数表示

$\Gamma$ の対数微分(対数導関数) $\psi(z):=\Gamma'(z)/\Gamma(z)$ を digamma関数という。$z\notin\{0,-1,-2,\ldots\}$ に対し
$$ \psi(z)=-\gamma-\frac1z+\sum_{n=1}^\infty\frac{z}{n(z+n)}=-\gamma-\frac1z+\sum_{n=1}^\infty\Bigl(\frac1n-\frac1{z+n}\Bigr) $$
が成り立つ。特に $\Gamma'(1)=\psi(1)=-\gamma$ であり、$\gamma=-\int_0^\infty(\log t)e^{-t}\,dt$ である。

rem-gamma-weierstrass-entire により、$F(z):=z\,e^{\gamma z}\prod_{n\ge1}(1+z/n)e^{-z/n}=1/\Gamma(z)$ は $\mathbb{C}$ 上正則で、有限積 $F_N(z):=z\,e^{\gamma z}\prod_{n=1}^N(1+z/n)e^{-z/n}$ の局所一様収束極限である。Weierstrass の定理(Ahl79 Chapter 5 §1.1)により $F_N'\to F'$ も局所一様収束する。$z\notin\{0,-1,\ldots,-N\}$ で
$$ \frac{F_N'(z)}{F_N(z)}=\frac1z+\gamma+\sum_{n=1}^N\Bigl(\frac1{z+n}-\frac1n\Bigr)=\frac1z+\gamma-\sum_{n=1}^N\frac{z}{n(z+n)} $$
であり、$F(z)\ne0$ なる $z$、すなわち $z\notin\{0,-1,-2,\ldots\}$ で $F_N'(z)/F_N(z)\to F'(z)/F(z)=-\Gamma'(z)/\Gamma(z)=-\psi(z)$ である。よって $\psi(z)=-\gamma-\frac1z+\sum_{n\ge1}\frac{z}{n(z+n)}$ を得る。$z=1$ では $\sum_{n\ge1}\frac1{n(n+1)}=\sum_{n\ge1}\bigl(\frac1n-\frac1{n+1}\bigr)=1$ なので $\psi(1)=-\gamma-1+1=-\gamma$ であり、$\Gamma(1)=1$ から $\Gamma'(1)=-\gamma$ である。最後の式は prop-gamma-holomorphic の導関数の積分表示で $z=1$ とおけば従う。$\square$

$\Gamma$ は $(0,\infty)$ 上で対数凸、すなわち $\log\Gamma$ は凸関数である。逆に、$(0,\infty)$ 上の正値関数 $f$ で $f(1)=1$、$f(x+1)=xf(x)$ を満たし $\log f$ が凸なものは $\Gamma$ に限る。これらの証明と、$\Gamma$ についての等式の証明への応用は Bohr–Mollerupの定理 にある。

相反公式・Beta関数・Legendreの倍数公式

次が成り立つ。

  1. (相反公式)整数でない複素数 $z$ に対し $\Gamma(z)\Gamma(1-z)=\dfrac{\pi}{\sin\pi z}$。
  2. (Beta関数)$\operatorname{Re}x>0$, $\operatorname{Re}y>0$ に対し $B(x,y):=\displaystyle\int_0^1t^{x-1}(1-t)^{y-1}\,dt=\frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}$。
  3. (Legendreの倍数公式、ベータ関数)$z\notin\{0,-\tfrac12,-1,-\tfrac32,\ldots\}$ に対し $\Gamma(z)\Gamma\bigl(z+\tfrac12\bigr)=2^{1-2z}\sqrt{\pi}\,\Gamma(2z)$。
古典的公式の証明の所在

thm-gamma-classical-formulas の証明はそれぞれ次の記事にある。(1) は 無限積の定理「Γ関数の相反公式」で、正弦関数の無限積表示と thm-gamma-weierstrass を突き合わせて示される(WW96 §12.14, Ahl79 Chapter 5 §2.4, Art15 §4)。(2) は Beta関数の定理「Γ関数による表示」(WW96 §12.41)で、二重積分の順序を入れ替えて示される。(3) は Beta関数の定理「Legendre の倍数公式」(WW96 §12.15, Art15 §3)で、$B(z,z)$ を 2 通りに計算して $\operatorname{Re}z>0$ の場合を示し、一致の定理で広げる。相反公式は Beta 関数の積分を留数定理で計算しても得られ、Beta 関数との関係は対数凸性と Bohr–Mollerupの定理 からも導ける(Rud76 Theorem 8.20, Art15 §2)。

$\Gamma$ 関数の零点

$\Gamma$ は $\mathbb{C}\setminus\{0,-1,-2,\ldots\}$ 上で零点を持たない。したがって $1/\Gamma$ は $\mathbb{C}$ 全体で正則な関数であり、その零点は $z=0,-1,-2,\ldots$ に限る。

rem-gamma-weierstrass-entire により $\Gamma=1/F$ であり、$F$ は $\mathbb{C}$ 上の整関数でその零点は $z=0,-1,-2,\ldots$ に限る。したがって $\Gamma$ は定義域のどこでも零にならず、$1/\Gamma=F$ の零点も $z=0,-1,-2,\ldots$ に限る。相反公式(thm-gamma-classical-formulas (1))からも同じ結論が得られる:$z$ が整数でないとき $\Gamma(z)\Gamma(1-z)=\pi/\sin\pi z\ne0$ より $\Gamma(z)\ne0$ であり、$z$ が正の整数なら $\Gamma(z)=(z-1)!\ne0$ である。$\square$

Stirling の漸近公式

実数 $x\to+\infty$ のとき
$$ \Gamma(x+1)\sim\sqrt{2\pi x}\Bigl(\frac xe\Bigr)^{x},\qquad\text{特に}\qquad n!\sim\sqrt{2\pi n}\Bigl(\frac ne\Bigr)^{n}\quad(n\to\infty) $$
が成り立つ($\sim$ は比が $1$ に収束することを表す。Landau記号)。さらに $|\arg z|\le\pi-\delta$($\delta>0$)を満たす $z$ について $|z|\to\infty$ のとき、$\log\Gamma(z)=(z-\tfrac12)\log z-z+\tfrac12\log(2\pi)+O(1/|z|)$ という漸近展開が成り立つ。

Stirling の公式の証明の所在

thm-gamma-stirling の証明は本記事では扱わない。$n!$ についての形は Stirlingの公式 で証明している。実数 $x\to\infty$ での主部 $\Gamma(x+1)\sim\sqrt{2\pi x}\,(x/e)^x$ は Rud76 8.22 の (103)、複素数 $z$ での $\log\Gamma(z)$ の漸近展開は WW96 §12.33 と AAR99 Chapter 1 を引用する。$n!$ についての形の誤差を上下から押さえる初等的な評価は Rob55 にある。

計算と応用

半整数での値

置換 $t=u^2$($dt=2u\,du$)により
$$ \Gamma\Bigl(\frac12\Bigr)=\int_0^\infty t^{-1/2}e^{-t}\,dt=\int_0^\infty\frac1u\,e^{-u^2}\,2u\,du=2\int_0^\infty e^{-u^2}\,du=\int_{-\infty}^\infty e^{-u^2}\,du=\sqrt\pi $$
である(最後の等号は Gauss積分の値であり、Rud76 8.21 の (101) に従う)。関数等式 prop-gamma-functional-equation により $\Gamma(3/2)=\frac12\Gamma(1/2)=\frac{\sqrt\pi}2$、一般に非負整数 $n$ に対し
$$ \Gamma\Bigl(n+\frac12\Bigr)=\frac{(2n)!}{4^n\,n!}\sqrt\pi $$
である($n$ に関する帰納法:$\Gamma(n+\tfrac32)=(n+\tfrac12)\Gamma(n+\tfrac12)=\frac{2n+1}{2}\cdot\frac{(2n)!}{4^nn!}\sqrt\pi=\frac{(2n+2)!}{4^{n+1}(n+1)!}\sqrt\pi$)。負の半整数では thm-gamma-continuation (1) により $\Gamma(-\tfrac12)=\Gamma(\tfrac12)/(-\tfrac12)=-2\sqrt\pi$ である。

尺度変換と Laplace 型積分

$a>0$, $\operatorname{Re}z>0$ に対し、置換 $t=au$ により
$$ \int_0^\infty u^{z-1}e^{-au}\,du=\frac{\Gamma(z)}{a^{z}} $$
である。$z=n+1$($n$ は非負整数)なら右辺は $n!/a^{n+1}$ であり、これは $u^n$ のLaplace変換である。また $p>0$ に対し置換 $t=u^p$ により $\int_0^\infty e^{-u^p}\,du=\Gamma(1+1/p)$ であり、$p=2$ で Gauss 積分 $\int_0^\infty e^{-u^2}\,du=\Gamma(3/2)=\sqrt\pi/2$ を再現する。

原点近くでの展開

thm-gamma-continuation により $\Gamma(z)=\Gamma(z+1)/z$ であり、$\Gamma$ は $z=1$ で正則で $\Gamma(1)=1$, $\Gamma'(1)=-\gamma$(prop-gamma-digamma)なので、$\Gamma(z+1)=1-\gamma z+O(z^2)$($z\to0$)である。よって
$$ \Gamma(z)=\frac1z-\gamma+O(z)\quad(z\to0) $$
が成り立ち、$z=0$ が留数 $1$ の1位の極であることと、その定数項が $-\gamma$ であることがわかる。

反例:階乗を延長する関数の非一意性

$f(x):=\Gamma(x)\bigl(1+\tfrac12\sin2\pi x\bigr)$($x>0$)は、正値であるという性質、$f(1)=1$ という性質、$f(x+1)=xf(x)$ という性質を満たす($\sin2\pi(x+1)=\sin2\pi x$ による)が、対数凸であるという性質を満たさない。実際、もし対数凸なら Bohr–Mollerupの定理 により $f=\Gamma$ となるが、$f(1/4)=\Gamma(1/4)\cdot\tfrac32\ne\Gamma(1/4)$ である。したがって「正値で $f(1)=1$, $f(x+1)=xf(x)$ を満たす $\Rightarrow$ $f=\Gamma$」という含意は成り立たず、Bohr–Mollerup の定理の対数凸性の仮定は省けない。

反例:積分表示の収束域

$z=-1/2$ に対し $\Gamma(-1/2)=-2\sqrt\pi$ は ex-gamma-half-integer のとおり解析接続によって定まるが、$\int_0^\infty t^{-3/2}e^{-t}\,dt$ は thm-gamma-convergence (2) により $+\infty$ に発散する。したがって $\operatorname{Re}z\le0$ では「$\Gamma(z)$ は積分 $\int_0^\infty t^{z-1}e^{-t}\,dt$ に等しい」という主張は成り立たず、定義の積分表示を $\operatorname{Re}z>0$ の外で使ってはならない。

補足

歴史と記法。 階乗を整数以外へ延長する問題は Euler が1729年に Goldbach への手紙で扱い、fml-gamma-euler-product の形の乗積と第二種 Euler 積分を与えた。記号 $\Gamma$ と名称は Legendre による。Gauss は $\Pi(z):=\Gamma(z+1)$($\Pi(n)=n!$)を用い、Weierstrass は thm-gamma-weierstrass の乗積を出発点に採った(以上 WW96 §12.1, AAR99 Chapter 1)。$\Gamma(z+1)=z\Gamma(z)$ のため $\Gamma(n)=(n-1)!$ と添字がずれることに注意する。
定義の流儀。 本記事は第二種 Euler 積分を定義とし、Gauss の極限公式・Weierstrass の乗積表示を定理として導いた。Weierstrass の乗積表示を定義とする流儀(WW96)では、$\Gamma$ が最初から $\mathbb{C}\setminus\{0,-1,-2,\ldots\}$ で定義され、Euler 積分は定理になる(WW96 §12.2)。$\Gamma$ を $(0,\infty)$ 上の関数として、Bohr–Mollerupの定理 の三条件で特徴づける流儀(Art15, Rud76 Theorem 8.19)もある。いずれの流儀でも得られる関数は同じである。
関連する関数。 Beta関数 $B(x,y)$ は thm-gamma-classical-formulas (2) により $\Gamma$ で表される。digamma関数 $\psi=\Gamma'/\Gamma$ とその導関数(polygamma関数)は prop-gamma-digamma の級数から得られる。Riemannゼータ関数の関数等式 $\pi^{-s/2}\Gamma(s/2)\zeta(s)=\pi^{-(1-s)/2}\Gamma((1-s)/2)\zeta(1-s)$ には $\Gamma$ が現れ、$\zeta$ の自明な零点は $\Gamma(s/2)$ の極から生じる(WW96 §13.15)。確率論ではガンマ分布の密度 $x^{k-1}e^{-x/\theta}/(\Gamma(k)\theta^{k})$ の正規化定数として現れる。$n$ 次元単位球の体積が $\pi^{n/2}/\Gamma(n/2+1)$ であることも $\Gamma$ の代表的な応用である(AAR99 Chapter 1)。

関連項目

参考文献

[1]
E. T. Whittaker and G. N. Watson, A Course of Modern Analysis, Cambridge Mathematical Library, Cambridge University Press, 1996, Chapter XII(§12.1 Weierstrass の定義と Euler の乗積、§12.14 相反公式、§12.15 倍数公式、§12.16 対数微分、§12.2 Euler 積分、§12.33 Stirling の級数、§12.41 Beta 関数)、§13.15(ζ の関数等式)
[2]
Emil Artin (translated by Michael Butler), The Gamma Function, Dover Books on Mathematics, Dover Publications, 2015, §2(Euler 積分と Gauss の公式、Beta 関数)、§3(倍数公式)、§4(sin との関係:相反公式)
[3]
George E. Andrews, Richard Askey and Ranjan Roy, Special Functions, Encyclopedia of Mathematics and its Applications 71, Cambridge University Press, 1999, Chapter 1(The Gamma and Beta Functions:歴史、Bohr–Mollerup の定理、Stirling の公式、球の体積)
[4]
[6]
Walter Rudin, Principles of Mathematical Analysis, 3rd ed., McGraw-Hill, 1976, Chapter 8(8.19 Bohr–Mollerup の定理、8.20 Beta 関数、8.21 (101) Gauss 積分、8.22 Stirling の公式(実数 x→∞ の主部のみ))
[7]
Lars V. Ahlfors, Complex Analysis, 3rd ed., McGraw-Hill, 1979, Chapter 5 §1.1(Weierstrass の定理:正則関数列の一様収束)、§2.2(無限積)、§2.4(Γ関数)

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する