logo

オイラーのガンマ関数に対する極限公式の導出 📂関数

オイラーのガンマ関数に対する極限公式の導出

1

ガンマ関数について Γ:(0,)R\Gamma : (0, \infty) \to \mathbb{R} が成立する。 Γ(x)=limnnxn!x(x+1)(x+2)(x+n) \Gamma (x) = \lim_{n \to \infty} {{n^x n!} \over {x(x+1)(x+2) \cdots (x+n) }}

説明

以前知っていたガンマ関数は積分形 Γ(x)=0tx1etdt \Gamma (x) = \int_{0}^{\infty} t^{x-1} e^{-t} dt で全く異なる形をしていたが、1729年にオイラーが二つの表現が完全に同じであることを証明した。この記事で紹介する導出は元のものより少し簡略化されているが、理解する上で本質的な問題はないだろう。

導出

Γn(x):=0ntx1(1tn)ndt\displaystyle \Gamma_{n}(x) := \int_{0}^{n} t^{x-1} \left( 1 - { t \over n } \right) ^{n} dtet=limn(1tn)n\displaystyle e^{-t} = \lim_{n \to \infty } \left( 1 - { t \over n } \right) ^{-n}とすると limnΓn(x)=limn0ntx1(1tn)ndt=0tx1etdt=Γ(x) \begin{align*} \lim_{n \to \infty} \Gamma_{n}(x) =& \lim _{n \to \infty} \int_{0}^{n} t^{x-1} \left( 1 - { t \over n } \right) ^{n} dt \\ =& \int_{0}^{\infty} t^{x-1} e^{-t} dt \\ =& \Gamma (x) \end{align*} となる。一方で、Γn(x)=0ntx1(1tn)ndt\displaystyle \Gamma_{n}(x) = \int_{0}^{n} t^{x-1} \left( 1 - { t \over n } \right) ^{n} dt から u=tnu = {t \over n}への置換を行うと Γn(x)=0ntx1(1tn)ndt=01(nu)x1(1u)nndu=nx01ux1(1u)ndu \begin{align*} \Gamma_{n}(x) =& \int_{0}^{n} t^{x-1} \left( 1 - { t \over n } \right) ^{n} dt \\ =& \int_{0}^{1} (nu) ^{x-1} ( 1 - u ) ^{n} n du \\ =& n^{x} \int_{0}^{1} u^{x-1} ( 1 - u ) ^{n} du \end{align*} となる。部分積分法により 01ux1(1u)ndu=[(1x)ux(1u)n]0101(1x)uxn(1u)n1du=(nx)01ux(1u)n1du=(nx)(n1x+1)01ux+1(1u)n2du=(nx)(2x+n2)01ux+n2(1u)1du=(nx)(2x+n2)(1x+n1)01ux+n1(1u)0du=(nx)(2x+n2)(1x+n1)01ux+n1du=(nx)(2x+n2)(1x+n1)[1x+nux+n]01=(nx)(2x+n2)(1x+n1)(1x+n)=n(n1)(n2)211x(x+1)(x+2)(x+n)=n!x(x+1)(x+2)(x+n) \begin{align*} & \int_{0}^{1} u^{x-1} ( 1 - u ) ^{n} du \\ =& \left[ \left( { 1 \over x} \right) u^x (1-u)^n \right] _{0}^{1} - \int_{0}^{1} - \left( { 1 \over x} \right) u^x n ( 1 - u ) ^{n-1} du \\ =& \left( { n \over x} \right) \int_{0}^{1} u^x ( 1 - u ) ^{n-1} du \\ =& \left( { n \over x} \right) \left( { {n-1} \over {x+1}} \right) \int_{0}^{1} u^{x+1} ( 1 - u ) ^{n-2} du \\ \vdots& \\ =& \left( { n \over x} \right) \cdots \left( { 2 \over {x+n-2}} \right) \int_{0}^{1} u^{x+n-2} ( 1 - u ) ^{1} du \\ =& \left( { n \over x} \right) \cdots \left( { 2 \over {x+n-2}} \right) \left( { 1 \over {x+n-1}} \right) \int_{0}^{1} u^{x+n-1} ( 1 - u ) ^{0} du \\ =& \left( { n \over x} \right) \cdots \left( { 2 \over {x+n-2}} \right) \left( { 1 \over {x+n-1}} \right) \int_{0}^{1} u^{x+n-1} du \\ =& \left( { n \over x} \right) \cdots \left( { 2 \over {x+n-2}} \right) \left( { 1 \over {x+n-1}} \right) \left[ {{ 1 } \over { x+n }} u^{x+n} \right]_{0}^{1} \\ =& \left( { n \over x} \right) \cdots \left( { 2 \over {x+n-2}} \right) \left( { 1 \over {x+n-1}} \right) \left( { 1 \over {x+n}} \right) \\ =& { {n(n-1)(n-2) \cdots 2 \cdot 1 \cdot 1 } \over {x(x+1)(x+2) \cdots (x + n ) } } \\ =& { { n! } \over {x(x+1)(x+2) \cdots (x + n ) } } \end{align*} を得て、したがって Γn(x)=nxn!x(x+1)(x+2)(x+n) \Gamma_{n}(x) = n^{x} { { n! } \over {x(x+1)(x+2) \cdots (x + n ) } } を得る。前に Γ(x)=limnΓn(x)\displaystyle \Gamma (x) = \lim_{n \to \infty} \Gamma_{n}(x) であることを示したので、次が成立する。 Γ(x)=limnΓn(x)=limnnxn!x(x+1)(x+2)(x+n) \begin{align*} \Gamma (x) =& \lim_{n \to \infty} \Gamma_{n}(x) \\ =& \lim_{n \to \infty} { { n^{x} n! } \over {x(x+1)(x+2) \cdots (x + n ) } } \end{align*}

参考