ガンマ関数・ベータ関数の公式|Γ(n+1)=n!・B(p,q)=Γ(p)Γ(q)/Γ(p+q)の証明と例題【大学院入試・微分積分】

ガンマ関数・ベータ関数の公式|Γ(n+1)=n!・B(p,q)=Γ(p)Γ(q)/Γ(p+q)の証明と例題【大学院入試・微分積分】

合格大作戦 算数・数学公式辞典 › 関数

大学院入試

Γ(s)=∫0∞xs−1e−x dx,B(p,q)=Γ(p)Γ(q)Γ(p+q)\Gamma(s)=\int_0^\infty x^{s-1}e^{-x}\,dx,\qquad B(p,q)=\frac{\Gamma(p)\Gamma(q)}{\Gamma(p+q)}

Γ(n+1)=n!\Gamma(n+1)=n!、Γ(12)=π\Gamma(\frac12)=\sqrt\pi、B(p,q)=∫01xp−1(1−x)q−1dxB(p,q)=\displaystyle\int_0^1x^{p-1}(1-x)^{q-1}dx

この公式のポイント

  • Γ(s+1)=sΓ(s)\Gamma(s+1)=s\Gamma(s)(部分積分)で、Γ(n+1)=n!\Gamma(n+1)=n!
  • Γ(12)=π\Gamma(\frac12)=\sqrt\pi はガウス積分から
  • B(p,q)=Γ(p)Γ(q)Γ(p+q)B(p,q)=\dfrac{\Gamma(p)\Gamma(q)}{\Gamma(p+q)} はヤコビアン uu の変数変換で示す
  • 指数は 11 ずれるので注意
うかるくん

うかるくん

ガンマ関数って、階乗とどういう関係なの?

合格先生

合格先生

部分積分で Γ(s+1)=sΓ(s)\Gamma(s+1)=s\Gamma(s) が出る。Γ(1)=1\Gamma(1)=1 からくり返すと Γ(n+1)=n!\Gamma(n+1)=n!。階乗をなめらかにつないだ関数だ。

ガンマ関数・ベータ関数の公式

ガンマ関数は、階乗 n!n! を整数でない数にまで広げた関数です。ベータ関数は ∫01xa(1−x)b dx\displaystyle\int_0^1x^a(1-x)^b\,dx の形の積分で、ガンマ関数で表せます。

ガンマ関数とベータ関数

Γ(s)=∫0∞xs−1e−x dx(s>0)\Gamma(s)=\int_0^\infty x^{s-1}e^{-x}\,dx\quad(s>0)Γ(s+1)=s Γ(s),Γ(n+1)=n!,Γ ⁣(12)=π\Gamma(s+1)=s\,\Gamma(s),\qquad\Gamma(n+1)=n!,\qquad\Gamma\!\left(\tfrac12\right)=\sqrt\piB(p,q)=∫01xp−1(1−x)q−1 dx=Γ(p)Γ(q)Γ(p+q)(p,q>0)B(p,q)=\int_0^1x^{p-1}(1-x)^{q-1}\,dx=\frac{\Gamma(p)\Gamma(q)}{\Gamma(p+q)}\quad(p,q>0)

三角関数の形:B(p,q)=2∫0π2sin⁡2p−1θcos⁡2q−1θ dθ\displaystyle B(p,q)=2\int_0^{\frac\pi2}\sin^{2p-1}\theta\cos^{2q-1}\theta\,d\theta

xyO12341260!=11!=12!=23!=6Γ(1/2)=√π≒1.77y=Γ(x)
y=Γ(x)y=\Gamma(x) は、整数のところで階乗の値を通る。Γ(12)=π\Gamma(\frac12)=\sqrt\pi

証明

① 収束(s>0s>0 のとき)

0<x≤10<x\le1 では 0<xs−1e−x≤xs−10<x^{s-1}e^{-x}\le x^{s-1} で、∫01xs−1dx=1s\displaystyle\int_0^1x^{s-1}dx=\frac1s は s>0s>0 なら有限。x≥1x\ge1 では、xs−1e−x2x^{s-1}e^{-\frac x2} が有界(指数関数は多項式より速く小さくなる)なので xs−1e−x≤Ce−x2x^{s-1}e^{-x}\le Ce^{-\frac x2} となり有限です。

② Γ(s+1)=sΓ(s)\Gamma(s+1)=s\Gamma(s)(部分積分)

Γ(s+1)=∫0∞xse−x dx=[−xse−x]0∞+s∫0∞xs−1e−x dx=s Γ(s)\Gamma(s+1)=\int_0^\infty x^se^{-x}\,dx=\Big[-x^se^{-x}\Big]_0^\infty+s\int_0^\infty x^{s-1}e^{-x}\,dx=s\,\Gamma(s)

s>0s>0 なので x=0x=0 でも x→∞x\to\infty でも xse−x→0x^se^{-x}\to0。Γ(1)=∫0∞e−xdx=1\Gamma(1)=\displaystyle\int_0^\infty e^{-x}dx=1 から、Γ(n+1)=n⋅(n−1)⋯1⋅Γ(1)=n!\Gamma(n+1)=n\cdot(n-1)\cdots1\cdot\Gamma(1)=n! です。

