Γ関数(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$ も成り立つ。
前提知識: 階乗, 広義積分, 部分積分, 正則関数, 凸関数
以下、複素数 $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}$ である。複素数値関数の広義積分は、実部・虚部それぞれの広義積分の和として定める。
$\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$
複素数 $z$ について次が成り立つ。
$|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$
$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$
$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$
$\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$
$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 に譲る。
$\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$ について次が成り立つ。
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$
区間 $I$ 上の正値関数 $f\colon I\to(0,\infty)$ が対数凸(logarithmically convex)であるとは、$\log f$ が $I$ 上の凸関数であること、すなわち任意の $x,y\in I$ と $\lambda\in[0,1]$ に対し
$$
f(\lambda x+(1-\lambda)y)\le f(x)^{\lambda}f(y)^{1-\lambda}
$$
が成り立つことをいう。
$a,b\ge0$, $\lambda\in(0,1)$ に対し $a^{\lambda}b^{1-\lambda}\le\lambda a+(1-\lambda)b$ が成り立つ(相加相乗平均の不等式の重みつき一般化)。
$a=0$ または $b=0$ なら左辺は $0$ で自明である。$a,b>0$ とし $s:=a/b$ とおくと、$b$ で割って $s^{\lambda}\le\lambda s+(1-\lambda)$ を示せばよい。$g(s):=\lambda s+(1-\lambda)-s^{\lambda}$ は $g'(s)=\lambda(1-s^{\lambda-1})$ であり、$\lambda-1<0$ なので $0< s<1$ で $g'<0$、$s>1$ で $g'>0$ である。よって $g$ は $s=1$ で最小値 $g(1)=0$ をとり、$g\ge0$ である。$\square$
$F,G\colon(0,\infty)\to[0,\infty)$ を連続関数とし、$A:=\int_0^\infty F$, $B:=\int_0^\infty G$ がともに収束するとする。$\lambda\in(0,1)$ に対し $\int_0^\infty F^{\lambda}G^{1-\lambda}$ は収束し、
$$
\int_0^\infty F^{\lambda}G^{1-\lambda}\le A^{\lambda}B^{1-\lambda}
$$
が成り立つ(Hölderの不等式の特別な場合)。
$H:=F^{\lambda}G^{1-\lambda}$ とおく。lem-gamma-weighted-am-gm($a=F(t)$, $b=G(t)$)により $0\le H\le\lambda F+(1-\lambda)G$ であり、右辺の広義積分は収束するので、比較判定法により $\int_0^\infty H$ は収束する。$A=0$ なら、$F$ は連続かつ非負なので $F\equiv0$ となり $H\equiv0$ で不等式は自明である。$B=0$ も同様である。$A,B>0$ とする。lem-gamma-weighted-am-gm を $a=F(t)/A$, $b=G(t)/B$ に適用すると
$$
\frac{H(t)}{A^{\lambda}B^{1-\lambda}}\le\lambda\frac{F(t)}{A}+(1-\lambda)\frac{G(t)}{B}
$$
であり、$0<\varepsilon< R$ で積分すると $\int_\varepsilon^RH\le A^{\lambda}B^{1-\lambda}\bigl(\lambda\int_\varepsilon^RF/A+(1-\lambda)\int_\varepsilon^RG/B\bigr)\le A^{\lambda}B^{1-\lambda}(\lambda+(1-\lambda))=A^{\lambda}B^{1-\lambda}$ である。$\varepsilon\to0$, $R\to\infty$ として結論を得る。$\square$
$\Gamma$ は $(0,\infty)$ 上正値で対数凸である。すなわち $x,y>0$, $\lambda\in[0,1]$ に対し $\Gamma(\lambda x+(1-\lambda)y)\le\Gamma(x)^{\lambda}\Gamma(y)^{1-\lambda}$ が成り立つ。
$x>0$ なら被積分関数 $t^{x-1}e^{-t}$ は正なので $\Gamma(x)>0$ である。$\lambda\in(0,1)$ とし($\lambda=0,1$ は自明)、$F(t):=t^{x-1}e^{-t}$, $G(t):=t^{y-1}e^{-t}$ とおくと、thm-gamma-convergence により $\int_0^\infty F=\Gamma(x)$, $\int_0^\infty G=\Gamma(y)$ は収束する。指数法則により $F^{\lambda}G^{1-\lambda}=t^{(\lambda x+(1-\lambda)y)-1}e^{-t}$ なので、lem-gamma-holder により $\Gamma(\lambda x+(1-\lambda)y)=\int_0^\infty F^{\lambda}G^{1-\lambda}\le\Gamma(x)^{\lambda}\Gamma(y)^{1-\lambda}$ である。$\square$
区間 $I$ 上の凸関数 $\varphi$ と $I$ の点 $a< b< c$ に対し
$$
\frac{\varphi(b)-\varphi(a)}{b-a}\le\frac{\varphi(c)-\varphi(a)}{c-a}\le\frac{\varphi(c)-\varphi(b)}{c-b}
$$
が成り立つ。
$\mu:=(c-b)/(c-a)\in(0,1)$ とおくと $\mu a+(1-\mu)c=b$ であり、凸性から $\varphi(b)\le\mu\varphi(a)+(1-\mu)\varphi(c)$ である。両辺から $\varphi(a)$ を引くと $\varphi(b)-\varphi(a)\le(1-\mu)(\varphi(c)-\varphi(a))=\frac{b-a}{c-a}(\varphi(c)-\varphi(a))$ となり、第1の不等式を得る。両辺を $\varphi(c)$ から引くと $\varphi(c)-\varphi(b)\ge\mu(\varphi(c)-\varphi(a))=\frac{c-b}{c-a}(\varphi(c)-\varphi(a))$ となり、第2の不等式を得る。$\square$
$f\colon(0,\infty)\to(0,\infty)$ が (1) $f(1)=1$、(2) すべての $x>0$ で $f(x+1)=xf(x)$、(3) $f$ は対数凸、を満たすならば、$x>0$ で $f(x)=\Gamma(x)$ である。すなわち、この三条件は $(0,\infty)$ 上の $\Gamma$ を特徴づける。
$\varphi:=\log f$ は凸で、$\varphi(1)=0$, $\varphi(x+1)=\varphi(x)+\log x$ を満たす。帰納法により正の整数 $n$ に対し $\varphi(n+1)=\log(n!)$ である。
段階1:$0< x\le1$、$n$ を正の整数とする。三点 $n< n+1< n+1+x$ に lem-gamma-convex-slope を適用すると $\log n=\varphi(n+1)-\varphi(n)\le\dfrac{\varphi(n+1+x)-\varphi(n+1)}{x}$ であり、$x<1$ のとき三点 $n+1< n+1+x< n+2$ に適用すると $\dfrac{\varphi(n+1+x)-\varphi(n+1)}{x}\le\varphi(n+2)-\varphi(n+1)=\log(n+1)$ である($x=1$ では等号で成立)。よって
$$
x\log n\le\varphi(n+1+x)-\log(n!)\le x\log(n+1),\qquad\text{すなわち}\qquad n!\,n^{x}\le f(n+1+x)\le n!\,(n+1)^{x}
$$
である。
段階2:条件 (2) を繰り返すと $f(n+1+x)=f(x)\,x(x+1)\cdots(x+n)$ である。これを段階1に代入して $x(x+1)\cdots(x+n)>0$ で割ると
$$
a_n(x)\le f(x)\le a_n(x)\Bigl(1+\frac1n\Bigr)^{x},\qquad a_n(x):=\frac{n!\,n^{x}}{x(x+1)\cdots(x+n)}
$$
を得る。
段階3:$(1+1/n)^{x}\to1$ なので、はさみうちにより $a_n(x)\to f(x)$($0< x\le1$)である。$a_n(x)$ は $f$ によらないので、条件 (1)(2)(3) を満たす関数はすべて $(0,1]$ 上で一致する。$\Gamma$ は prop-gamma-functional-equation と thm-gamma-log-convexity により (1)(2)(3) を満たすので、$0< x\le1$ で $f(x)=\Gamma(x)$ である(この $a_n(x)$ の極限は thm-gamma-gauss-limit とも整合する)。
段階4:一般の $x>0$ を $x=y+m$($y\in(0,1]$, $m$ は非負整数)と書き、$m$ に関する帰納法で $f(y+m)=\Gamma(y+m)$ を示す。$m=0$ は段階3である。$f(y+m)=\Gamma(y+m)$ なら $f(y+m+1)=(y+m)f(y+m)=(y+m)\Gamma(y+m)=\Gamma(y+m+1)$ である。$\square$
次が成り立つ。
thm-gamma-classical-formulas の証明は本記事では割愛し、標準的な教科書を引用する。相反公式は、$\sin\pi z$ の無限積展開と thm-gamma-weierstrass を突き合わせるか、Beta 関数の積分を留数定理で計算して得られる(WW96 §12.14, Ahl79 Chapter 5 §2.4, Art15 §4)。Beta 関数との関係は Rud76 Chapter 8 および WW96 §12.41、Legendre の倍数公式は WW96 §12.15 と Art15 §3 による。Beta 関数との関係は、対数凸性と Bohr–Mollerup の定理からも導ける(Art15 §2)。
$\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$
実数 $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|)$ という漸近展開が成り立つ。
thm-gamma-stirling(Stirlingの公式)は教科書レベルを超える精密な評価を含むので、本記事では証明を割愛し、Rud76 Chapter 8、WW96 §12.33、AAR99 Chapter 1 を引用する。実数の場合の初等的な証明は 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 Chapter 8 に従う)。関数等式 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$ である。
$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$ による)が、対数凸であるという性質を満たさない。実際、もし対数凸なら thm-gamma-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 Chapter 8)もある。いずれの流儀でも得られる関数は同じである。
関連する関数。 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)。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する