階乗関数の拡張
ガンマ 実軸の一部に沿ったガンマ関数
一般的な定義
Γ
(
ず
)
=
∫
0
∞
t
ず
−
1
e
−
t
d
t
{\displaystyle \Gamma (z)=\int _{0}^{\infty }t^{z-1}e^{-t}\,dt}
応用分野 微積分、数学的解析、統計学、物理学
数学 において 、 ガンマ関数 (Γ、 ギリシャ 文字の大文字 ガンマで表される)は、 階乗関数の 複素数 への 最も一般的な拡張です 。 ダニエル・ベルヌーイ によって導出されたガンマ関数は、 非正の整数を除く すべての複素数に対して定義され、すべての 正の整数 に対して定義されます。ガンマ関数は 、実部が正の複素数の
収束する 仮積分によって定義できます。
Γ
(
ず
)
{\displaystyle \Gamma (z)}
ず
{\displaystyle z}
ず
=
ん
{\displaystyle z=n}
Γ
(
ん
)
=
(
ん
−
1
)
!
。
{\displaystyle \ガンマ (n)=(n-1)!\,.}
Γ
(
ず
)
=
∫
0
∞
t
ず
−
1
e
−
t
d
t
、
ℜ
(
ず
)
>
0
。
{\displaystyle \Gamma (z)=\int _{0}^{\infty }t^{z-1}e^{-t}{\text{ d}}t,\ \qquad \Re (z)>0\,.}
ガンマ関数は、この積分関数の 解析接続 として複素平面上で定義されます。これは、ゼロと負 の整数を除いて正則な 有理 型関数 であり、ゼロと負の整数では単純な 極を 持ちます。
ガンマ関数にはゼロがないので、 逆ガンマ関数 は 1 / Γ( z ) は は 完全な関数 です 。実際、ガンマ関数は 負の 指数関数の メリン変換 に対応します。
Γ
(
ず
)
=
ま
{
e
−
x
}
(
ず
)
。
{\displaystyle \Gamma (z)={\mathcal {M}}\{e^{-x}\}(z)\,.}
階乗関数の他の拡張も存在しますが、ガンマ関数が最も一般的で便利です。これは、 確率 、 統計 、 解析的数論 、 組合せ論 の分野におけるさまざまな確率分布関数やその他の式の因子として使用されます。
モチベーション
Γ
(
x
+
1
)
{\displaystyle \Gamma (x+1)}
階乗関数を非整数値に補間します。
ガンマ関数は、階乗列の点を結ぶ 滑らかな曲線 を見つける 補間 問題の解として考えることができます。 のすべての正の整数値に対して 。階乗の簡単な式 x ! = 1 × 2 × ⋯ × x は、 x が 正の整数の場合にのみ有効であり 、この特性を持つ 基本関数は ありませんが、ガンマ関数は優れた解です 。 [1]
ええ
=
ふ
(
x
)
{\displaystyle y=f(x)}
(
x
、
ええ
)
=
(
ん
、
ん
!
)
{\displaystyle (x,y)=(n,n!)}
ん
{\displaystyle n}
ふ
(
x
)
=
Γ
(
x
+
1
)
{\displaystyle f(x)=\Gamma (x+1)}
ガンマ関数は滑らかであるだけでなく 解析的 (非正の整数を除く)であり、いくつかの明示的な方法で定義できます。ただし、これは階乗を拡張する唯一の解析関数ではありません。 整数 など、正の整数に対してゼロになる任意の解析関数を追加できます 。 [1] このような関数は 擬似ガンマ関数 と呼ばれ、最も有名なのは アダマール 関数です。 [2]
け
罪
(
メートル
π
x
)
{\displaystyle k\sin(m\pi x)}
メートル
{\displaystyle m}
ガンマ関数 Γ( z ) (青) が、 Γ( z ) + sin(π z ) (緑)とともにプロットされています 。正の整数での交差に注意してください。どちらも、複素平面上の有理型関数への階乗の有効な拡張です。
より厳しい要件は、 シフトされた階乗を補間する 関数方程式 である : [3] [4]
ふ
(
ん
)
=
(
ん
−
1
)
!
{\displaystyle f(n)=(n{-}1)!}
ふ
(
x
+
1
)
=
x
ふ
(
x
)
いかなる
x
>
0
、
ふ
(
1
)
=
1.
{\displaystyle f(x+1)=xf(x)\ {\text{ 任意のx>0に対して\qquad f(1)=1.}
しかし、および を含む任意の周期関数 ( など) による乗算が可能であるため、これはまだ一意の解を与えません 。
グ
(
x
)
{\displaystyle g(x)}
グ
(
x
)
=
グ
(
x
+
1
)
{\displaystyle g(x)=g(x+1)}
グ
(
0
)
=
1
{\displaystyle g(0)=1}
グ
(
x
)
=
e
け
罪
(
メートル
π
x
)
{\displaystyle g(x)=e^{k\sin(m\pi x)}}
この曖昧さを解決する一つの方法は ボーア・モレルプ定理 であり、これは 階乗の唯一の補間関数であり、正の実数上で定義され、 対数的に凸であることを示しています [5] 、 つまり 凸 です 。 [6]
ふ
(
x
)
=
Γ
(
x
)
{\displaystyle f(x)=\Gamma (x)}
ええ
=
ログ
ふ
(
x
)
{\displaystyle y=\log f(x)}
意味
主な定義
この表記は ルジャンドル によるものです 。 [1] 複素数 z の実部が厳密に正 ( ) の場合、 積分は 絶対収束し 、 第二種オイラー積分 として知られています。 (第一種オイラー積分は ベータ関数 です。 [1] ) 部分積分を 使用すると、次のことがわかります。
Γ
(
ず
)
{\displaystyle \Gamma (z)}
ℜ
(
ず
)
>
0
{\displaystyle \Re (z)>0}
Γ
(
ず
)
=
∫
0
∞
t
ず
−
1
e
−
t
d
t
{\displaystyle \Gamma (z)=\int _{0}^{\infty }t^{z-1}e^{-t}\,dt}
複素平面上のガンマ関数の絶対値(垂直)と引数(色)
Γ
(
ず
+
1
)
=
∫
0
∞
t
ず
e
−
t
d
t
=
[
−
t
ず
e
−
t
]
0
∞
+
∫
0
∞
ず
t
ず
−
1
e
−
t
d
t
=
リム
t
→
∞
(
−
t
ず
e
−
t
)
−
(
−
0
ず
e
−
0
)
+
ず
∫
0
∞
t
ず
−
1
e
−
t
d
t
。
{\displaystyle {\begin{aligned}\Gamma (z+1)&=\int _{0}^{\infty }t^{z}e^{-t}\,dt\\&={\Bigl [}-t^{z}e^{-t}{\Bigr ]}_{0}^{\infty }+\int _{0}^{\infty }zt^{z-1}e^{-t}\,dt\\&=\lim _{t\to \infty }\left(-t^{z}e^{-t}\right)-\left(-0^{z}e^{-0}\right)+z\int _{0}^{\infty }t^{z-1}e^{-t}\,dt.\end{aligned}}}
それを認識し て
−
t
ず
e
−
t
→
0
{\displaystyle -t^{z}e^{-t}\to 0}
t
→
∞
、
{\displaystyle t\to \infty ,}
Γ
(
ず
+
1
)
=
ず
∫
0
∞
t
ず
−
1
e
−
t
d
t
=
ず
Γ
(
ず
)
。
{\displaystyle {\begin{aligned}\Gamma (z+1)&=z\int _{0}^{\infty }t^{z-1}e^{-t}\,dt\\&=z\Gamma (z).\end{aligned}}}
次の ように計算できます。
Γ
(
1
)
{\displaystyle \Gamma (1)}
Γ
(
1
)
=
∫
0
∞
t
1
−
1
e
−
t
d
t
=
∫
0
∞
e
−
t
d
t
=
1.
{\displaystyle {\begin{aligned}\Gamma (1)&=\int _{0}^{\infty }t^{1-1}e^{-t}\,dt\\&=\int _{0}^{\infty }e^{-t}\,dt\\&=1.\end{aligned}}}
したがって、任意 の正の整数 n に対して、帰納法によって が 成り立つことが示せます 。具体的には、基本ケースは であり 、帰納法のステップは です。
Γ
(
ん
)
=
(
ん
−
1
)
!
{\displaystyle \Gamma (n)=(n-1)!}
Γ
(
1
)
=
1
=
0
!
{\displaystyle \ガンマ (1)=1=0!}
Γ
(
ん
+
1
)
=
ん
Γ
(
ん
)
=
ん
(
ん
−
1
)
!
=
ん
!
。
{\displaystyle \Gamma (n+1)=n\Gamma (n)=n(n-1)!=n!.}
この恒等式 を使うと(あるいは同じ結果になる 解析接続 を使うと) 、 0以下の整数を除く すべての複素数 zに対して定義される 有理型関数 への積分定式化を一意に拡張することができる。 [1] この拡張版は一般にガンマ関数と呼ばれている。 [1]
Γ
(
ず
)
=
Γ
(
ず
+
1
)
ず
{\textstyle \Gamma (z)={\frac {\Gamma (z+1)}{z}}}
Γ
(
ず
)
{\displaystyle \Gamma (z)}
代替定義
同等の定義は多数あります。
オイラーの無限積の定義
固定された整数に対して 、整数が 増加すると、 [7]
メートル
{\displaystyle m}
ん
{\displaystyle n}
リム
ん
→
∞
ん
!
(
ん
+
1
)
メートル
(
ん
+
メートル
)
!
=
1
。
{\displaystyle \lim _{n\to \infty}{\frac {n!\,\left(n+1\right)^{m}}{(n+m)!}}=1\,.}
が整数でない 場合、この式は意味がありません。このセクションでは非整数の階乗がまだ定義されていないためです。ただし、 非整数のガンマ関数を定義するために、
を任意の複素数 に置き換えてもこの式が成立すると仮定します。
メートル
{\displaystyle m}
メートル
{\displaystyle m}
ず
{\displaystyle z}
リム
ん
→
∞
ん
!
(
ん
+
1
)
ず
(
ん
+
ず
)
!
=
1
。
{\displaystyle \lim _{n\to \infty}{\frac {n!\,\left(n+1\right)^{z}}{(n+z)!}}=1\,.}
両辺に を掛けると が 得られ
ます。 この 無限積 はオイラー [8] によるもので、非正の整数を除く すべての複素数に対して収束します 。非正の整数はゼロ除算のため収束しません。したがって、上記の仮定により の一意の定義が得られます 。
(
ず
−
1
)
!
{\displaystyle (z-1)!}
Γ
(
ず
)
=
(
ず
−
1
)
!
=
1
ず
リム
ん
→
∞
ん
!
ず
!
(
ん
+
ず
)
!
(
ん
+
1
)
ず
=
1
ず
リム
ん
→
∞
(
1
⋅
2
⋯
ん
)
1
(
1
+
ず
)
⋯
(
ん
+
ず
)
(
2
1
⋅
3
2
⋯
ん
+
1
ん
)
ず
=
1
ず
∏
ん
=
1
∞
[
1
1
+
ず
ん
(
1
+
1
ん
)
ず
]
。
{\displaystyle {\begin{aligned}\Gamma (z)&=(z-1)!\\[8pt]&={\frac {1}{z}}\lim _{n\to \infty }n!{\frac {z!}{(n+z)!}}(n+1)^{z}\\[8pt]&={\frac {1}{z}}\lim _{n\to \infty }(1\cdot 2\cdots n){\frac {1}{(1+z)\cdots (n+z)}}\left({\frac {2}{1}}\cdot {\frac {3}{2}}\cdots {\frac {n+1}{n}}\right)^{z}\\[8pt]&={\frac {1}{z}}\prod _{n=1}^{\infty }\left[{\frac {1}{1+{\frac {z}{n}}}}\left(1+{\frac {1}{n}}\right)^{z}\right].\end{aligned}}}
ず
{\displaystyle z}
ず
!
{\displaystyle z!}
直感的には、この式は、 が、 ある大きな整数 について 計算し 、 を乗じて を近似し 、関係を 逆に回使用して の 近似値を得た結果に近似していること 、さらに、 が無限大に増加するにつれてこの近似値が正確になることを示唆しています 。
Γ
(
ず
)
{\displaystyle \Gamma (z)}
Γ
(
ん
+
1
)
=
ん
!
{\displaystyle \Gamma (n+1)=n!}
ん
{\displaystyle n}
(
ん
+
1
)
ず
{\displaystyle (n+1)^{z}}
Γ
(
ん
+
ず
+
1
)
{\displaystyle \ガンマ (n+z+1)}
Γ
(
x
+
1
)
=
x
Γ
(
x
)
{\displaystyle \Gamma (x+1)=x\Gamma (x)}
ん
+
1
{\displaystyle n+1}
Γ
(
ず
)
{\displaystyle \Gamma (z)}
ん
{\displaystyle n}
逆数の無限積は
完全な関数
であり 、すべての複素数 z に対して収束します。
1
Γ
(
ず
)
=
ず
∏
ん
=
1
∞
[
(
1
+
ず
ん
)
/
(
1
+
1
ん
)
ず
]
{\displaystyle {\frac {1}{\Gamma (z)}}=z\prod _{n=1}^{\infty }\left[\left(1+{\frac {z}{n}}\right)/{\left(1+{\frac {1}{n}}\right)^{z}}\right]}
ワイエルシュトラスの定義
ワイエルシュトラス によるガンマ関数の定義は、 非正整数を除く
すべての複素数に対しても有効です 。
ここで は オイラー・マスケロニ定数 です 。 [1] これはの アダマール積 を書き直したものです 。この定義は、円周率を含む重要な恒等式に登場します。 [ 要出典 ]
ず
{\displaystyle z}
Γ
(
ず
)
=
e
−
γ
ず
ず
∏
ん
=
1
∞
(
1
+
ず
ん
)
−
1
e
ず
/
ん
、
{\displaystyle \Gamma (z)={\frac {e^{-\gamma z}}{z}}\prod _{n=1}^{\infty }\left(1+{\frac {z}{n}}\right)^{-1}e^{z/n},}
γ
≈
0.577216
{\displaystyle \gamma \approx 0.577216}
1
/
Γ
(
z
)
{\displaystyle 1/\Gamma (z)}
プロパティ
一般的な
上で述べた基本的な性質の他に、
ガンマ関数の他の重要な関数方程式は、 オイラーの反射公式
と
ルジャンドル
の複製公式である。
Γ
(
z
+
1
)
=
z
Γ
(
z
)
{\displaystyle \Gamma (z+1)=z\ \Gamma (z)}
Γ
(
1
−
z
)
Γ
(
z
)
=
π
sin
π
z
,
z
∉
Z
{\displaystyle \Gamma (1-z)\Gamma (z)={\frac {\pi }{\sin \pi z}},\qquad z\not \in \mathbb {Z} }
Γ
(
z
−
n
)
=
(
−
1
)
n
−
1
Γ
(
−
z
)
Γ
(
1
+
z
)
Γ
(
n
+
1
−
z
)
,
n
∈
Z
{\displaystyle \Gamma (z-n)=(-1)^{n-1}\;{\frac {\Gamma (-z)\Gamma (1+z)}{\Gamma (n+1-z)}},\qquad n\in \mathbb {Z} }
Γ
(
z
)
Γ
(
z
+
1
2
)
=
2
1
−
2
z
π
Γ
(
2
z
)
.
{\displaystyle \Gamma (z)\Gamma \left(z+{\tfrac {1}{2}}\right)=2^{1-2z}\;{\sqrt {\pi }}\;\Gamma (2z).}
複製公式は 乗法定理の特殊なケースである( [9] 式5.5.6を
参照 )。
∏
k
=
0
m
−
1
Γ
(
z
+
k
m
)
=
(
2
π
)
m
−
1
2
m
1
2
−
m
z
Γ
(
m
z
)
.
{\displaystyle \prod _{k=0}^{m-1}\Gamma \left(z+{\frac {k}{m}}\right)=(2\pi )^{\frac {m-1}{2}}\;m^{{\frac {1}{2}}-mz}\;\Gamma (mz).}
限界の定義からわかる、単純だが便利な特性は次のとおりです。
Γ
(
z
)
¯
=
Γ
(
z
¯
)
⇒
Γ
(
z
)
Γ
(
z
¯
)
∈
R
.
{\displaystyle {\overline {\Gamma (z)}}=\Gamma ({\overline {z}})\;\Rightarrow \;\Gamma (z)\Gamma ({\overline {z}})\in \mathbb {R} .}
特に、 z = a + bi の場合、この積は
|
Γ
(
a
+
b
i
)
|
2
=
|
Γ
(
a
)
|
2
∏
k
=
0
∞
1
1
+
b
2
(
a
+
k
)
2
{\displaystyle |\Gamma (a+bi)|^{2}=|\Gamma (a)|^{2}\prod _{k=0}^{\infty }{\frac {1}{1+{\frac {b^{2}}{(a+k)^{2}}}}}}
実部が整数または半整数の場合、これは 閉じた形式 で有限に表現できます。
|
Γ
(
b
i
)
|
2
=
π
b
sinh
π
b
|
Γ
(
1
2
+
b
i
)
|
2
=
π
cosh
π
b
|
Γ
(
1
+
b
i
)
|
2
=
π
b
sinh
π
b
|
Γ
(
1
+
n
+
b
i
)
|
2
=
π
b
sinh
π
b
∏
k
=
1
n
(
k
2
+
b
2
)
,
n
∈
N
|
Γ
(
−
n
+
b
i
)
|
2
=
π
b
sinh
π
b
∏
k
=
1
n
(
k
2
+
b
2
)
−
1
,
n
∈
N
|
Γ
(
1
2
±
n
+
b
i
)
|
2
=
π
cosh
π
b
∏
k
=
1
n
(
(
k
−
1
2
)
2
+
b
2
)
±
1
,
n
∈
N
{\displaystyle {\begin{aligned}|\Gamma (bi)|^{2}&={\frac {\pi }{b\sinh \pi b}}\\[1ex]\left|\Gamma \left({\tfrac {1}{2}}+bi\right)\right|^{2}&={\frac {\pi }{\cosh \pi b}}\\[1ex]\left|\Gamma \left(1+bi\right)\right|^{2}&={\frac {\pi b}{\sinh \pi b}}\\[1ex]\left|\Gamma \left(1+n+bi\right)\right|^{2}&={\frac {\pi b}{\sinh \pi b}}\prod _{k=1}^{n}\left(k^{2}+b^{2}\right),\quad n\in \mathbb {N} \\[1ex]\left|\Gamma \left(-n+bi\right)\right|^{2}&={\frac {\pi }{b\sinh \pi b}}\prod _{k=1}^{n}\left(k^{2}+b^{2}\right)^{-1},\quad n\in \mathbb {N} \\[1ex]\left|\Gamma \left({\tfrac {1}{2}}\pm n+bi\right)\right|^{2}&={\frac {\pi }{\cosh \pi b}}\prod _{k=1}^{n}\left(\left(k-{\tfrac {1}{2}}\right)^{2}+b^{2}\right)^{\pm 1},\quad n\in \mathbb {N} \\[-1ex]&\end{aligned}}}
おそらく、非整数引数におけるガンマ関数の最もよく知られた値は です。
これは、 反射または複製の公式で を設定するか、 以下に示す ベータ関数 との関係を で使用するか、ガンマ関数の積分定義で を 代入して ガウス 積分 を得ることによって見つけることができます。一般に、 の非負の整数値に対して、次の式 が成り立ちます。
ここで、 二重階乗 です。計算値については、 「ガンマ関数の特定の値」 を参照してください。
Γ
(
1
2
)
=
π
,
{\displaystyle \Gamma \left({\tfrac {1}{2}}\right)={\sqrt {\pi }},}
z
=
1
2
{\textstyle z={\frac {1}{2}}}
z
1
=
z
2
=
1
2
{\textstyle z_{1}=z_{2}={\frac {1}{2}}}
u
=
z
{\displaystyle u={\sqrt {z}}}
n
{\displaystyle n}
Γ
(
1
2
+
n
)
=
(
2
n
)
!
4
n
n
!
π
=
(
2
n
−
1
)
!
!
2
n
π
=
(
n
−
1
2
n
)
n
!
π
Γ
(
1
2
−
n
)
=
(
−
4
)
n
n
!
(
2
n
)
!
π
=
(
−
2
)
n
(
2
n
−
1
)
!
!
π
=
π
(
−
1
/
2
n
)
n
!
{\displaystyle {\begin{aligned}\Gamma \left({\tfrac {1}{2}}+n\right)&={(2n)! \over 4^{n}n!}{\sqrt {\pi }}={\frac {(2n-1)!!}{2^{n}}}{\sqrt {\pi }}={\binom {n-{\frac {1}{2}}}{n}}n!{\sqrt {\pi }}\\[8pt]\Gamma \left({\tfrac {1}{2}}-n\right)&={(-4)^{n}n! \over (2n)!}{\sqrt {\pi }}={\frac {(-2)^{n}}{(2n-1)!!}}{\sqrt {\pi }}={\frac {\sqrt {\pi }}{{\binom {-1/2}{n}}n!}}\end{aligned}}}
(
2
n
−
1
)
!
!
=
(
2
n
−
1
)
(
2
n
−
3
)
⋯
(
3
)
(
1
)
{\displaystyle (2n-1)!!=(2n-1)(2n-3)\cdots (3)(1)}
が有理数である 他の個々の値に対する式を探すことによって、 という結果を一般化することは魅力的かもしれない。特に、 ガウスのディガンマ定理 によれば、 すべての有理数で密接に関連する ディガンマ関数 に対してそうすることが可能だからである。しかし、これらの数は、 基本関数でそれ自体で表現できることは知られていない。は 超越数 であり 、任意の 整数 および各分数に対して から 代数的に独立である ことが証明されている 。 [10] 一般に、ガンマ関数の値を計算するときは、数値近似で妥協しなければならない。
Γ
(
1
2
)
=
π
{\textstyle \Gamma \left({\frac {1}{2}}\right)={\sqrt {\pi }}}
Γ
(
r
)
{\displaystyle \Gamma (r)}
r
{\displaystyle r}
Γ
(
r
)
{\displaystyle \Gamma (r)}
Γ
(
n
+
r
)
{\displaystyle \Gamma (n+r)}
π
{\displaystyle \pi }
n
{\displaystyle n}
r
=
1
6
,
1
4
,
1
3
,
2
3
,
3
4
,
5
6
{\textstyle r={\frac {1}{6}},{\frac {1}{4}},{\frac {1}{3}},{\frac {2}{3}},{\frac {3}{4}},{\frac {5}{6}}}
ガンマ関数の導関数は、 ポリガンマ関数 ψ (0) ( z ) で表されます 。
正の整数 m に対して、ガンマ関数の導関数は次のように計算できます。
Γ
′
(
z
)
=
Γ
(
z
)
ψ
(
0
)
(
z
)
.
{\displaystyle \Gamma '(z)=\Gamma (z)\psi ^{(0)}(z).}
−2 − 2 i から6 + 2 i までの 複素平面におけるガンマ関数の引数を示す色
Γ
′
(
m
+
1
)
=
m
!
(
−
γ
+
∑
k
=
1
m
1
k
)
=
m
!
(
−
γ
+
H
(
m
)
)
,
{\displaystyle \Gamma '(m+1)=m!\left(-\gamma +\sum _{k=1}^{m}{\frac {1}{k}}\right)=m!\left(-\gamma +H(m)\right)\,,}
ここでH(m)はm次の 高調波数 、 γは オイラー・マスケロニ定数 です 。
ガンマ関数の 番目の導関数は、次のようになります。
( これ
は、ガンマ関数の積分形式を について微分し 、 積分記号 の下での微分 法を使用することで導出できます。)
ℜ
(
z
)
>
0
{\displaystyle \Re (z)>0}
n
{\displaystyle n}
d
n
d
z
n
Γ
(
z
)
=
∫
0
∞
t
z
−
1
e
−
t
(
log
t
)
n
d
t
.
{\displaystyle {\frac {d^{n}}{dz^{n}}}\Gamma (z)=\int _{0}^{\infty }t^{z-1}e^{-t}(\log t)^{n}\,dt.}
z
{\displaystyle z}
は リーマン ゼータ関数 、は -次の ベル多項式 である 恒等式を用いると
、 ガンマ関数の ローラン級数 展開が得られる [11]。
Γ
(
n
)
(
1
)
=
(
−
1
)
n
B
n
(
γ
,
1
!
ζ
(
2
)
,
…
,
(
n
−
1
)
!
ζ
(
n
)
)
{\displaystyle \Gamma ^{(n)}(1)=(-1)^{n}B_{n}(\gamma ,1!\zeta (2),\ldots ,(n-1)!\zeta (n))}
ζ
(
z
)
{\displaystyle \zeta (z)}
B
n
{\displaystyle B_{n}}
n
{\displaystyle n}
Γ
(
z
)
=
1
z
−
γ
+
1
2
(
γ
2
+
π
2
6
)
z
−
1
6
(
γ
3
+
γ
π
2
2
+
2
ζ
(
3
)
)
z
2
+
O
(
z
3
)
.
{\displaystyle \Gamma (z)={\frac {1}{z}}-\gamma +{\frac {1}{2}}\left(\gamma ^{2}+{\frac {\pi ^{2}}{6}}\right)z-{\frac {1}{6}}\left(\gamma ^{3}+{\frac {\gamma \pi ^{2}}{2}}+2\zeta (3)\right)z^{2}+O(z^{3}).}
不平等
正の実数に制限される場合、ガンマ関数は厳密に 対数的に凸な関数 です。この特性は、次の 3 つの同等な方法のいずれかで表すことができます。
任意の2つの正の実数 と 、および任意の に対して 、
x
1
{\displaystyle x_{1}}
x
2
{\displaystyle x_{2}}
t
∈
[
0
,
1
]
{\displaystyle t\in [0,1]}
Γ
(
t
x
1
+
(
1
−
t
)
x
2
)
≤
Γ
(
x
1
)
t
Γ
(
x
2
)
1
−
t
.
{\displaystyle \Gamma (tx_{1}+(1-t)x_{2})\leq \Gamma (x_{1})^{t}\Gamma (x_{2})^{1-t}.}
任意の2つの正の実数 とに対して 、 >
x
1
{\displaystyle x_{1}}
x
2
{\displaystyle x_{2}}
x
2
{\displaystyle x_{2}}
x
1
{\displaystyle x_{1}}
(
Γ
(
x
2
)
Γ
(
x
1
)
)
1
x
2
−
x
1
>
exp
(
Γ
′
(
x
1
)
Γ
(
x
1
)
)
.
{\displaystyle \left({\frac {\Gamma (x_{2})}{\Gamma (x_{1})}}\right)^{\frac {1}{x_{2}-x_{1}}}>\exp \left({\frac {\Gamma '(x_{1})}{\Gamma (x_{1})}}\right).}
任意の正の実数に対して 、
x
{\displaystyle x}
Γ
″
(
x
)
Γ
(
x
)
>
Γ
′
(
x
)
2
.
{\displaystyle \Gamma ''(x)\Gamma (x)>\Gamma '(x)^{2}.}
これらのステートメントの最後のものは、本質的には、 というステートメントと同じです。 ここで、は 1 次 ポリガンマ関数 です。したがって、ガンマ関数の対数凸性を証明するには、 が、正の実数 x に対して、正の項のみで構成される
級数表現を持つこと を観察するだけで十分です。
ψ
(
1
)
(
x
)
>
0
{\displaystyle \psi ^{(1)}(x)>0}
ψ
(
1
)
{\displaystyle \psi ^{(1)}}
ψ
(
1
)
{\displaystyle \psi ^{(1)}}
対数凸性と ジェンセンの不等式 を合わせると、任意の正の実数 とに対して 、
x
1
,
…
,
x
n
{\displaystyle x_{1},\ldots ,x_{n}}
a
1
,
…
,
a
n
{\displaystyle a_{1},\ldots ,a_{n}}
Γ
(
a
1
x
1
+
⋯
+
a
n
x
n
a
1
+
⋯
+
a
n
)
≤
(
Γ
(
x
1
)
a
1
⋯
Γ
(
x
n
)
a
n
)
1
a
1
+
⋯
+
a
n
.
{\displaystyle \Gamma \left({\frac {a_{1}x_{1}+\cdots +a_{n}x_{n}}{a_{1}+\cdots +a_{n}}}\right)\leq {\bigl (}\Gamma (x_{1})^{a_{1}}\cdots \Gamma (x_{n})^{a_{n}}{\bigr )}^{\frac {1}{a_{1}+\cdots +a_{n}}}.}
ガンマ関数の比にも上限がある。最もよく知られているのは ガウチの不等式 で、任意の正の実数 x と任意の s∈ (0, 1) に対して、
x
1
−
s
<
Γ
(
x
+
1
)
Γ
(
x
+
s
)
<
(
x
+
1
)
1
−
s
.
{\displaystyle x^{1-s}<{\frac {\Gamma (x+1)}{\Gamma (x+s)}}<\left(x+1\right)^{1-s}.}
複素平面におけるガンマ関数の表現。各点は の引数に応じて色分けされます 。 係数の等高線図 も表示されます。
z
{\displaystyle z}
Γ
(
z
)
{\displaystyle \Gamma (z)}
|
Γ
(
z
)
|
{\displaystyle |\Gamma (z)|}
複素ガンマ関数の絶対値の3次元プロット
増加する正の実変数に対する の挙動は スターリングの公式 で与えられ、
ここで記号は 漸近収束を意味します。つまり、両辺の比は極限 で 1 に収束します 。 [1] の 任意の固定値に対して、 この増加は指数関数的な増加よりも速くなります 。
Γ
(
x
)
{\displaystyle \Gamma (x)}
Γ
(
x
+
1
)
∼
2
π
x
(
x
e
)
x
,
{\displaystyle \Gamma (x+1)\sim {\sqrt {2\pi x}}\left({\frac {x}{e}}\right)^{x},}
∼
{\displaystyle \sim }
x
→
+
∞
{\textstyle x\to +\infty }
exp
(
β
x
)
{\displaystyle \exp(\beta x)}
β
{\displaystyle \beta }
漸近近似のもう一つの有用な限界は次の通り です。
x
→
+
∞
{\displaystyle x\to +\infty }
Γ
(
x
+
α
)
∼
Γ
(
x
)
x
α
,
α
∈
C
.
{\displaystyle {\Gamma (x+\alpha )}\sim {\Gamma (x)x^{\alpha }},\qquad \alpha \in \mathbb {C} .}
誤差項を無限積として表す場合、スターリングの公式を使用してガンマ関数を定義することができる: [12]
Γ
(
x
)
=
2
π
x
(
x
e
)
x
∏
n
=
0
∞
e
−
1
(
1
+
1
x
+
n
)
1
2
+
x
+
n
{\displaystyle \Gamma \left(x\right)={\sqrt {\frac {2\pi }{x}}}\left({\frac {x}{e}}\right)^{x}\prod _{n=0}^{\infty }e^{-1}\left(1+{\frac {1}{x+n}}\right)^{{\frac {1}{2}}+x+n}}
負の非整数値への拡張
ガンマ関数の主な定義である第二種オイラー積分は(実軸上で)正の引数に対してのみ有効であるが、負の引数を正の値にシフトすることにより、オイラーの反射公式
または基本的性質
のいずれかを使用して、
解析接続 [13] によって負の引数にその定義域を拡張することが
できる 。たとえば、
Γ
(
−
x
)
=
1
Γ
(
x
+
1
)
π
sin
(
π
(
x
+
1
)
)
,
{\displaystyle \Gamma (-x)={\frac {1}{\Gamma (x+1)}}{\frac {\pi }{\sin {\big (}\pi (x+1){\big )}}},}
Γ
(
−
x
)
:=
1
−
x
Γ
(
−
x
+
1
)
,
{\displaystyle \Gamma (-x):={\frac {1}{-x}}\Gamma (-x+1),}
x
∉
Z
{\displaystyle x\not \in \mathbb {Z} }
Γ
(
−
1
2
)
=
−
2
Γ
(
1
2
)
.
{\displaystyle \Gamma \left(-{\frac {1}{2}}\right)=-2\Gamma \left({\frac {1}{2}}\right).}
残留物
非正の場合には、 より複雑になります。オイラー積分は に対して収束しませんが 、 正の複素半平面で定義される関数は、負の半平面への一意の 解析接続を 持ちます。解析接続を見つける 1 つの方法は、正の引数に対してオイラー積分を使用し、漸化式 [1] を繰り返し適用してドメインを負の数に拡張することです。このとき、が正になるように
選択します。 が 整数 のいずれかに等しい 場合、分母の積はゼロになります。したがって、 ゼロ除算を 回避するには、ガンマ関数をそれらの点で未定義にする必要があります。これは、 非正の整数で 単純な極 を持つ 有理型関数 です。 [1]
z
{\displaystyle z}
ℜ
(
z
)
≤
0
{\displaystyle \Re (z)\leq 0}
Γ
(
z
)
=
Γ
(
z
+
n
+
1
)
z
(
z
+
1
)
⋯
(
z
+
n
)
,
{\displaystyle \Gamma (z)={\frac {\Gamma (z+n+1)}{z(z+1)\cdots (z+n)}},}
n
{\displaystyle n}
z
+
n
{\displaystyle z+n}
z
{\displaystyle z}
0
,
−
1
,
−
2
,
…
{\displaystyle 0,-1,-2,\ldots }
複素変数の 関数の場合 、 単極 におけるの 留 数 は次のように与えられます。
f
{\displaystyle f}
z
{\displaystyle z}
c
{\displaystyle c}
f
{\displaystyle f}
Res
(
f
,
c
)
=
lim
z
→
c
(
z
−
c
)
f
(
z
)
.
{\displaystyle \operatorname {Res} (f,c)=\lim _{z\to c}(z-c)f(z).}
単極 の場合 、漸化式は次のように書き直すことができます。
における分子は で
、分母は
です。したがって、これらの点におけるガンマ関数の留数は、次のようになります。 [14] ガンマ関数は実数線上のどこでもゼロではありませんが、 z → −∞ のときにゼロに任意に近づきます。 となる複素数は実際には存在しない ため 、 逆ガンマ関数 はでゼロとなる 整関数 です 。 [1]
z
=
−
n
{\displaystyle z=-n}
(
z
+
n
)
Γ
(
z
)
=
Γ
(
z
+
n
+
1
)
z
(
z
+
1
)
⋯
(
z
+
n
−
1
)
.
{\displaystyle (z+n)\Gamma (z)={\frac {\Gamma (z+n+1)}{z(z+1)\cdots (z+n-1)}}.}
z
=
−
n
,
{\displaystyle z=-n,}
Γ
(
z
+
n
+
1
)
=
Γ
(
1
)
=
1
{\displaystyle \Gamma (z+n+1)=\Gamma (1)=1}
z
(
z
+
1
)
⋯
(
z
+
n
−
1
)
=
−
n
(
1
−
n
)
⋯
(
n
−
1
−
n
)
=
(
−
1
)
n
n
!
.
{\displaystyle z(z+1)\cdots (z+n-1)=-n(1-n)\cdots (n-1-n)=(-1)^{n}n!.}
Res
(
Γ
,
−
n
)
=
(
−
1
)
n
n
!
.
{\displaystyle \operatorname {Res} (\Gamma ,-n)={\frac {(-1)^{n}}{n!}}.}
z
{\displaystyle z}
Γ
(
z
)
=
0
{\displaystyle \Gamma (z)=0}
1
Γ
(
z
)
{\textstyle {\frac {1}{\Gamma (z)}}}
z
=
0
,
−
1
,
−
2
,
…
{\displaystyle z=0,-1,-2,\ldots }
最小値と最大値
実数直線上では、ガンマ関数はz min ≈ +1.46163 21449 68362 34126 [15] で極小 値を持ち、 Γ( z min ) ≈ +0.88560 31944 10888 70027 の 値に達します。 [16]ガンマ関数はこの最小値のどちら側でも上昇します。 Γ( z − 0.5) = Γ( z + 0.5) の解は z = +1.5 であり 、共通の値は Γ(1) = Γ(2) = +1です。 Γ( z − 1) = Γ( z + 1) の正の解は z = φ ≈ +1.618 、 黄金比 であり 、共通の値は Γ( φ − 1) = Γ( φ + 1) = φ です。 ≈ +1.44922 96022 69896 60037 。 [17]
ガンマ関数は、非正の整数における極の符号が交互に変わる必要がある。これは 、との間の極の数 が奇数の場合は順方向回帰の積に奇数の負の因子が含まれ、極の数が偶数の場合は偶数が含まれるためである。 [14] 非正の整数間の実軸に沿ったガンマ関数の局所的極値の値は次のようになる。
z
{\displaystyle z}
z
+
n
{\displaystyle z+n}
Γ( −0.50408 30082 64455 40925... [18] ) = −3.54464 36111 55005 08912... 、
Γ( −1.57349 84731 62390 45877... [19] ) = 2.30240 72583 39680 13582... 、
Γ( −2.61072 08684 44144 65000... [20] ) = −0.88813 63584 01241 92009... 、
Γ( −3.63529 33664 36901 09783... [21] ) = 0.24512 75398 34366 25043... 、
Γ( −4.65323 77617 43142 44171... [22] ) = −0.05277 96395 87319 40076... など
積分表現
ガンマ関数を積分として表す公式は、第二種オイラー積分の他にもたくさんある。例えば、 z の実部が正のとき、 [23]
と [24] である。ここで、3つの積分は、それぞれ、オイラーの第二積分における、 [25] と [26]の
代入から導かれる。特に最後の積分は 、
半整数引数におけるガンマ関数と ガウス積分 との関係を明確にする。
Γ
(
z
)
=
∫
−
∞
∞
e
z
t
−
e
t
d
t
{\displaystyle \Gamma (z)=\int _{-\infty }^{\infty }e^{zt-e^{t}}\,dt}
Γ
(
z
)
=
∫
0
1
(
log
1
t
)
z
−
1
d
t
,
{\displaystyle \Gamma (z)=\int _{0}^{1}\left(\log {\frac {1}{t}}\right)^{z-1}\,dt,}
Γ
(
z
)
=
2
c
z
∫
0
∞
t
2
z
−
1
e
−
c
t
2
d
t
,
c
>
0
{\displaystyle \Gamma (z)=2c^{z}\int _{0}^{\infty }t^{2z-1}e^{-ct^{2}}\,dt\,,\;c>0}
t
=
e
−
x
{\displaystyle t=e^{-x}}
t
=
−
log
x
{\displaystyle t=-\log x}
t
=
c
x
2
{\displaystyle t=cx^{2}}
z
=
1
/
2
,
c
=
1
{\displaystyle z=1/2,\;c=1}
Γ
(
1
/
2
)
=
2
∫
0
∞
e
−
t
2
d
t
=
π
.
{\displaystyle \Gamma (1/2)=2\int _{0}^{\infty }e^{-t^{2}}\,dt={\sqrt {\pi }}\;.}
ビネのガンマ関数の第一積分公式は、 z の実部が正のとき、次のようになると述べている。 [27]
右辺の積分は ラプラス変換 として解釈することができる。つまり、
l
o
g
Γ
(
z
)
=
(
z
−
1
2
)
log
z
−
z
+
1
2
log
(
2
π
)
+
∫
0
∞
(
1
2
−
1
t
+
1
e
t
−
1
)
e
−
t
z
t
d
t
.
{\displaystyle \operatorname {log\Gamma } (z)=\left(z-{\frac {1}{2}}\right)\log z-z+{\frac {1}{2}}\log(2\pi )+\int _{0}^{\infty }\left({\frac {1}{2}}-{\frac {1}{t}}+{\frac {1}{e^{t}-1}}\right){\frac {e^{-tz}}{t}}\,dt.}
log
(
Γ
(
z
)
(
e
z
)
z
z
2
π
)
=
L
(
1
2
t
−
1
t
2
+
1
t
(
e
t
−
1
)
)
(
z
)
.
{\displaystyle \log \left(\Gamma (z)\left({\frac {e}{z}}\right)^{z}{\sqrt {\frac {z}{2\pi }}}\right)={\mathcal {L}}\left({\frac {1}{2t}}-{\frac {1}{t^{2}}}+{\frac {1}{t(e^{t}-1)}}\right)(z).}
ビネの第二積分公式によれば、 z の実部 が正のとき、次の式が成り立つ。 [28]
l
o
g
Γ
(
z
)
=
(
z
−
1
2
)
log
z
−
z
+
1
2
log
(
2
π
)
+
2
∫
0
∞
arctan
(
t
/
z
)
e
2
π
t
−
1
d
t
.
{\displaystyle \operatorname {log\Gamma } (z)=\left(z-{\frac {1}{2}}\right)\log z-z+{\frac {1}{2}}\log(2\pi )+2\int _{0}^{\infty }{\frac {\arctan(t/z)}{e^{2\pi t}-1}}\,dt.}
C を ハンケル路、つまり リーマン球面 上の点 ∞ で始まり点 ∞ で終わる経路 とし 、 その単位接線ベクトルは 経路の始めで −1 に収束し、終わりで 1に収束し、 0の周りの 巻数 1を持ち、 [0, ∞) を横切らないものとする。 [0, ∞) に沿って分枝切断を取り、 t が 負の実軸上にある ときに を実数として、 の 枝を固定する 。z は整数でないと仮定する。すると、ガンマ関数のハンケルの公式は次のようになる: [29]
ここで は と解釈される。反射公式から、 z が整数でない
場合は常に有効な、 密接に関連した式が得られる
。
log
(
−
t
)
{\displaystyle \log(-t)}
log
(
−
t
)
{\displaystyle \log(-t)}
Γ
(
z
)
=
−
1
2
i
sin
π
z
∫
C
(
−
t
)
z
−
1
e
−
t
d
t
,
{\displaystyle \Gamma (z)=-{\frac {1}{2i\sin \pi z}}\int _{C}(-t)^{z-1}e^{-t}\,dt,}
(
−
t
)
z
−
1
{\displaystyle (-t)^{z-1}}
exp
(
(
z
−
1
)
log
(
−
t
)
)
{\displaystyle \exp((z-1)\log(-t))}
1
Γ
(
z
)
=
i
2
π
∫
C
(
−
t
)
−
z
e
−
t
d
t
,
{\displaystyle {\frac {1}{\Gamma (z)}}={\frac {i}{2\pi }}\int _{C}(-t)^{-z}e^{-t}\,dt,}
連分数表現
ガンマ関数は2つの 連分数 の和で表すこともできる: [30] [31]
ここで 。
Γ
(
z
)
=
e
−
1
2
+
0
−
z
+
1
z
−
1
2
+
2
−
z
+
2
z
−
2
2
+
4
−
z
+
3
z
−
3
2
+
6
−
z
+
4
z
−
4
2
+
8
−
z
+
5
z
−
5
2
+
10
−
z
+
⋱
+
e
−
1
z
+
0
−
z
+
0
z
+
1
+
1
z
+
2
−
z
+
1
z
+
3
+
2
z
+
4
−
z
+
2
z
+
5
+
3
z
+
6
−
⋱
{\displaystyle {\begin{aligned}\Gamma (z)&={\cfrac {e^{-1}}{2+0-z+1{\cfrac {z-1}{2+2-z+2{\cfrac {z-2}{2+4-z+3{\cfrac {z-3}{2+6-z+4{\cfrac {z-4}{2+8-z+5{\cfrac {z-5}{2+10-z+\ddots }}}}}}}}}}}}\\&+\ {\cfrac {e^{-1}}{z+0-{\cfrac {z+0}{z+1+{\cfrac {1}{z+2-{\cfrac {z+1}{z+3+{\cfrac {2}{z+4-{\cfrac {z+2}{z+5+{\cfrac {3}{z+6-\ddots }}}}}}}}}}}}}}\end{aligned}}}
z
∈
C
{\displaystyle z\in \mathbb {C} }
フーリエ級数展開
ガンマ関数の対数には次のような フーリエ級数 展開があり
、これは長い間 エルンスト・クンマー が1847年に導出したと考えられていた。 [32] [33]しかし、イアロスラフ・ブラグーシンは、 カール・ヨハン・マルムステンが 1842年に初めてこの級数を導出したこと を発見した。 [34] [35]
0
<
z
<
1
:
{\displaystyle 0<z<1:}
l
o
g
Γ
(
z
)
=
(
1
2
−
z
)
(
γ
+
log
2
)
+
(
1
−
z
)
log
π
−
1
2
log
sin
(
π
z
)
+
1
π
∑
n
=
1
∞
log
n
n
sin
(
2
π
n
z
)
,
{\displaystyle \operatorname {log\Gamma } (z)=\left({\frac {1}{2}}-z\right)(\gamma +\log 2)+(1-z)\log \pi -{\frac {1}{2}}\log \sin(\pi z)+{\frac {1}{\pi }}\sum _{n=1}^{\infty }{\frac {\log n}{n}}\sin(2\pi nz),}
1840
年に ヨーゼフ・ルートヴィヒ・ラーベ は
、特に、
∫
a
a
+
1
log
Γ
(
z
)
d
z
=
1
2
log
2
π
+
a
log
a
−
a
,
a
>
0.
{\displaystyle \int _{a}^{a+1}\log \Gamma (z)\,dz={\tfrac {1}{2}}\log 2\pi +a\log a-a,\quad a>0.}
a
=
0
{\displaystyle a=0}
∫
0
1
log
Γ
(
z
)
d
z
=
1
2
log
2
π
.
{\displaystyle \int _{0}^{1}\log \Gamma (z)\,dz={\tfrac {1}{2}}\log 2\pi .}
後者は、上記の乗算公式の対数を取ることで導出でき、これは被積分関数のリーマン和の表現を与えます。 の極限を取ると、 式が得られます。
a
→
∞
{\displaystyle a\to \infty }
π関数
ガウス によって導入された別の表記法は 、ガンマ関数のシフトバージョンである -関数
です。
つまり、 すべての非負の整数 に対して となります 。
Π
{\displaystyle \Pi }
Π
(
z
)
=
Γ
(
z
+
1
)
=
z
Γ
(
z
)
=
∫
0
∞
e
−
t
t
z
d
t
,
{\displaystyle \Pi (z)=\Gamma (z+1)=z\Gamma (z)=\int _{0}^{\infty }e^{-t}t^{z}\,dt,}
Π
(
n
)
=
n
!
{\displaystyle \Pi (n)=n!}
n
{\displaystyle n}
π 関数を使用すると、反射式は次のようになります。
正規化された sinc 関数 を使用すると、乗法定理は次のようになります。
Π
(
z
)
Π
(
−
z
)
=
π
z
sin
(
π
z
)
=
1
sinc
(
z
)
{\displaystyle \Pi (z)\Pi (-z)={\frac {\pi z}{\sin(\pi z)}}={\frac {1}{\operatorname {sinc} (z)}}}
Π
(
z
m
)
Π
(
z
−
1
m
)
⋯
Π
(
z
−
m
+
1
m
)
=
(
2
π
)
m
−
1
2
m
−
z
−
1
2
Π
(
z
)
.
{\displaystyle \Pi \left({\frac {z}{m}}\right)\,\Pi \left({\frac {z-1}{m}}\right)\cdots \Pi \left({\frac {z-m+1}{m}}\right)=(2\pi )^{\frac {m-1}{2}}m^{-z-{\frac {1}{2}}}\Pi (z)\ .}
シフトされた 逆ガンマ関数は 、
完全 関数 と呼ばれることもあります 。
π
(
z
)
=
1
Π
(
z
)
,
{\textstyle \pi (z)={\frac {1}{\Pi (z)}}\ ,}
半径 r 1 , …, r nの n 楕円体 の体積は 次 のように表される。
V
n
(
r
1
,
…
,
r
n
)
=
π
n
2
Π
(
n
2
)
∏
k
=
1
n
r
k
.
{\displaystyle V_{n}(r_{1},\dotsc ,r_{n})={\frac {\pi ^{\frac {n}{2}}}{\Pi \left({\frac {n}{2}}\right)}}\prod _{k=1}^{n}r_{k}.}
他の機能との関係
ガンマ関数を定義する最初の積分では、積分の限界は固定されています。上側の 不完全ガンマ関数は、 積分の下限を可変にすることで得られます。 同様の下側の不完全ガンマ関数もあります。
Γ
(
z
,
x
)
=
∫
x
∞
t
z
−
1
e
−
t
d
t
.
{\displaystyle \Gamma (z,x)=\int _{x}^{\infty }t^{z-1}e^{-t}dt.}
ガンマ関数はオイラーの ベータ関数 と次の式で関係している。
B
(
z
1
,
z
2
)
=
∫
0
1
t
z
1
−
1
(
1
−
t
)
z
2
−
1
d
t
=
Γ
(
z
1
)
Γ
(
z
2
)
Γ
(
z
1
+
z
2
)
.
{\displaystyle \mathrm {B} (z_{1},z_{2})=\int _{0}^{1}t^{z_{1}-1}(1-t)^{z_{2}-1}\,dt={\frac {\Gamma (z_{1})\,\Gamma (z_{2})}{\Gamma (z_{1}+z_{2})}}.}
ガンマ関数の 対数微分 は ディガンマ関数 と呼ばれ、高次の微分は ポリガンマ関数 と呼ばれます。
有限体 または 有限環 上のガンマ関数の類似物は、 指数和 の一種である ガウス和 です 。
逆 ガンマ関数は 完全な関数 であり 、特定のトピックとして研究されてきました。
ガンマ関数は、 リーマンゼータ関数 との重要な関係にも現れます 。また 、 の場合にのみ有効な 次の式にも現れます 。
ζ
(
z
)
{\displaystyle \zeta (z)}
π
−
z
2
Γ
(
z
2
)
ζ
(
z
)
=
π
−
1
−
z
2
Γ
(
1
−
z
2
)
ζ
(
1
−
z
)
.
{\displaystyle \pi ^{-{\frac {z}{2}}}\;\Gamma \left({\frac {z}{2}}\right)\zeta (z)=\pi ^{-{\frac {1-z}{2}}}\;\Gamma \left({\frac {1-z}{2}}\right)\;\zeta (1-z).}
ζ
(
z
)
Γ
(
z
)
=
∫
0
∞
u
z
e
u
−
1
d
u
u
,
{\displaystyle \zeta (z)\Gamma (z)=\int _{0}^{\infty }{\frac {u^{z}}{e^{u}-1}}\,{\frac {du}{u}},}
ℜ
(
z
)
>
1
{\displaystyle \Re (z)>1}
ガンマ関数の対数は、レルヒによる次の式を満たします。 ここで 、は フルヴィッツのゼータ関数 、 はリーマンのゼータ関数、プライム ( ′ ) は最初の変数における微分を表します。
l
o
g
Γ
(
z
)
=
ζ
H
′
(
0
,
z
)
−
ζ
′
(
0
)
,
{\displaystyle \operatorname {log\Gamma } (z)=\zeta _{H}'(0,z)-\zeta '(0),}
ζ
H
{\displaystyle \zeta _{H}}
ζ
{\displaystyle \zeta }
ガンマ関数は 引き伸ばされた指数関数 と関連している。例えば、その関数のモーメントは
⟨
τ
n
⟩
≡
∫
0
∞
t
n
−
1
e
−
(
t
τ
)
β
d
t
=
τ
n
β
Γ
(
n
β
)
.
{\displaystyle \langle \tau ^{n}\rangle \equiv \int _{0}^{\infty }t^{n-1}\,e^{-\left({\frac {t}{\tau }}\right)^{\beta }}\,\mathrm {d} t={\frac {\tau ^{n}}{\beta }}\Gamma \left({n \over \beta }\right).}
特定の価値観
小数点以下20桁まで含めると、ガンマ関数の特定の値は次のようになる。
(これらの数値は OEIS
に掲載されている 。 [36] [37] [38] [39] [40 ] [41] ここで示す値は四捨五入ではなく切り捨てである。) 複素数値ガンマ関数は非正整数に対しては定義されていないが、これらの場合には リーマン球面上で値を ∞ として定義できる 。 逆ガンマ関数は これらの値(および 複素平面全体 )
で 明確に定義され 、 解析的である 。
Γ
(
−
3
2
)
=
4
π
3
≈
+
2.36327
18012
07354
70306
Γ
(
−
1
2
)
=
−
2
π
≈
−
3.54490
77018
11032
05459
Γ
(
1
2
)
=
π
≈
+
1.77245
38509
05516
02729
Γ
(
1
)
=
0
!
=
+
1
Γ
(
3
2
)
=
π
2
≈
+
0.88622
69254
52758
01364
Γ
(
2
)
=
1
!
=
+
1
Γ
(
5
2
)
=
3
π
4
≈
+
1.32934
03881
79137
02047
Γ
(
3
)
=
2
!
=
+
2
Γ
(
7
2
)
=
15
π
8
≈
+
3.32335
09704
47842
55118
Γ
(
4
)
=
3
!
=
+
6
{\displaystyle {\begin{array}{rcccl}\Gamma \left(-{\tfrac {3}{2}}\right)&=&{\tfrac {4{\sqrt {\pi }}}{3}}&\approx &+2.36327\,18012\,07354\,70306\\\Gamma \left(-{\tfrac {1}{2}}\right)&=&-2{\sqrt {\pi }}&\approx &-3.54490\,77018\,11032\,05459\\\Gamma \left({\tfrac {1}{2}}\right)&=&{\sqrt {\pi }}&\approx &+1.77245\,38509\,05516\,02729\\\Gamma (1)&=&0!&=&+1\\\Gamma \left({\tfrac {3}{2}}\right)&=&{\tfrac {\sqrt {\pi }}{2}}&\approx &+0.88622\,69254\,52758\,01364\\\Gamma (2)&=&1!&=&+1\\\Gamma \left({\tfrac {5}{2}}\right)&=&{\tfrac {3{\sqrt {\pi }}}{4}}&\approx &+1.32934\,03881\,79137\,02047\\\Gamma (3)&=&2!&=&+2\\\Gamma \left({\tfrac {7}{2}}\right)&=&{\tfrac {15{\sqrt {\pi }}}{8}}&\approx &+3.32335\,09704\,47842\,55118\\\Gamma (4)&=&3!&=&+6\end{array}}}
1
Γ
(
−
3
)
=
1
Γ
(
−
2
)
=
1
Γ
(
−
1
)
=
1
Γ
(
0
)
=
0.
{\displaystyle {\frac {1}{\Gamma (-3)}}={\frac {1}{\Gamma (-2)}}={\frac {1}{\Gamma (-1)}}={\frac {1}{\Gamma (0)}}=0.}
対数ガンマ関数
解析関数 logΓ( z )
ガンマ関数と階乗関数は、中程度の大きさの引数に対して非常に急速に増大するため、多くのコンピューティング環境には、 ガンマ関数の 自然対数 を返す関数が含まれており、プログラミング環境やスプレッドシートでは、この関数の lgamma名前が付け られていることが多い。この関数は増大がはるかに遅く、組み合わせ計算では、非常に大きな値の乗算や除算の代わりに、対数値の加算や減算を行うことができます。この関数は、次のように定義されることが多い [42]。 lngammagammaln
l
o
g
Γ
(
z
)
=
−
γ
z
−
log
z
+
∑
k
=
1
∞
[
z
k
−
log
(
1
+
z
k
)
]
.
{\displaystyle \operatorname {log\Gamma } (z)=-\gamma z-\log z+\sum _{k=1}^{\infty }\left[{\frac {z}{k}}-\log \left(1+{\frac {z}{k}}\right)\right].}
この関数の導関数であるディガンマ関数 も よく見られる。技術的および物理的応用の文脈では、例えば波の伝播では、関数方程式
l
o
g
Γ
(
z
)
=
l
o
g
Γ
(
z
+
1
)
−
log
z
{\displaystyle \operatorname {log\Gamma } (z)=\operatorname {log\Gamma } (z+1)-\log z}
−2 − 2iから2 + 2iまでの複素平面上の対数ガンマ関数(色付き)
は、隣接するストリップから z の幅 1 のストリップ内の関数値を決定できるため、よく使用されます。特に、大きな実部を持つ zの良好な近似から始めて、目的の z まで段階的に下げることができます 。 カール・フリードリヒ・ガウス の指示に従って、ロックテシェル (1922) は、 logΓ( z )に対して、大きな Re( z ) の近似を 提案しました 。
l
o
g
Γ
(
z
)
≈
(
z
−
1
2
)
log
z
−
z
+
1
2
log
(
2
π
)
.
{\displaystyle \operatorname {log\Gamma } (z)\approx (z-{\tfrac {1}{2}})\log z-z+{\tfrac {1}{2}}\log(2\pi ).}
これは、より小さな Re( z )を用いて、 z の logΓ( z ) を正確に近似するのに使用できる (PEBöhmer, 1939)。
l
o
g
Γ
(
z
−
m
)
=
l
o
g
Γ
(
z
)
−
∑
k
=
1
m
log
(
z
−
k
)
.
{\displaystyle \operatorname {log\Gamma } (z-m)=\operatorname {log\Gamma } (z)-\sum _{k=1}^{m}\log(z-k).}
スターリング近似に基づく
logΓ( z ) と Γ( z ) の漸近展開からのより多くの項を使用することで、より正確な近似を得ることができます。
Γ
(
z
)
∼
z
z
−
1
2
e
−
z
2
π
(
1
+
1
12
z
+
1
288
z
2
−
139
51
840
z
3
−
571
2
488
320
z
4
)
{\displaystyle \Gamma (z)\sim z^{z-{\frac {1}{2}}}e^{-z}{\sqrt {2\pi }}\left(1+{\frac {1}{12z}}+{\frac {1}{288z^{2}}}-{\frac {139}{51\,840z^{3}}}-{\frac {571}{2\,488\,320z^{4}}}\right)}
定数 | arg( z ) | < πで | z | → ∞ となる。( OEIS のシーケンス A001163 および A001164 を参照 。)
より「自然な」表現では:
l
o
g
Γ
(
z
)
=
z
log
z
−
z
−
1
2
log
z
+
1
2
log
2
π
+
1
12
z
−
1
360
z
3
+
1
1260
z
5
+
o
(
1
z
5
)
{\displaystyle \operatorname {log\Gamma } (z)=z\log z-z-{\tfrac {1}{2}}\log z+{\tfrac {1}{2}}\log 2\pi +{\frac {1}{12z}}-{\frac {1}{360z^{3}}}+{\frac {1}{1260z^{5}}}+o\left({\frac {1}{z^{5}}}\right)}
定数 | arg( z ) | < πで | z | → ∞ となる。( OEIS のシーケンス A046968 および A046969 を参照 。)
最後の展開における
z 1− k の k > 1 の項の係数は単に、 B k がベルヌーイ数
で ある 。
B
k
k
(
k
−
1
)
{\displaystyle {\frac {B_{k}}{k(k-1)}}}
ガンマ関数 には、1900年に シャルル・エルミートによって導出されたスターリング級数も存在する。 [43]
l
o
g
Γ
(
1
+
x
)
=
x
(
x
−
1
)
2
!
log
(
2
)
+
x
(
x
−
1
)
(
x
−
2
)
3
!
(
log
(
3
)
−
2
log
(
2
)
)
+
⋯
,
ℜ
(
x
)
>
0.
{\displaystyle \operatorname {log\Gamma } (1+x)={\frac {x(x-1)}{2!}}\log(2)+{\frac {x(x-1)(x-2)}{3!}}(\log(3)-2\log(2))+\cdots ,\quad \Re (x)>0.}
プロパティ
ボーア ・モレルプ定理は、 階乗関数を正の実数に拡張したすべての関数の中で、ガンマ関数だけが対数 凸 、つまりその 自然対数 が正の実軸上で 凸であることを述べています。別の特徴付けは、 ヴィーラント定理 によって与えられます 。
ガンマ関数は、次の式を同時に満たす唯一の関数である。
Γ
(
1
)
=
1
{\displaystyle \Gamma (1)=1}
、
Γ
(
z
+
1
)
=
z
Γ
(
z
)
{\displaystyle \Gamma (z+1)=z\Gamma (z)}
非正整数を除くすべての複素数に対して 、
z
{\displaystyle z}
整数 n に対して、 すべての複素数に対して 。 [1]
lim
n
→
∞
Γ
(
n
+
z
)
Γ
(
n
)
n
z
=
1
{\textstyle \lim _{n\to \infty }{\frac {\Gamma (n+z)}{\Gamma (n)\;n^{z}}}=1}
z
{\displaystyle z}
ある意味では、対数ガンマ関数はより自然な形式です。関数のいくつかの固有の属性がより明確になります。顕著な例は、 1 付近の
logΓ の テイラー級数です。ζ ( k )
は k における リーマン ゼータ関数 を表します 。
l
o
g
Γ
(
z
+
1
)
=
−
γ
z
+
∑
k
=
2
∞
ζ
(
k
)
k
(
−
z
)
k
∀
|
z
|
<
1
{\displaystyle \operatorname {log\Gamma } (z+1)=-\gamma z+\sum _{k=2}^{\infty }{\frac {\zeta (k)}{k}}\,(-z)^{k}\qquad \forall \;|z|<1}
したがって、次の特性を使用すると、
対数ガンマ関数の積分表現は次のようになります。
または、 z = 1を設定して γ の積分を取得すると 、 γ 項をその積分に置き換えて上記の式に組み込むと、次のようになります。
ζ
(
s
)
Γ
(
s
)
=
∫
0
∞
t
s
e
t
−
1
d
t
t
{\displaystyle \zeta (s)\Gamma (s)=\int _{0}^{\infty }{\frac {t^{s}}{e^{t}-1}}\,{\frac {dt}{t}}}
l
o
g
Γ
(
z
+
1
)
=
−
γ
z
+
∫
0
∞
e
−
z
t
−
1
+
z
t
t
(
e
t
−
1
)
d
t
{\displaystyle \operatorname {log\Gamma } (z+1)=-\gamma z+\int _{0}^{\infty }{\frac {e^{-zt}-1+zt}{t\left(e^{t}-1\right)}}\,dt}
l
o
g
Γ
(
z
+
1
)
=
∫
0
∞
e
−
z
t
−
z
e
−
t
−
1
+
z
t
(
e
t
−
1
)
d
t
.
{\displaystyle \operatorname {log\Gamma } (z+1)=\int _{0}^{\infty }{\frac {e^{-zt}-ze^{-t}-1+z}{t\left(e^{t}-1\right)}}\,dt\,.}
有理数z に対するガンマ関数の対数に対する特別な公式も存在します 。例えば、および が整数で 、 の場合 、 [44] この公式は、積分関数が非常に急速に減少するため、数値計算に使用されることがあります。
k
{\displaystyle k}
n
{\displaystyle n}
k
<
n
{\displaystyle k<n}
k
≠
n
/
2
,
{\displaystyle k\neq n/2\,,}
l
o
g
Γ
(
k
n
)
=
(
n
−
2
k
)
log
2
π
2
n
+
1
2
{
log
π
−
log
sin
π
k
n
}
+
1
π
∑
r
=
1
n
−
1
γ
+
log
r
r
⋅
sin
2
π
r
k
n
−
1
2
π
sin
2
π
k
n
⋅
∫
0
∞
e
−
n
x
⋅
log
x
cosh
x
−
cos
(
2
π
k
/
n
)
d
x
.
{\displaystyle {\begin{aligned}\operatorname {log\Gamma } \left({\frac {k}{n}}\right)={}&{\frac {\,(n-2k)\log 2\pi \,}{2n}}+{\frac {1}{2}}\left\{\,\log \pi -\log \sin {\frac {\pi k}{n}}\,\right\}+{\frac {1}{\pi }}\!\sum _{r=1}^{n-1}{\frac {\,\gamma +\log r\,}{r}}\cdot \sin {\frac {\,2\pi rk\,}{n}}\\&{}-{\frac {1}{2\pi }}\sin {\frac {2\pi k}{n}}\cdot \!\int _{0}^{\infty }\!\!{\frac {\,e^{-nx}\!\cdot \log x\,}{\,\cosh x-\cos(2\pi k/n)\,}}\,{\mathrm {d} }x.\end{aligned}}}
対数ガンマ積分
この積分は
Barnes G 関数 [45] [46]
を用いて表すことができる ( 証明については
Barnes G 関数 を参照)。
ここで Re( z )>−1 である。
∫
0
z
l
o
g
Γ
(
x
)
d
x
{\displaystyle \int _{0}^{z}\operatorname {log\Gamma } (x)\,dx}
∫
0
z
l
o
g
Γ
(
x
)
d
x
=
z
2
log
(
2
π
)
+
z
(
1
−
z
)
2
+
z
l
o
g
Γ
(
z
)
−
log
G
(
z
+
1
)
{\displaystyle \int _{0}^{z}\operatorname {log\Gamma } (x)\,dx={\frac {z}{2}}\log(2\pi )+{\frac {z(1-z)}{2}}+z\operatorname {log\Gamma } (z)-\log G(z+1)}
これはフルヴィッツゼータ関数 で表すこともできる : [47] [48]
∫
0
z
l
o
g
Γ
(
x
)
d
x
=
z
2
log
(
2
π
)
+
z
(
1
−
z
)
2
−
ζ
′
(
−
1
)
+
ζ
′
(
−
1
,
z
)
.
{\displaystyle \int _{0}^{z}\operatorname {log\Gamma } (x)\,dx={\frac {z}{2}}\log(2\pi )+{\frac {z(1-z)}{2}}-\zeta '(-1)+\zeta '(-1,z).}
となると
、これも ラーベの公式
の結果である 。O.エスピノサとV.モルは、の平方積分に対する同様の公式を導出した : [49]
ここで 、はである 。
z
=
1
{\displaystyle z=1}
∫
0
1
l
o
g
Γ
(
x
)
d
x
=
1
2
log
(
2
π
)
,
{\displaystyle \int _{0}^{1}\operatorname {log\Gamma } (x)\,dx={\frac {1}{2}}\log(2\pi ),}
l
o
g
Γ
{\displaystyle \operatorname {log\Gamma } }
∫
0
1
log
2
Γ
(
x
)
d
x
=
γ
2
12
+
π
2
48
+
1
3
γ
L
1
+
4
3
L
1
2
−
(
γ
+
2
L
1
)
ζ
′
(
2
)
π
2
+
ζ
′
′
(
2
)
2
π
2
,
{\displaystyle \int _{0}^{1}\log ^{2}\Gamma (x)dx={\frac {\gamma ^{2}}{12}}+{\frac {\pi ^{2}}{48}}+{\frac {1}{3}}\gamma L_{1}+{\frac {4}{3}}L_{1}^{2}-\left(\gamma +2L_{1}\right){\frac {\zeta ^{\prime }(2)}{\pi ^{2}}}+{\frac {\zeta ^{\prime \prime }(2)}{2\pi ^{2}}},}
L
1
{\displaystyle L_{1}}
1
2
log
(
2
π
)
{\displaystyle {\frac {1}{2}}\log(2\pi )}
DH Baileyと共著者 [50] は、Tornheim-Wittenゼータ関数とその導関数を用いての
評価 を行った
。
L
n
:=
∫
0
1
log
n
Γ
(
x
)
d
x
{\displaystyle L_{n}:=\int _{0}^{1}\log ^{n}\Gamma (x)\,dx}
n
=
1
,
2
{\displaystyle n=1,2}
さらに、 [51]
lim
n
→
∞
L
n
n
!
=
1.
{\displaystyle \lim _{n\to \infty }{\frac {L_{n}}{n!}}=1.}
近似値
ガンマ(青線)と階乗(青点)およびスターリング近似(赤線)の比較
ガンマ関数の複素数値は、 スターリング近似 または ランチョス近似 を使用して近似できます。
これは、 | z | が 無限大に近づく
につれて、近似値と真の値の比が極限で 1 に近づくという意味で 正確です。
Γ
(
z
)
∼
2
π
z
z
−
1
/
2
e
−
z
as
z
→
∞
in
|
arg
(
z
)
|
<
π
.
{\displaystyle \Gamma (z)\sim {\sqrt {2\pi }}z^{z-1/2}e^{-z}\quad {\hbox{as }}z\to \infty {\hbox{ in }}\left|\arg(z)\right|<\pi .}
ガンマ関数は、オイラー積分に 部分積分を 適用することで、 固定精度で計算できます。任意の正の数 x に対して、ガンマ関数は次のように表すことができます。
Re
(
z
)
∈
[
1
,
2
]
{\displaystyle \operatorname {Re} (z)\in [1,2]}
Γ
(
z
)
=
∫
0
x
e
−
t
t
z
d
t
t
+
∫
x
∞
e
−
t
t
z
d
t
t
=
x
z
e
−
x
∑
n
=
0
∞
x
n
z
(
z
+
1
)
⋯
(
z
+
n
)
+
∫
x
∞
e
−
t
t
z
d
t
t
.
{\displaystyle {\begin{aligned}\Gamma (z)&=\int _{0}^{x}e^{-t}t^{z}\,{\frac {dt}{t}}+\int _{x}^{\infty }e^{-t}t^{z}\,{\frac {dt}{t}}\\&=x^{z}e^{-x}\sum _{n=0}^{\infty }{\frac {x^{n}}{z(z+1)\cdots (z+n)}}+\int _{x}^{\infty }e^{-t}t^{z}\,{\frac {dt}{t}}.\end{aligned}}}
Re( z )∈[1,2] かつ のとき 、 最後の積分の絶対値は より小さくなります 。 を十分に大きい を選択すると、この最後の式は 任意の値 に対して より小さくすることができます 。したがって、ガンマ関数は 上記の級数を使用してビットの精度で評価できます。
x
≥
1
{\displaystyle x\geq 1}
(
x
+
1
)
e
−
x
{\displaystyle (x+1)e^{-x}}
x
{\displaystyle x}
2
−
N
{\displaystyle 2^{-N}}
N
{\displaystyle N}
N
{\displaystyle N}
EA Karatsubaは、任意の代数引数(有理数を含む)に対するオイラーガンマ関数の高速計算アルゴリズムを構築した。 [52] [53] [54]
の整数倍の引数の場合 1 / 24 、ガンマ関数は 算術幾何平均 反復法を使って素早く評価することもできる( ガンマ関数の特定の値 を参照)。 [55]
実践的な実装
正規分布 などの他の多くの関数とは異なり 、ガンマ関数に対して簡単に実装できる高速で正確な実装は 簡単には見つかりません。したがって、潜在的な解決策を調査する価値があります。速度が精度よりも重要な場合は、Online Wiley Libraryなどのインターネット検索で、公開されている の表を簡単に見つけることができます。このような表は 線形補間 で使用できます。3 次 補間を使用すると、 計算オーバーヘッドが増えますが、より高い精度が得られます。 表は通常、1から2までの引数値に対して公開されているため、プロパティを 使用して、すべての実数値 とを 範囲 に迅速かつ簡単に変換し 、1から2までの表の値のみ を使用する必要があります。 [56]
Γ
(
z
)
{\displaystyle \Gamma (z)}
Γ
(
z
)
{\displaystyle \Gamma (z)}
Γ
(
z
)
{\displaystyle \Gamma (z)}
Γ
(
z
+
1
)
=
z
Γ
(
z
)
{\displaystyle \Gamma (z+1)=z\ \Gamma (z)}
z
<
1
{\displaystyle z<1}
z
>
2
{\displaystyle z>2}
1
≤
z
≤
2
{\displaystyle 1\leq z\leq 2}
z
{\displaystyle z}
補間テーブルが望ましくない場合は、上記の Lanczos 近似が、一般的に使用される小さな z 値に対して 1 ~ 2 桁の精度でうまく機能します。Lanczos 近似の精度が十分でない場合は、 ガンマ関数の Stirling の式 を使用できます。
アプリケーション
ある著者はガンマ関数を「おそらく最も一般的な特殊関数、あるいは最も「特殊」でない関数である。他の超越関数は、多くの専門的な数学的話題を避けることで、それらのいくつかを避けることができると考えられるため、「特殊」と呼ばれている。一方、ガンマ関数 Γ( z )は 避けるのが最も難しい」と述べている。 [57]
統合の問題
ガンマ関数は 量子物理学 、 天体物理学 、 流体力学 など多様な分野で応用されています。 [58] ガンマ関数に基づいて定式化されたガンマ分布は、統計学において、地震の発生間隔など、さまざまなプロセスをモデル化するために使用され ます 。 [ 59 ]
このような状況でガンマ関数が有用である主な理由は、 時間や空間で指数関数的に減少する過程を記述するタイプの表現が広く普及していることである。そのような表現の積分は、基本的な解が存在しない場合にガンマ関数で解くことができることがある。例えば、 f がべき乗関数で gが 線形関数である場合、変数を単純に変更するだけで、 評価は次のようになる。
f
(
t
)
e
−
g
(
t
)
{\displaystyle f(t)e^{-g(t)}}
u
:=
a
⋅
t
{\displaystyle u:=a\cdot t}
∫
0
∞
t
b
e
−
a
t
d
t
=
1
a
b
∫
0
∞
u
b
e
−
u
d
(
u
a
)
=
Γ
(
b
+
1
)
a
b
+
1
.
{\displaystyle \int _{0}^{\infty }t^{b}e^{-at}\,dt={\frac {1}{a^{b}}}\int _{0}^{\infty }u^{b}e^{-u}d\left({\frac {u}{a}}\right)={\frac {\Gamma (b+1)}{a^{b+1}}}.}
積分が正の実数直線全体に沿って実行されるという事実は、ガンマ関数が無期限に続く時間依存プロセスの累積を記述していることを意味している可能性があり、またはその値は無限空間内の分布の合計である可能性があります。
もちろん、有限過程の累積を記述するために0 と ∞ 以外の積分極限を取ることがしばしば有用であり、その場合、通常のガンマ関数はもはや解ではなくなり、その解は 不完全ガンマ関数と呼ばれる。(正の実数直線全体にわたって積分することによって得られる通常のガンマ関数は、対比のために 完全ガンマ関数 と呼ばれることもある 。)
指数関数的に減少する関数の重要なカテゴリは、誤差関数
などの ガウス関数
とその積分です。これらの関数とガンマ関数の間には多くの相互関係があり、特に、 を評価することによって得られる 係数は、 誤差関数および 正規分布 の正規化係数に見られる係数と「同じ」です。
a
e
−
(
x
−
b
)
2
c
2
{\displaystyle ae^{-{\frac {(x-b)^{2}}{c^{2}}}}}
π
{\displaystyle {\sqrt {\pi }}}
Γ
(
1
2
)
{\textstyle \Gamma \left({\frac {1}{2}}\right)}
これまで説明した積分は 超越関数 を含んでいますが、ガンマ関数は純粋に代数的な関数の積分からも生じます。特に、代数方程式によって定義される曲線 である 楕円 と レムニスケートの 弧の長さは 、特殊な場合にはガンマ関数で評価できる 楕円積分 によって与えられます。ガンマ関数は、 n 次元の 超球 の 「体積」と「面積」を計算する ためにも使用できます。
積の計算
ガンマ関数の階乗積を一般化する能力は、数学の多くの分野、例えば 組合せ論 、さらには 確率論や べき級数 の計算などの分野への応用に直ちにつながります。連続する整数の積を含む多くの式は、階乗の組み合わせとして表すことができます。最も重要な例は、おそらく 二項係数 です 。たとえば、 | z | < 1 である任意の複素数 z と nに対して、 n が 負でない整数である
場合の二項係数に非常によく似た式を 書くことができます
。
(
1
+
z
)
n
=
∑
k
=
0
∞
Γ
(
n
+
1
)
k
!
Γ
(
n
−
k
+
1
)
z
k
,
{\displaystyle (1+z)^{n}=\sum _{k=0}^{\infty }{\frac {\Gamma (n+1)}{k!\Gamma (n-k+1)}}z^{k},}
(
1
+
z
)
n
=
∑
k
=
0
n
n
!
k
!
(
n
−
k
)
!
z
k
=
∑
k
=
0
n
(
n
k
)
z
k
.
{\displaystyle (1+z)^{n}=\sum _{k=0}^{n}{\frac {n!}{k!(n-k)!}}z^{k}=\sum _{k=0}^{n}{\binom {n}{k}}z^{k}.}
二項係数の例は、ガンマ関数を負の数に拡張した場合の特性が自然である理由を示しています。二項係数は、 n 個の要素 のセットから k個の要素を選択する方法の数を示します。 k > n の場合、もちろん方法はありません。 k > n の場合、 ( n − k )! は 負の整数の階乗であり、階乗のガンマ関数の定義を使用すると無限大になります。つまり、無限大で割ると、期待値は 0 になります。
階乗をガンマ関数に置き換えることで、このような式を複素数に拡張することができます。一般に、この方法は、各因数がインデックス変数の 有理関数 である積であれば、有理関数を線形式に因数分解することで適用できます。P と Qが それぞれ 根 p 1 、…、 p m と q 1 、… 、
q nを 持つ次数 m と nの単項多項式である場合、
∏
i
=
a
b
P
(
i
)
Q
(
i
)
=
(
∏
j
=
1
m
Γ
(
b
−
p
j
+
1
)
Γ
(
a
−
p
j
)
)
(
∏
k
=
1
n
Γ
(
a
−
q
k
)
Γ
(
b
−
q
k
+
1
)
)
.
{\displaystyle \prod _{i=a}^{b}{\frac {P(i)}{Q(i)}}=\left(\prod _{j=1}^{m}{\frac {\Gamma (b-p_{j}+1)}{\Gamma (a-p_{j})}}\right)\left(\prod _{k=1}^{n}{\frac {\Gamma (a-q_{k})}{\Gamma (b-q_{k}+1)}}\right).}
ガンマ関数を数値的に計算する方法があれば、そのような積の数値を計算するのは非常に簡単です。右辺のガンマ関数の数は多項式の次数のみに依存するため、 b − a が 5 に等しいか 10 5 に等しいかは問題ではありません。適切な制限を取ることで、左辺の積にゼロや極が含まれている場合でも方程式が成立するようにすることができます。
極限をとることで、無限数の因数を持つ特定の有理積もガンマ関数で評価できます。 ワイエルシュトラスの因数分解定理 により、解析関数は無限積として表すことができ、これらはガンマ関数の有限積または商として表されることもあります。すでに 1 つの顕著な例を見ました。反射式は、基本的に正弦関数を 2 つのガンマ関数の積として表します。この式から始めて、指数関数だけでなく、すべての三角関数と双曲関数をガンマ関数で表すことができます。
超幾何 関数とその特殊なケースを含むさらに多くの関数は、 メリン・バーンズ積分 と呼ばれるガンマ関数の 積 と商の複素積分によって表すことができます 。
解析的数論
ガンマ関数の応用として、リーマンゼータ関数 の研究があります 。リーマンゼータ関数の基本的な性質は、その 関数方程式 です。
Γ
(
s
2
)
ζ
(
s
)
π
−
s
2
=
Γ
(
1
−
s
2
)
ζ
(
1
−
s
)
π
−
1
−
s
2
.
{\displaystyle \Gamma \left({\frac {s}{2}}\right)\zeta (s)\pi ^{-{\frac {s}{2}}}=\Gamma \left({\frac {1-s}{2}}\right)\zeta (1-s)\pi ^{-{\frac {1-s}{2}}}.}
とりわけ、これは ゼータ関数の複素平面上の有理型関数への 解析接続の明示的な形式を提供し、ゼータ関数が実数直線上に無限に多くのいわゆる「自明な」零点を持つことを即座に証明する。Borwein らは この公式を「数学における最も美しい発見の1つ」と呼んでいる。 [60] この称号のもう1つの候補は、
ζ
(
s
)
Γ
(
s
)
=
∫
0
∞
t
s
e
t
−
1
d
t
t
.
{\displaystyle \zeta (s)\;\Gamma (s)=\int _{0}^{\infty }{\frac {t^{s}}{e^{t}-1}}\,{\frac {dt}{t}}.}
どちらの公式も、 ベルンハルト・リーマンが 1859 年に発表した画期的な論文「与えられた絶対値より小さい素数の個数 について 」で導き出したもので、数学的解析のツールを使用して素数を 研究する数学の分野である 解析的数論 の発展における画期的な出来事の 1 つです 。
歴史
ガンマ関数は、歴史上最も著名な数学者の何人かの関心を集めてきました。1963年に ショーヴネ賞を受賞した フィリップ J. デイビス の論文に記録されているガンマ関数の歴史は 、18 世紀以降の数学の主要な発展の多くを反映しています。デイビスの言葉を借りれば、「各世代はガンマ関数について何か興味深いことを発見してきました。おそらく次の世代もそうするでしょう。」 [1]
18世紀: オイラーとスターリング
ダニエル・ベルヌーイから クリスチャン・ゴールドバッハ への手紙 、1729年10月6日
階乗を非整数引数に拡張する問題は、 1720年代に ダニエル・ベルヌーイ と クリスチャン・ゴールドバッハ によって初めて検討されたようです。特に、1729年10月6日のベルヌーイからゴールドバッハへの手紙の中で、ベルヌーイは、負の整数以外
の実数値 x
に対して明確に定義された 積表現 [61]を導入しました。
x
!
=
lim
n
→
∞
(
n
+
1
+
x
2
)
x
−
1
∏
k
=
1
n
k
+
1
k
+
x
{\displaystyle x!=\lim _{n\to \infty }\left(n+1+{\frac {x}{2}}\right)^{x-1}\prod _{k=1}^{n}{\frac {k+1}{k+x}}}
レオンハルト・オイラーは 後に2つの異なる定義を与えた。最初のものは彼の積分ではなく、 負の整数以外の
すべての複素数 n に対して明確に定義された無限積 であり、
彼は1729年10月13日の手紙でゴールドバッハにそのことを伝えた。彼は1730年1月8日に再びゴールドバッハに手紙を書き、
複素数 n の実部が厳密に−1 より大きい場合 (つまり、 )に有効な積分表現を発見したことを発表した。変数変換 t = −ln s により、これはよく知られたオイラー積分になる。オイラーは、 1729年11月28日に サンクトペテルブルクアカデミー に提出した論文「超越数列、つまり一般項を代数的に与えることができないものについて」でその結果を発表しました。 [62] オイラーはさらに、反射公式を含むガンマ関数の重要な機能特性のいくつかを発見しました。
n
!
=
∏
k
=
1
∞
(
1
+
1
k
)
n
1
+
n
k
,
{\displaystyle n!=\prod _{k=1}^{\infty }{\frac {\left(1+{\frac {1}{k}}\right)^{n}}{1+{\frac {n}{k}}}}\,,}
n
!
=
∫
0
1
(
−
log
s
)
n
d
s
,
{\displaystyle n!=\int _{0}^{1}(-\log s)^{n}\,ds\,,}
ℜ
(
n
)
>
−
1
{\displaystyle \Re (n)>-1}
オイラーと同時代のジェームズ・スターリング も階乗の連続的な表現を見つけようとし、現在 スターリングの公式として知られるものを考案しました。スターリングの公式は非整数についても n ! の良好な推定値を与えますが、正確な値を与えるものではありません。この誤りを修正する彼の公式の拡張は、スターリング自身と ジャック・フィリップ・マリー・ビネ によって与えられました 。
19世紀: ガウス、ヴァイエルシュトラス、ルジャンドル
オイラーの論文の最初のページ
カール・フリードリヒ・ガウスは オイラーの積を と書き直し
、この公式を使ってガンマ関数の新しい性質を発見した。オイラーは複素変数理論の先駆者であったが、複素数の階乗についてはガウスが初めて考慮したようには見えなかった。 [63]ガウスはガンマ関数の 乗法定理 も証明し、ガンマ関数と 楕円積分 の関係を調査した 。
Γ
(
z
)
=
lim
m
→
∞
m
z
m
!
z
(
z
+
1
)
(
z
+
2
)
⋯
(
z
+
m
)
{\displaystyle \Gamma (z)=\lim _{m\to \infty }{\frac {m^{z}m!}{z(z+1)(z+2)\cdots (z+m)}}}
カール・ワイエルシュトラスは、 さらに別の積表現から始めて、
複素解析 におけるガンマ関数の役割を確立しました。
ここで、 γ はオイラー・マスケローニ定数 です。ワイエルシュトラスはもともと、この積を の積として書きました。
Γ
(
z
)
=
e
−
γ
z
z
∏
k
=
1
∞
(
1
+
z
k
)
−
1
e
z
k
,
{\displaystyle \Gamma (z)={\frac {e^{-\gamma z}}{z}}\prod _{k=1}^{\infty }\left(1+{\frac {z}{k}}\right)^{-1}e^{\frac {z}{k}},}
1 / Γ 、その場合、関数の極ではなく零点に作用します。この結果に触発されて、彼は ワイエルシュトラスの因数分解定理 として知られるものを証明しました。これは、任意の関数全体が複素平面上の零点の積として表されるというものであり、 代数の基本定理 の一般化です。
ガンマ関数という名前と記号 Γ は、1811 年頃に アドリアン マリー ルジャンドル によって導入されました 。ルジャンドルは、オイラーの積分定義を現代的な形に書き直しました。記号はギリシャ語の大文字のガンマですが、関数名を「ガンマ関数」と書くべきか「ガンマ関数」と書くべきかについての標準は認められていません (単に「 Γ 関数」と書く著者もいます)。ガウスによる代替の「π 関数」表記 Π( z ) = z ! は、古い文献で時々見られますが、ルジャンドルの表記は現代の著作では主流です。
なぜ「通常の階乗」とガンマ関数を別々の記号で区別するのか、特にガンマ関数を単に「 Γ( n ) = n ! 」ではなく Γ( n + 1) = n !に正規化する必要があるのか、という疑問は当然生じる。指数 x n の表記が 整数から複素数 x z にそのまま一般化されていることを考えてみよう。ルジャンドルが正規化を行った動機は不明であり、扱いにくいと批判する者もいる(たとえば、20世紀の数学者 コルネリウス・ランチョスは これを「合理性を欠く」と呼び、代わりに z ! を 使用した)。 [64] ルジャンドルの正規化はいくつかの式を単純化するが、他の式を複雑にする。現代的な観点から見ると、ガンマ関数のルジャンドル正規化は、 リー群 R + 上の ハール測度 に関して、加法 特性 e − x と乗法特性 x z の積分である。したがって、この正規化により、ガンマ関数がガウス和 の連続的な類似物であることがより明確になる 。 [65]
d
x
x
{\textstyle {\frac {dx}{x}}}
19世紀~20世紀: ガンマ関数の特徴づけ
ガンマ関数には多数の定義が与えられていることに多少問題がある。それらは同じ関数を記述しているにもかかわらず、その等価性を証明するのは簡単ではない。スターリングは彼の拡張された公式がオイラーのガンマ関数に正確に対応することを証明しなかった。証明は 1900年に シャルル・エルミートによって初めて与えられた。 [66] 各公式に特化した証明を見つける代わりに、ガンマ関数を識別する一般的な方法があることが望ましい。
同等性を証明する 1 つの方法は、ガンマ関数を特徴付ける 微分方程式 を見つけることです。応用数学のほとんどの特殊関数は、微分方程式の解として現れ、その解は一意です。しかし、ガンマ関数は単純な微分方程式を満たさないようです。 オットー ヘルダーは 1887 年に、ガンマ関数が少なくとも 代数 微分方程式を満たさないことを証明しました。これは、そのような方程式の解がガンマ関数の漸化式を満たすことができず、ガンマ関数が 超越的超越関数になること を示しました。この結果は ヘルダーの定理 として知られています 。
ガンマ関数の明確かつ一般的に適用可能な特徴付けは、1922 年まで示されていませんでした。その後、 ハラルド ボーア と ヨハネス モレルプが、 ボーア - モレルプ定理 として知られる定理を証明しました 。これは、ガンマ関数は、正の z に対して正かつ 対数的に凸で あり、1 での値が 1 である (関数の対数が凸である場合、その関数は対数的に凸である) 階乗再帰関係の唯一の解であるという定理です。別の特徴付けは、ヴィーラント定理 によって与えられます 。
ボーア・モレルプ定理は、ガンマ関数を定義するために使用されるさまざまな式のいずれに対しても、対数凸性を証明するのが比較的簡単なので便利です。さらに、ガンマ関数を特定の式で定義する代わりに、ボーア・モレルプ定理の条件を定義として選択し、ガンマ関数の研究の出発点として、条件を満たす任意の式を選択できます。このアプローチは、ブルバキ グループ によって使用されました。
ボーウェイン とコーレスはガンマ関数に関する3世紀にわたる研究をレビューする。 [67]
参照表とソフトウェア
ガンマ関数は、現代のコンピュータ(プログラム可能なポケット電卓さえも)を使えば、数学的に単純な関数とほぼ同じくらい簡単に計算できるが、もちろん昔からそうだったわけではない。20世紀半ばまで、数学者は手作りの表に頼っていた。ガンマ関数の場合、1813年にガウスが計算した表や、1825年にルジャンドルが計算した表が有名である。 [68]
Jahnke とEmde 著 『Tables of Higher Functions』 [de] より、複素ガンマ関数の絶対値を手描きしたグラフ 。
ガンマ関数の複素数値の表と手書きのグラフは、 1909年にドイツで最初に出版された ヤーンケ とエムデ [de] による 「関数表、数式、曲線」 に掲載されている。マイケル・ベリー によると、「複素平面におけるガンマ関数の極を示す3次元グラフがJ&Eに掲載され、ほぼ象徴的な地位を獲得した。」 [69]
1930年代まで、ガンマ関数の実数値以外の実用的なニーズはほとんどありませんでしたが、その年に複素ガンマ関数の応用が理論物理学で発見されました。1950年代に電子計算機が表の作成に利用できるようになると、需要を満たすために複素ガンマ関数の詳細な表がいくつか出版されました。その中には、米国 国立標準技術局 の12桁の精度の表も含まれていました。 [1]
Janhke と Emde による有名な複素グラフ (Tables of Functions with Formulas and Curves、第 4 版、Dover、1945 年) の再現。ガンマ関数の -4.5 - 2.5i から 4.5 + 2.5i までのグラフです。
ガンマ関数とその対数の倍精度浮動小数点実装は、 TK Solver 、 Matlab 、 GNU Octave 、 GNU Scientific Library など、ほとんどの科学計算ソフトウェアと特殊関数ライブラリで利用できるようになりました。ガンマ関数は C 標準ライブラリ( math.h )にも追加されました。任意精度実装は、 Mathematica や Maple など、ほとんどの コンピュータ代数システム で利用できます。PARI /GP 、 MPFR 、MPFUNには、無料の任意精度実装が含まれています。Windows Calculator や GNOME Calculatorなどの一部の ソフトウェア計算機 では、入力 x が整数値以外の場合、 階乗関数はΓ( x + 1)を返します。 [70] [71]
参照
注記
^ abcdefghijklmn Davis, PJ (1959). 「レオンハルト・オイラーの積分:ガンマ関数の歴史的プロファイル」. American Mathematical Monthly . 66 (10): 849–869. doi :10.2307/2309786. JSTOR 2309786. 2012年11月7日時点のオリジナルよりアーカイブ 。 2016年 12月3日 閲覧。
^ 「ガンマ関数は誤って定義されていますか? または: アダマール対オイラー — より優れたガンマ関数を見つけたのは誰か?」
^ ビールズ、リチャード、ウォン、ロデリック(2010)。特殊機能:大学院テキスト。ケンブリッジ大学出版局。p. 28。ISBN 978-1-139-49043-6 。 28ページの抜粋
^ Ross, Clay C. (2013). 微分方程式: Mathematica による入門 (イラスト入り) Springer Science & Business Media. p. 293. ISBN 978-1-4757-3949-7 。 293ページの式G.2
^ Kingman, JFC (1961). 「正行列の凸性特性」. 季刊数学誌 . 12 (1): 283–284. Bibcode :1961QJMat..12..283K. doi :10.1093/qmath/12.1.283.
^ Weisstein, Eric W. 「ボーア・モレルプ定理」 。MathWorld 。
^ Davis, Philip. 「レオンハルト・オイラーの積分 : ガンマ関数の歴史的概要」 (PDF) 。maa.org 。
^ マルコ、ボンヴィーニ (2010 年 10 月 9 日)。 「ガンマ関数」 (PDF) 。 Roma1.infn.it 。
^ Askey, RA ; Roy, R. (2010)、「Series Expansions」、 Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.)、 NIST Handbook of Mathematical Functions 、Cambridge University Press、 ISBN 978-0-521-19225-5 、 MR 2723248 。
^ Waldschmidt, M. (2006). 「周期の超越: 最新技術」 (PDF) . Pure Appl. Math. Quart . 2 (2): 435–463. doi : 10.4310/pamq.2006.v2.n2.a3 . 2006年5月6日時点のオリジナルよりアーカイブ (PDF) 。
^ 「$z=0$付近のガンマ関数のローラン展開を得るにはどうすればいいですか?」 Mathematics Stack Exchange 。 2022年 8月17日 閲覧 。
^ Artin, Emil (2015). ガンマ関数 . ドーバー. p. 24.
^ オールダム、キース、マイランド、ヤン、スパニアー、ジェローム (2010)。関数 の アトラス (第 2 版)。第 43 章: Springer Science & Business Media。ISBN 9780387488073 。 {{cite book}}: CS1 maint: location (link)
^ ab ワイスタイン、エリック W. 「ガンマ関数」。 マスワールド 。
^ Sloane, N. J. A. (編)。「シーケンス A030169 (y = Gamma(x) が最小となるような実数 x の 10 進展開)」。 整数シーケンスのオンライン百科事典 。OEIS Foundation。
^ Sloane, N. J. A. (編)。「シーケンス A030171 (y = Gamma(x) が最小となるような実数 y の 10 進展開)」。 整数シーケンスのオンライン百科事典 。OEIS Foundation。
^ Sloane, N. J. A. (編)。「シーケンス A178840 (黄金比の階乗の 10 進展開)」。 整数シーケンスのオンライン百科事典 。OEIS 財団。
^ Sloane, N. J. A. (編)。「数列 A175472 (区間 [ -1,0] におけるガンマ関数の極大値の横座標の絶対値の 10 進展開)」。 整数数列のオンライン百科事典 。OEIS Foundation。
^ Sloane, N. J. A. (編)。「シーケンス A175473 (区間 [ -2,-1] におけるガンマ関数の極小値の横座標の絶対値の 10 進展開)」。 整数シーケンスのオンライン百科事典 。OEIS Foundation。
^ Sloane, N. J. A. (編)。「シーケンス A175474 (区間 [ -3,-2] におけるガンマ関数の極大値の横座標の絶対値の 10 進展開)」。 オンライン整数シーケンス百科事典 。OEIS Foundation。
^ Sloane, N. J. A. (編)。「シーケンス A256681 (区間 [-4,-3] におけるガンマ関数の極小値の [負の] 横座標の 10 進展開)」。 整数シーケンスのオンライン百科事典 。OEIS Foundation。
^ Sloane, N. J. A. (編)。「シーケンス A256682 (区間 [-5,-4] におけるガンマ関数の極大値の [負の] 横座標の 10 進展開)」。 オンライン整数シーケンス百科事典 。OEIS Foundation。
^ Gradshteyn, IS; Ryzhik, IM (2007). 積分、級数、積の表 (第 7 版). Academic Press. p. 893. ISBN 978-0-12-373637-6 。
^ WhittakerとWatson、12.2例1。
^ Detlef, Gronau. 「ガンマ関数がなぜそうなっているのか?」 (PDF) . Imsc.uni-graz.at .
^ Pascal Sebah、Xavier Gourdon。「ガンマ関数入門」 (PDF) 。 数値計算 。 2023年1月30日時点の オリジナル (PDF)からアーカイブ。 2023年 1月30日 閲覧 。
^ WhittakerとWatson、12.31。
^ WhittakerとWatson、12.32。
^ WhittakerとWatson、12.22。
^ 「指数積分E: 連分数表現(式06.34.10.0005)」。
^ 「指数積分E: 連分数表現(式06.34.10.0003)」。
^ ベイトマン、ハリー; エルデーリ、アーサー (1955)。 高等超越関数 。 マグロウヒル。OCLC 627135。
^ Srivastava, HM; Choi, J. (2001). ゼータと関連関数に関連するシリーズ 。オランダ:Kluwer Academic。ISBN 0-7923-7054-6 。
^ Blagouchine, Iaroslav V. (2014). 「マルムステン積分の再発見、等高線積分法による評価、および関連する結果」 Ramanujan J. 35 ( 1): 21–110. doi :10.1007/s11139-013-9528-5. S2CID 120943474.
^ Blagouchine, Iaroslav V. (2016). 「マルムステン積分の再発見、等高線積分法による評価、および関連する結果」に対する訂正と補遺 「 ラマヌジャンJ.42 ( 3 ):777–781.doi : 10.1007 / s11139-015-9763-z.S2CID 125198685 .
^ Sloane, N. J. A. (編)。「数列 A245886 (ガンマ(-3/2) の 10 進展開。ガンマはオイラーのガンマ関数)」。 整数数列のオンライン百科事典 。OEIS 財団。
^ Sloane, N. J. A. (編)。「シーケンス A019707 (sqrt(Pi)/5 の 10 進展開)」。 整数シーケンスのオンライン百科事典 。OEIS Foundation。
^ Sloane, N. J. A. (編)。「シーケンス A002161 (円周率の平方根の 10 進展開)」。 整数シーケンスのオンライン百科事典 。OEIS Foundation。
^ Sloane, N. J. A. (編)。「シーケンス A019704 (sqrt(Pi)/2 の 10 進展開)」。 整数シーケンスのオンライン百科事典 。OEIS Foundation。
^ Sloane, N. J. A. (編)。「数列 A245884 (ガンマ(5/2) の 10 進展開、ガンマはオイラーのガンマ関数)」。 整数数列のオンライン百科事典 。OEIS Foundation。
^ Sloane, N. J. A. (編)。「数列 A245885 (ガンマ(7/2) の 10 進展開、ガンマはオイラーのガンマ関数)」。 整数数列のオンライン百科事典 。OEIS 財団。
^ 「対数ガンマ関数」。Wolfram MathWorld 。 2019年 1月3日 閲覧 。
^ 「レオンハルト・オイラーの積分:ガンマ関数の歴史的プロファイル」 (PDF) 。 2014年9月12日時点のオリジナルより アーカイブ (PDF) 。 2022年 4月11日 閲覧。
^ Blagouchine , Iaroslav V. (2015). 「有理数引数における最初の一般化スティルチェス定数の閉形式の評価と関連するいくつかの合計に関する定理」。Journal of Number Theory。148 : 537–592。arXiv : 1401.3724。doi : 10.1016 / j.jnt.2014.08.009 。
^ Alexejewsky、WP (1894)。 「Über eine Classe von Funktionen, die der Gammafunktion analog sind」[ガンマ関数に類似した関数のクラスについて]。 ライプツィヒ ワイドマンシェ ブッフハンドルング 。 46 : 268–275。
^ Barnes, EW (1899). 「 G 関数の理論 」. Quart. J. Math . 31 : 264–314.
^ Adamchik, Victor S. (1998). 「負の順序のポリガンマ関数」 J. Comput. Appl. Math . 100 (2): 191–199. doi : 10.1016/S0377-0427(98)00192-7 .
^ Gosper, RW (1997). 「 特殊関数における q 級数と関連トピック」 J. Am. Math. Soc . 14 .
∫
n
/
4
m
/
6
log
F
(
z
)
d
z
{\displaystyle \textstyle \int _{n/4}^{m/6}\log F(z)\,dz}
^ Espinosa, Olivier; Moll, Victor H. (2002). 「Hurwitzゼータ関数を含むいくつかの積分について: パート1」. The Ramanujan Journal . 6 (2): 159–188. doi :10.1023/A:1015706300169. S2CID 128246166.
^ Bailey, David H.; Borwein, David; Borwein, Jonathan M. (2015). 「オイラー対数ガンマ積分とTornheim-Wittenゼータ関数について」. The Ramanujan Journal . 36 (1–2): 43–68. doi :10.1007/s11139-012-9427-1. S2CID 7335291.
^ Amdeberhan, T.; Coffey, Mark W.; Espinosa, Olivier; Koutschan, Christoph; Manna, Dante V.; Moll, Victor H. (2011). 「loggamma のべき乗の積分」. Proc. Amer. Math. Soc . 139 (2): 535–545. doi : 10.1090/S0002-9939-2010-10589-0 .
^ EA Karatsuba、「超越関数の高速評価」。Probl. Inf. Transm. Vol.27、No.4、pp. 339–360 (1991)。
^ EA Karatsuba、「超越関数の高速評価のための新しい方法について」、Russ. Math. Surv. Vol.46、No.2、pp. 246–247 (1991)。
^ EA Karatsuba「高速アルゴリズムと FEE メソッド」。
^ Borwein, JM; Zucker, IJ (1992). 「第一種完全楕円積分を用いた小さな有理数に対するガンマ関数の高速評価」 IMA Journal of Numerical Analysis . 12 (4): 519–526. doi :10.1093/IMANUM/12.4.519.
^ Werner, Helmut; Collinge, Robert (1961). 「ガンマ関数のチェビシェフ近似」. Math. Comput . 15 (74): 195–197. doi :10.1090/S0025-5718-61-99220-1. JSTOR 2004230.
^ Michon, GP 「三角法と基本関数」Wayback Machine で 2010 年 1 月 9 日にアーカイブ 。Numericana。2007 年 5 月 5 日閲覧。
^ Chaudry, MA; Zubair, SM (2001). 不完全ガンマ関数のクラスとその応用について。 ボカラトン : CRC プレス。p. 37。ISBN 1-58488-143-7 。
^ ライス、JA(1995)。 数理統計とデータ分析 ( 第2版)。ベルモント:ダックスベリープレス。pp.52-53。ISBN 0-534-20934-3 。
^ Borwein, J.; Bailey, DH & Girgensohn, R. (2003). 数学における実験 . AK Peters. p. 133. ISBN 978-1-56881-136-9 。
^ 「自然階乗 n! の補間、または実階乗関数の誕生 (1729 - 1826)」。
^ オイラーの論文は、 『Commentarii academiae scientiarum Petropolitanae 5』、1738 年、36–57 に掲載されました。 The Euler Archive の「E19 -- De progressionibus transcentibus seu quarum termini generales algebraice dari nequeunt」を参照してください。これには、元の記事のスキャン コピーが含まれています。
^ Remmert, R. (2006). 複素関数理論における古典的な話題 。 翻訳:Kay, LD Springer。ISBN 978-0-387-98221-2 。
^ Lanczos, C. (1964). 「ガンマ関数の高精度近似」. Journal of the Society for Industrial and Applied Mathematics, Series B: Numerical Analysis . 1 (1): 86. Bibcode :1964SJNA....1...86L. doi :10.1137/0701008.
^ イルケル・イナム、エンギン・ビュユカシュク(2019)。国際計算数論秋季学校のノート。シュプリンガー。p. 205。ISBN 978-3-030-12558-5 。 205ページの抜粋
^ Knuth, DE (1997). コンピュータプログラミングの芸術 . 第 1 巻 (基本アルゴリズム). Addison-Wesley. ISBN 0-201-89683-4 。
^ Borwein, Jonathan M. ; Corless, Robert M. (2017). 「月刊誌におけるガンマと階乗」. American Mathematical Monthly . 125 (5). Mathematical Association of America: 400–24. arXiv : 1703.05349 . Bibcode :2017arXiv170305349B. doi :10.1080/00029890.2018.1420983. S2CID 119324101.
^ 「Gamma_functionの歴史は何ですか?」 yearis.com . 2022年 11月5日 閲覧 。
^ Berry, M. (2001 年 4 月). 「特殊関数はなぜ特別なのか?」 Physics Today .
^ 「microsoft/calculator」 。GitHub 。 2020年 12月25日 閲覧 。
^ "gnome-calculator". GNOME.org . 2023年 3月3日 閲覧 。
この記事には、 Creative Commons Attribution-ShareAlike 3.0 Unported Licenseに基づいてライセンスされているが、 GFDL に基づいてライセンスされていない Citizendium の 記事「ガンマ関数」の資料が組み込まれています 。
さらに読む
アブラモウィッツ、ミルトン、ステグン、アイリーン A. 編 (1972)。「第 6 章」。 数式、グラフ、および数表付き数学関数ハンドブック 。ニューヨーク: ドーバー。
Andrews, GE ; Askey, R.; Roy, R. (1999) 「第 1 章 (ガンマ関数とベータ関数)」。 特殊関数 。ニューヨーク: Cambridge University Press。ISBN 978-0-521-78988-2 。
アルティン、エミル (2006)。「ガンマ関数」。ローゼン、マイケル (編)。 エミル・アルティンの解説: 選集 。数学史。第 30 巻。プロビデンス、ロードアイランド州: アメリカ数学協会。
Askey, R. ; Roy, R. (2010)、「ガンマ関数」、 Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (編)、 NIST Handbook of Mathematical Functions 、Cambridge University Press、 ISBN 978-0-521-19225-5 、 MR 2723248 。
バーコフ、ジョージ D. (1913)。「ガンマ関数に関する注記」。Bull . Amer. Math. Soc . 20 (1): 1–10. doi : 10.1090/s0002-9904-1913-02429-7 . MR 1559418。
PE ベーマー (1939 年)。 Differenzengleichungen und bestimmte Integrale [ 微分方程式と定積分 ]。ライプツィヒ:ケーラーフェルラーク。
デイビス、フィリップ J. (1959)。「レオンハルト オイラーの積分: ガンマ関数の歴史的概要」。 アメリカ数学月刊誌 。66 ( 10): 849–869。doi :10.2307/2309786。JSTOR 2309786。
ポスト、エミール(1919)。「一般化されたガンマ関数」。 数学年報 。第2シリーズ。20 ( 3):202–217。doi :10.2307/1967871。JSTOR 1967871。 2021年3 月 5 日 閲覧 。
Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007)。「セクション 6.1. ガンマ関数」。 数値レシピ: 科学計算の技術 (第 3 版)。ニューヨーク: Cambridge University Press。ISBN 978-0-521-88068-8 。
オレゴン州ロックテッシェル (1922 年)。 複雑な引数のガンマ関数の計算方法 [ 複素引数のガンマ関数を計算する方法 ]。ドレスデン: ドレスデン工科大学 。
テム、ニコ・M. (1996)。 特殊関数:数理物理学の古典関数入門 。ニューヨーク:ジョン・ワイリー・アンド・サンズ 。ISBN 978-0-471-11313-3 。
ウィテカー、ET ; ワトソン、GN (1927)。 現代分析の講座 。ケンブリッジ大学出版局 。ISBN 978-0-521-58807-2 。
外部リンク
ウィキメディア・コモンズには、ガンマと関連機能 に関連するメディアがあります 。
NIST 数学関数デジタルライブラリ:ガンマ関数
Pascal Sebah と Xavier Gourdon。 ガンマ関数入門 。PostScript および HTML 形式。
std::tgamma の C++ リファレンス
ガンマ関数に関する問題の例は、Exampleproblems.com で見つかります。
「ガンマ関数」、 数学百科事典 、 EMS Press 、2001 [1994]
Wolfram ガンマ関数評価器 (任意精度)
「ガンマ 」。Wolfram 関数サイト。
MathPages の n 球の体積とガンマ関数