ここで行き詰まった? AI で解決しよう
Calculus · real student question

0 ≤ x ≤ k、0 ≤ y ≤ k − x、0 ≤ z ≤ k − x − y で定まる領域上で、xyz に e の −(x + y + z) 乗を掛けたものの三重積分を求めよ。

Question

次を求めよ。

I=∫0k ⁣ ⁣∫0k−x ⁣ ⁣∫0k−x−yxyz e−x−y−z dz dy dxI = \int_0^k\!\!\int_0^{k-x}\!\!\int_0^{k-x-y} xyz\,e^{-x-y-z}\,dz\,dy\,dx

Step-by-step solution

  1. 領域が四面体であることを見抜く。 入れ子の積分限界はちょうど次の条件

    x≥0,y≥0,z≥0,x+y+z≤kx \ge 0, \quad y \ge 0, \quad z \ge 0, \quad x + y + z \le k

    を表しており、原点と (k,0,0)(k,0,0)、(0,k,0)(0,k,0)、(0,0,k)(0,0,k) を頂点とする単体である。こう見ることで、x+y+zx+y+z にしか依存しない指数関数を等位面に沿ってくくり出せるようになる。

  2. 指数関数が各断面上で一定であることを使う。 u=x+y+zu = x + y + z とおく。x+y+z=ux+y+z = u という断面上では因子 e−ue^{-u} は一定なので、各断面上の xyzxyz の「量」が分かれば、積分は uu についての 1 つの積分に帰着する。

  3. 部分単体上での xyz の重みを計算する。 スケーリング則により

    ∫x+y+z≤uxyz dV=u6720\int_{x+y+z\le u} xyz\,dV = \frac{u^6}{720}

    となる(3 次元領域上の 3 次同次な被積分関数は u3+3=u6u^{3+3} = u^6 のようにスケールし、定数はディリクレ積分から 1720=1! 1! 1!6!\frac{1}{720} = \frac{1!\,1!\,1!}{6!} となる)。uu で微分すれば断面上の密度 u5120\frac{u^5}{120} が得られる。

  4. 1 次元の積分に帰着させる。

    I=1120∫0ku5e−u duI = \frac{1}{120}\int_0^k u^5 e^{-u}\,du

    これは不完全ガンマ積分であり、部分積分を繰り返せば(あるいは標準的な原始関数から)

    ∫u5e−u du=−e−u(u5+5u4+20u3+60u2+120u+120)\int u^5 e^{-u}\,du = -e^{-u}\left(u^5 + 5u^4 + 20u^3 + 60u^2 + 120u + 120\right)

    が得られる。

  5. 0 から k まで評価する。

    ∫0ku5e−u du=120−e−k(k5+5k4+20k3+60k2+120k+120)\int_0^k u^5 e^{-u}\,du = 120 - e^{-k}\left(k^5 + 5k^4 + 20k^3 + 60k^2 + 120k + 120\right)

    120120 で割ると

    I=1−e−k(1+k+k22+k36+k424+k5120)I = 1 - e^{-k}\left(1 + k + \frac{k^2}{2} + \frac{k^3}{6} + \frac{k^4}{24} + \frac{k^5}{120}\right)

    括弧の中は eke^{k} のマクローリン級数の最初の 6 項である。これは正規化不完全ガンマ関数 P(6,k)P(6, k) であり、平均 kk のポアソン確率変数が 66 以上となる確率にほかならない。

  6. 数値と極限で確認する。 k→∞k \to \infty のとき I→1I \to 1、k→0k \to 0 のとき I→0I \to 0 となる。四面体上で入れ子の数値求積を直接行うと、k=1k = 1 で 0.00059418480.0005941848、k=2.5k = 2.5 で 0.04202103820.0420210382、k=5k = 5 で 0.38403934520.3840393452 が得られ、いずれも閉じた式と小数第 10 位まで一致する。

Answer

1−e−k(1+k+k22+k36+k424+k5120)1 - e^{-k}\left(1 + k + \frac{k^2}{2} + \frac{k^3}{6} + \frac{k^4}{24} + \frac{k^5}{120}\right)

Need to solve a different problem like this? Open the solver →