③ Γ(12)=π\Gamma(\frac12)=\sqrt\pi

x=t2x=t^2(dx=2t dtdx=2t\,dt)とおくと Γ(12)=∫0∞e−t2t 2t dt=2∫0∞e−t2dt=π\Gamma(\frac12)=\displaystyle\int_0^\infty\frac{e^{-t^2}}{t}\,2t\,dt=2\int_0^\infty e^{-t^2}dt=\sqrt\pi(ガウス積分)。

④ B(p,q)=Γ(p)Γ(q)Γ(p+q)B(p,q)=\dfrac{\Gamma(p)\Gamma(q)}{\Gamma(p+q)}

Γ(p)Γ(q)=∫0∞ ⁣ ⁣∫0∞xp−1yq−1e−(x+y) dx dy\Gamma(p)\Gamma(q)=\displaystyle\int_0^\infty\!\!\int_0^\infty x^{p-1}y^{q-1}e^{-(x+y)}\,dx\,dy で、x=uv, y=u(1−v)x=uv,\ y=u(1-v)(u>0, 0<v<1u>0,\ 0<v<1)と変数変換します。

∂(x,y)∂(u,v)=det⁡(vu1−v−u)=−uv−u(1−v)=−u\frac{\partial(x,y)}{\partial(u,v)}=\det\begin{pmatrix}v&u\\1-v&-u\end{pmatrix}=-uv-u(1-v)=-u

∣J∣=u|J|=u なので

Γ(p)Γ(q)=∫0∞up+q−1e−u du∫01vp−1(1−v)q−1 dv=Γ(p+q) B(p,q)\Gamma(p)\Gamma(q)=\int_0^\infty u^{p+q-1}e^{-u}\,du\int_0^1v^{p-1}(1-v)^{q-1}\,dv=\Gamma(p+q)\,B(p,q)

(被積分関数が正なので、広義の重積分を累次積分に分けてよいことを使っています。)三角関数の形は、x=sin⁡2θx=\sin^2\theta(dx=2sin⁡θcos⁡θ dθdx=2\sin\theta\cos\theta\,d\theta)と置換すると出ます。

xyO2468100.51x e^(−x):面積 1x² e^(−x):面積 2x³ e^(−x):面積 6山の頂点は x=s−1面積 Γ(s)=(s−1)! は s とともに急に大きくなる
xs−1e−xx^{s-1}e^{-x} の下の面積が Γ(s)=(s−1)!\Gamma(s)=(s-1)!

大学院入試では

院試では、見た目のちがう積分を「ガンマ関数・ベータ関数の形」に置換して、一気に値を出す問題がよく出ます。

大学院入試

例題1(ガンマ関数の形にする)

∫0∞x3e−2x dx\displaystyle\int_0^\infty x^3e^{-2x}\,dx を求めなさい。

解き方

t=2xt=2x とおくと x=t2, dx=dt2x=\dfrac t2,\ dx=\dfrac{dt}2。

∫0∞t38e−tdt2=116Γ(4)=3!16=38\displaystyle\int_0^\infty\frac{t^3}{8}e^{-t}\frac{dt}2=\frac1{16}\Gamma(4)=\frac{3!}{16}=\frac38。

一般に ∫0∞xne−axdx=n!an+1\displaystyle\int_0^\infty x^ne^{-ax}dx=\frac{n!}{a^{n+1}} です。

答え 38\dfrac38

大学院入試

例題2(半整数での値)

Γ ⁣(52)\Gamma\!\left(\dfrac52\right) を求めなさい。

解き方

Γ(s+1)=sΓ(s)\Gamma(s+1)=s\Gamma(s) を2回使って Γ ⁣(52)=32Γ ⁣(32)=32⋅12Γ ⁣(12)=34π\Gamma\!\left(\dfrac52\right)=\dfrac32\Gamma\!\left(\dfrac32\right)=\dfrac32\cdot\dfrac12\Gamma\!\left(\dfrac12\right)=\dfrac34\sqrt\pi。

答え 3π4\dfrac{3\sqrt\pi}{4}

大学院入試

例題3(ベータ関数)

∫01x2(1−x)3 dx\displaystyle\int_0^1x^2(1-x)^3\,dx と ∫0π2sin⁡4θcos⁡2θ dθ\displaystyle\int_0^{\frac\pi2}\sin^4\theta\cos^2\theta\,d\theta を求めなさい。

解き方

1つめは B(3,4)=Γ(3)Γ(4)Γ(7)=2! 3!6!=12720=160B(3,4)=\dfrac{\Gamma(3)\Gamma(4)}{\Gamma(7)}=\dfrac{2!\,3!}{6!}=\dfrac{12}{720}=\dfrac1{60}。

2つめは 2p−1=4, 2q−1=22p-1=4,\ 2q-1=2 より p=52, q=32p=\dfrac52,\ q=\dfrac32。12B ⁣(52,32)=12⋅Γ(52)Γ(32)Γ(4)=12⋅34π⋅12π6=π32\dfrac12B\!\left(\dfrac52,\dfrac32\right)=\dfrac12\cdot\dfrac{\Gamma(\frac52)\Gamma(\frac32)}{\Gamma(4)}=\dfrac12\cdot\dfrac{\frac34\sqrt\pi\cdot\frac12\sqrt\pi}{6}=\dfrac{\pi}{32}。

答え 160\dfrac1{60}、π32\dfrac\pi{32}

学部の授業とのつながり・院試での注意

  • どの授業で習うか:大学1〜2年の微分積分学(広義積分)の応用として習います。確率統計(ガンマ分布・ベータ分布・カイ2乗分布)や物理(球の体積の一般化)でも使います。
  • 答案で示すこと:①広義積分が収束する範囲(s>0s>0、p,q>0p,q>0)、②どの置換でガンマ・ベータの形にしたか、③Γ(s+1)=sΓ(s)\Gamma(s+1)=s\Gamma(s) の部分積分では、はしの項が 00 になる理由。
  • よく出る形:∫0∞xne−axdx\displaystyle\int_0^\infty x^ne^{-ax}dx、∫0∞x2ne−x2dx\displaystyle\int_0^\infty x^{2n}e^{-x^2}dx(t=x2t=x^2)、∫0π2sin⁡mθcos⁡nθ dθ\displaystyle\int_0^{\frac\pi2}\sin^m\theta\cos^n\theta\,d\theta(ウォリス積分の一般化)、BB と Γ\Gamma の関係の証明そのもの。
  • スターリングの公式は、Γ(n+1)=n!\Gamma(n+1)=n! の大きな nn での近似です。

よくある間違い

腕でバツを作る合格先生

合格先生

Γ(n)=n!\Gamma(n)=n! ではない。正しくは Γ(n)=(n−1)!\Gamma(n)=(n-1)! だ。1つずれるから気をつけよう。

  • Γ(n)=n!\Gamma(n)=n! と書く:正しくは Γ(n)=(n−1)!\Gamma(n)=(n-1)!、Γ(n+1)=n!\Gamma(n+1)=n!。1つずれます。
  • ベータ関数の指数をそのまま入れる:∫01x2(1−x)3dx\displaystyle\int_0^1x^2(1-x)^3dx は B(2,3)B(2,3) ではなく B(3,4)B(3,4)。指数に 11 を足したものが p, qp,\ q です。
  • 三角関数の形で係数2を忘れる:∫0π2sin⁡2p−1θcos⁡2q−1θ dθ=12B(p,q)\displaystyle\int_0^{\frac\pi2}\sin^{2p-1}\theta\cos^{2q-1}\theta\,d\theta=\frac12B(p,q) です。

練習問題

ひらめいたサクラちゃん

サクラちゃん

指数に1を足して、p, qp,\ q を決めればいいんだね!

  1. 大学院入試 Γ ⁣(72)\Gamma\!\left(\dfrac72\right) を求めなさい。
    答えと解説を見る52⋅32⋅12π\dfrac52\cdot\dfrac32\cdot\dfrac12\sqrt\pi。答え 15π8\dfrac{15\sqrt\pi}8
  2. 大学院入試 ∫01x3(1−x)4 dx\displaystyle\int_0^1x^3(1-x)^4\,dx を求めなさい。
    答えと解説を見るB(4,5)=3! 4!8!=14440320B(4,5)=\dfrac{3!\,4!}{8!}=\dfrac{144}{40320}。答え 1280\dfrac1{280}
  3. 大学院入試 ∫0∞x4e−x2 dx\displaystyle\int_0^\infty x^4e^{-x^2}\,dx を、t=x2t=x^2 と置換してガンマ関数で求めなさい。
    答えと解説を見るx=t, dx=dt2tx=\sqrt t,\ dx=\dfrac{dt}{2\sqrt t} より ∫0∞t2e−tdt2t=12Γ ⁣(52)=12⋅3π4\displaystyle\int_0^\infty t^2e^{-t}\frac{dt}{2\sqrt t}=\frac12\Gamma\!\left(\frac52\right)=\frac12\cdot\frac{3\sqrt\pi}4。答え 3π8\dfrac{3\sqrt\pi}8

関連する公式

よくある質問

ガンマ関数とは?

Γ(s)=∫_0^∞ x^(s−1) e^(−x) dx(s>0)で定める関数です。Γ(n+1)=n! となるので、階乗を整数でない数まで広げたものと考えられます。

Γ(1/2)=√π になるのはなぜ?

x=t² と置換すると 2∫_0^∞ e^(−t²) dt になり、これはガウス積分で √π です。

ベータ関数とガンマ関数の関係は?

B(p,q)=Γ(p)Γ(q)/Γ(p+q) です。Γ(p)Γ(q) を重積分にして x=uv, y=u(1−v) と変数変換すると示せます。

大学院入試ではどう使いますか?

∫x^n e^(−ax)dx や ∫sin^m θ cos^n θ dθ などを、置換してガンマ関数・ベータ関数の形にし、値を一気に求めるのに使います。

あわせて読みたい

出典・参考

  • 大学初年級の微分積分学(広義積分)

最終更新:2026年10月11日/作成:大学受験合格大作戦