Lemma数学ロードマップ

01 微分積分学 · 第 9 章

重積分

目安 12〜16 時間定理など 12演習 8 問

この章の目標

  • 長方形上のリーマン積分を定義し、連続関数が可積分であることを示せる
  • フビニの定理を証明し、重積分を累次積分で計算できる
  • ジョルダン可測集合(面積・体積の定まる集合)を理解し、一般の領域上で積分できる
  • 変数変換公式の意味(ヤコビアンが体積の拡大率であること)を理解し、極座標・球座標で計算できる
  • 広義重積分を扱い、ガウス積分 ∫−∞∞e−x2dx=π\int_{-\infty}^{\infty}e^{-x^2}dx = \sqrt{\pi}、B=ΓΓ/ΓB = \Gamma\Gamma/\Gamma、nn 次元球の体積を導ける

前提:第5章 リーマン積分、第7章、02-linear-algebra 第4章 行列式

9.1 長方形上の積分

高校では立体の体積を「断面積を積分する」ことで求めた。これは後で述べるフビニの定理の特別な場合である。大学では、まず多変数関数の積分を 1 変数と同じくダルブー和で定義し、それを 1 変数の積分の繰り返し(累次積分)で計算できることを定理として証明する。

Rn\mathbb{R}^n の閉長方形(閉直方体)R=[a1,b1]×⋯×[an,bn]R = [a_1, b_1] \times \cdots \times [a_n, b_n] の体積を v(R)=∏i(bi−ai)v(R) = \prod_i(b_i - a_i) とする。各辺 [ai,bi][a_i, b_i] の分割 PiP_i の組 P=(P1,…,Pn)P = (P_1, \dots, P_n) を RR の分割といい、PP によって RR は有限個の小長方形 SS に分けられる。有界関数 f ⁣:R→Rf\colon R \to \mathbb{R} に対し、MS=sup⁡SfM_S = \sup_S f、mS=inf⁡Sfm_S = \inf_S f として

U(f,P)=∑SMSv(S),L(f,P)=∑SmSv(S)U(f, P) = \sum_S M_S v(S), \qquad L(f, P) = \sum_S m_S v(S)

と定め、第5章と全く同様に上積分・下積分・可積分性・積分 ∫Rf=∫Rf(x) dx\int_R f = \int_R f(x)\ dx を定義する。補題 5.2(細分で上和は減り下和は増える)、定理 5.5(リーマンの判定条件:U−LU - L をいくらでも小さくできる)、定理 5.10(線形性・単調性・∣∫f∣≤∫∣f∣\left\lvert \int f \right\rvert \leq \int\lvert f \rvert)は証明も含めてそのまま成り立つ。

定理 9.1 閉長方形 RR 上の連続関数は可積分である。

証明. RR はコンパクトなので ff は一様連続である(定理 7.7)。ε>0\varepsilon > 0 に対し、∥x−y∥<δ⇒∣f(x)−f(y)∣<ε\lVert x - y \rVert < \delta \Rightarrow \lvert f(x) - f(y) \rvert < \varepsilon となる δ\delta をとり、すべての小長方形の直径が δ\delta 未満となる分割 PP をとれば、MS−mS≤εM_S - m_S \leq \varepsilon なので U(f,P)−L(f,P)≤εv(R)U(f,P) - L(f,P) \leq \varepsilon v(R)。□\square

9.2 フビニの定理

重積分を実際に計算するための基本定理である。Rp+q\mathbb{R}^{p+q} の長方形を R=A×BR = A \times B(A⊂RpA \subset \mathbb{R}^p、B⊂RqB \subset \mathbb{R}^q は長方形)と書き、点を (x,y)(x, y) と書く。

定理 9.2(フビニの定理, Fubini's theorem)ff が R=A×BR = A \times B で可積分であり、各 x∈Ax \in A について y↦f(x,y)y \mapsto f(x, y) が BB で可積分であるとする。このとき F(x)=∫Bf(x,y) dyF(x) = \int_B f(x, y)\ dy は AA で可積分で

∫A×Bf(x,y) d(x,y)=∫A(∫Bf(x,y) dy)dx\int_{A \times B}f(x, y)\,d(x,y) = \int_A\left(\int_B f(x, y)\,dy\right)dx

特に ff が RR 上連続なら、仮定はすべて満たされ、積分の順序を入れ替えた等式も成り立つ:

∫A(∫Bf(x,y) dy)dx=∫A×Bf=∫B(∫Af(x,y) dx)dy\int_A\left(\int_B f(x, y)\,dy\right)dx = \int_{A \times B}f = \int_B\left(\int_A f(x, y)\,dx\right)dy

証明. AA の分割 PAP_A と BB の分割 PBP_B から RR の分割 P=PA×PBP = P_A \times P_B を作ると、その小長方形は S×TS \times T(SS は PAP_A の、TT は PBP_B の小長方形)である。x∈Sx \in S を固定すると、区間についての加法性より

F(x)=∑T∫Tf(x,y) dy≥∑TmS×T(f) v(T)F(x) = \sum_T\int_T f(x, y)\,dy \geq \sum_T m_{S \times T}(f)\,v(T)

(TT 上で f(x,y)≥mS×T(f)f(x, y) \geq m_{S \times T}(f) だから)。x∈Sx \in S について下限をとり、v(S)v(S) を掛けて SS について足すと

L(F,PA)≥∑S,TmS×T(f) v(S)v(T)=L(f,P)L(F, P_A) \geq \sum_{S, T}m_{S \times T}(f)\,v(S)v(T) = L(f, P)

同様に U(F,PA)≤U(f,P)U(F, P_A) \leq U(f, P)。したがって

L(f,P)≤L(F,PA)≤U(F,PA)≤U(f,P)L(f, P) \leq L(F, P_A) \leq U(F, P_A) \leq U(f, P)

ff は可積分なので、U(f,P)−L(f,P)U(f, P) - L(f, P) はいくらでも小さくでき(任意の分割は PA×PBP_A \times P_B の形に細分できる)、FF は可積分で ∫AF=∫Rf\int_A F = \int_R f。ff が連続なら、ff も各 f(x,⋅)f(x, \cdot) も連続なので可積分であり、xx と yy の役割を入れ替えた議論も成り立つ。□\square

例 9.3 ∫[0,1]×[0,2](x+y2)=∫01(∫02(x+y2) dy)dx=∫01(2x+83)dx=1+83=113\int_{[0,1] \times [0,2]}(x + y^2) = \int_0^1\left(\int_0^2(x + y^2)\ dy\right)dx = \int_0^1\left(2x + \frac{8}{3}\right)dx = 1 + \frac{8}{3} = \frac{11}{3}。順序を入れ替えても ∫02(12+y2)dy=1+83\int_0^2\left(\frac{1}{2} + y^2\right)dy = 1 + \frac{8}{3} で一致する。

注意

積分可能性の仮定がないと、積分の順序は交換できないことがある。f(x,y)=x2−y2(x2+y2)2f(x, y) = \frac{x^2 - y^2}{(x^2 + y^2)^2} を (0,1]2(0, 1]^2 で考えると、∂∂yyx2+y2=f\frac{\partial}{\partial y}\frac{y}{x^2 + y^2} = f より ∫01f(x,y) dy=11+x2\int_0^1 f(x, y)\ dy = \frac{1}{1 + x^2} なので

∫01∫01f dy dx=π4,∫01∫01f dx dy=−π4\int_0^1\int_0^1 f\,dy\,dx = \frac{\pi}{4}, \qquad \int_0^1\int_0^1 f\,dx\,dy = -\frac{\pi}{4}

(後者は f(x,y)=−f(y,x)f(x, y) = -f(y, x) による)。ff は原点の近くで有界でなく、∣f∣\lvert f \rvert の広義積分は発散する。

9.3 ジョルダン可測集合

長方形以外の集合、例えば円板の上で積分するには、関数を 00 で拡張して長方形上で積分する。そのためには集合の「境界が薄い」ことが必要になる。

定義 9.4 (1) Z⊂RnZ \subset \mathbb{R}^n が体積 0(ジョルダン零集合, content zero)であるとは、任意の ε>0\varepsilon > 0 に対して、体積の総和が ε\varepsilon 未満の有限個の閉長方形で ZZ を覆えることをいう。 (2) 有界集合 D⊂RnD \subset \mathbb{R}^n と DD 上の有界関数 ff に対して、DD を含む閉長方形 RR をとり、ff を R∖DR \setminus D で 00 として拡張した関数 f~\tilde{f} が RR で可積分なとき、ff は DD で可積分であるといい、∫Df=∫Rf~\int_D f = \int_R\tilde{f} と定める(RR のとり方によらない)。 (3) 定数関数 11 が DD で可積分なとき、DD はジョルダン可測 (Jordan measurable) であるといい、vol⁡(D)=∫D1\operatorname{vol}(D) = \int_D 1 を DD の体積(n=2n = 2 なら面積)という。

高校の「面積」は、ここで初めて定義されたことになる。次の補題が基本的である。

補題 9.5 閉長方形 RR 上の有界関数 ff が、体積 0 の集合 ZZ の外で連続ならば、ff は RR で可積分である。

証明. ∣f∣≤M\lvert f \rvert \leq M、ε>0\varepsilon > 0 とする。ZZ を覆う閉長方形 R1,…,RkR_1, \dots, R_k で ∑v(Ri)<ε\sum v(R_i) < \varepsilon となるものをとり、各 RiR_i を少し広げて、ZZ が Ri′R_i' たちの内部の和に含まれ ∑v(Ri′)<2ε\sum v(R_i') < 2\varepsilon となる閉長方形 Ri′R_i' をとる。K=R∖⋃i(Ri′)∘K = R \setminus \bigcup_i(R_i')^\circ はコンパクトで、ff は KK の各点で(RR 上の関数として)連続である。このとき、ある δ>0\delta > 0 があって「x∈Kx \in K、y∈Ry \in R、∥x−y∥<δ\lVert x - y \rVert < \delta なら ∣f(x)−f(y)∣<ε\lvert f(x) - f(y) \rvert < \varepsilon」となる(そうでなければ ∥xj−yj∥→0\lVert x_j - y_j \rVert \to 0、∣f(xj)−f(yj)∣≥ε\lvert f(x_j) - f(y_j) \rvert \geq \varepsilon となる xj∈Kx_j \in K、yj∈Ry_j \in R があり、収束部分列の極限 x∈Kx \in K での連続性に反する)。すべての Ri′R_i' の各辺の端点を分点に含み、小長方形の直径が δ\delta 未満の分割 PP をとる。PP の小長方形 SS は、ある Ri′R_i' に含まれるか、どの (Ri′)∘(R_i')^\circ とも交わらない(このとき S⊂KS \subset K)。前者の体積の和は 2ε2\varepsilon 未満で MS−mS≤2MM_S - m_S \leq 2M、後者では MS−mS≤εM_S - m_S \leq \varepsilon。よって U(f,P)−L(f,P)≤4Mε+εv(R)U(f,P) - L(f,P) \leq 4M\varepsilon + \varepsilon v(R)。□\square

命題 9.6 有界集合 DD がジョルダン可測であるための必要十分条件は、境界 ∂D\partial D が体積 0 であることである。

証明. DD を含む長方形 RR を、D‾\overline{D} が RR の内部に含まれるようにとる。χ=1~\chi = \tilde{1}(DD の定義関数)は D∘D^\circ と R∖D‾R \setminus \overline{D} で局所的に定数なので連続であり、不連続点は ∂D\partial D に含まれる。よって ∂D\partial D が体積 0 なら補題 9.5 より可積分。逆に χ\chi が可積分とし、U(χ,P)−L(χ,P)<εU(\chi, P) - L(\chi, P) < \varepsilon となる分割 PP をとる。MS−mSM_S - m_S は SS が DD と DcD^c の両方と交わるとき 11、そうでないとき 00 なので、そのような SS の体積の和は ε\varepsilon 未満である。x∈∂Dx \in \partial D とする。x∈Dx \in D なら xx のどんな近傍も DcD^c と交わるので、xx を含む小長方形(その和集合は xx の近傍)のどれかは DcD^c の点を含み、それは x∈Dx \in D も含むので両方と交わる。x∉Dx \notin D でも同様。よって ∂D\partial D はそれらの小長方形で覆われ、体積 0 である。□\square

命題 9.7 (1) A⊂Rn−1A \subset \mathbb{R}^{n-1} を閉長方形、φ ⁣:A→R\varphi\colon A \to \mathbb{R} を連続とすると、グラフ {(x,φ(x))∣x∈A}\lbrace (x, \varphi(x)) \mid x \in A \rbrace は Rn\mathbb{R}^n で体積 0 である。 (2) DD がジョルダン可測で、ff が DD 上有界かつ連続ならば、ff は DD で可積分である。

証明. (1) φ\varphi は一様連続なので、ε>0\varepsilon > 0 に対して AA を小長方形 SS に分け、各 SS 上で φ\varphi の振幅(sup⁡−inf⁡\sup - \inf)を ε\varepsilon 未満にできる。グラフは S×[inf⁡Sφ,sup⁡Sφ]S \times [\inf_S\varphi, \sup_S\varphi] の和で覆われ、その体積の和は εv(A)\varepsilon v(A) 以下。(2) f~\tilde{f} は D∘D^\circ では ff に一致して連続、R∖D‾R \setminus \overline{D} では 00 で連続なので、不連続点は ∂D\partial D(体積 0)に含まれる。補題 9.5 による。□\square

例 9.8 (1) 円板、球、多角形など、有限個の連続関数のグラフ(と長方形の面)で囲まれた有界集合はジョルダン可測である。 (2) D=Q2∩[0,1]2D = \mathbb{Q}^2 \cap [0, 1]^2 はジョルダン可測でない(∂D=[0,1]2\partial D = [0,1]^2 は体積 0 でない。例 5.4 (2) と同じ理由)。可算集合であっても、ジョルダン可測とは限らない。ルベーグ測度(06-measure-integration 第2章)では、可算集合はすべて測度 0 になる。

9.4 累次積分

実際の計算で最もよく使うのは、次の形の領域である。

定理 9.9(縦線集合上の累次積分)φ1≤φ2\varphi_1 \leq \varphi_2 を [a,b][a, b] 上の連続関数とし、D={(x,y)∣a≤x≤b, φ1(x)≤y≤φ2(x)}D = \lbrace (x, y) \mid a \leq x \leq b,\ \varphi_1(x) \leq y \leq \varphi_2(x) \rbrace とする。DD 上の連続関数 ff について

∫Df=∫ab(∫φ1(x)φ2(x)f(x,y) dy)dx\int_D f = \int_a^b\left(\int_{\varphi_1(x)}^{\varphi_2(x)}f(x, y)\,dy\right)dx

Rn\mathbb{R}^n の同様の領域(x∈Ax \in A を長方形上で動かし、最後の変数が二つの連続関数の間を動く)についても同様である。

証明. ∂D\partial D は二つのグラフと二本の線分に含まれるので体積 0 であり、DD はジョルダン可測、ff は DD で可積分(命題 9.7)。D⊂R=[a,b]×[c,d]D \subset R = [a, b] \times [c, d] として f~\tilde{f} を考えると、各 xx で y↦f~(x,y)y \mapsto \tilde{f}(x, y) は [φ1(x),φ2(x)][\varphi_1(x), \varphi_2(x)] で f(x,y)f(x, y)、その外で 00 であり、不連続点は高々 2 個なので可積分で、積分は ∫φ1(x)φ2(x)f(x,y) dy\int_{\varphi_1(x)}^{\varphi_2(x)}f(x,y)\ dy。フビニの定理(定理 9.2)を適用すればよい。□\square

例 9.10(積分の順序交換)∫01(∫x1ey2dy)dx\int_0^1\left(\int_x^1 e^{y^2}dy\right)dx は、内側の積分が初等関数で書けないのでこのままでは計算できない。積分領域 {0≤x≤y≤1}\lbrace 0 \leq x \leq y \leq 1 \rbrace を「yy を固定して xx を動かす」形に見直すと

∫01(∫0yey2dx)dy=∫01yey2dy=e−12\int_0^1\left(\int_0^y e^{y^2}dx\right)dy = \int_0^1 ye^{y^2}dy = \frac{e - 1}{2}

例 9.11(カヴァリエリの原理)ジョルダン可測な立体 K⊂R3K \subset \mathbb{R}^3 の、高さ zz での断面 KzK_z の面積を s(z)s(z) とすると、フビニの定理より vol⁡(K)=∫s(z) dz\operatorname{vol}(K) = \int s(z)\ dz である(KzK_z が各 zz で可測で、ss が可積分なら)。高校で「断面積を積分して体積を求めた」のはこの定理の応用である。例えば二つの円柱 x2+z2≤1x^2 + z^2 \leq 1、y2+z2≤1y^2 + z^2 \leq 1 の共通部分は、高さ zz での断面が一辺 21−z22\sqrt{1 - z^2} の正方形なので、体積は ∫−114(1−z2) dz=163\int_{-1}^{1}4(1 - z^2)\ dz = \frac{16}{3}。

9.5 変数変換公式

1 変数の置換積分 ∫φ(α)φ(β)f(x) dx=∫αβf(φ(t))φ′(t) dt\int_{\varphi(\alpha)}^{\varphi(\beta)}f(x)\ dx = \int_\alpha^\beta f(\varphi(t))\varphi'(t)\ dt の多変数版を考える。φ′(t)\varphi'(t) は「tt の近くの微小区間が φ\varphi で何倍の長さに引き伸ばされるか」を表していた。多変数では、微小な長方形は Φ\Phi によって近似的に DΦD\Phi による像、つまり平行体に写される。その体積の拡大率が ∣det⁡DΦ∣\lvert \det D\Phi \rvert である。

命題 9.12(線形変換と体積)TT を nn 次正則行列、EE をジョルダン可測集合とすると、T(E)T(E) もジョルダン可測で vol⁡(T(E))=∣det⁡T∣vol⁡(E)\operatorname{vol}(T(E)) = \lvert \det T \rvert\operatorname{vol}(E)。

証明の概略. 正則行列は基本行列の積であり、det⁡\det は乗法的なので、TT が基本行列で EE が長方形の場合に示せばよい(一般の EE は長方形の有限和で内側と外側から近似する)。(i) 一つの座標を cc 倍する行列なら、T(E)T(E) は長方形で体積は ∣c∣\lvert c \rvert 倍。(ii) 二つの座標を入れ替える行列なら、T(E)T(E) は同じ体積の長方形。(iii) xi↦xi+cxjx_i \mapsto x_i + cx_j(i≠ji \neq j)なら、T(E)T(E) は長方形ではないが、xix_i 以外の座標を固定した断面は同じ長さの区間を平行移動したものなので、フビニの定理(xix_i について先に積分する)より体積は変わらない。det⁡\det はそれぞれ cc、−1-1、11 である。□\square

定理 9.13(変数変換公式, change of variables formula)U⊂RnU \subset \mathbb{R}^n を開集合、Φ ⁣:U→Rn\Phi\colon U \to \mathbb{R}^n を C1C^1 級の単射で、UU 上 det⁡DΦ≠0\det D\Phi \neq 0 とする。D⊂UD \subset U をコンパクトなジョルダン可測集合、ff を Φ(D)\Phi(D) 上の連続関数とすると、Φ(D)\Phi(D) はジョルダン可測で

∫Φ(D)f(y) dy=∫Df(Φ(x)) ∣det⁡DΦ(x)∣ dx\int_{\Phi(D)}f(y)\,dy = \int_D f(\Phi(x))\,\lvert \det D\Phi(x) \rvert\,dx

さらに、Φ\Phi が DD を含む開集合上で C1C^1 級であれば、単射性と det⁡DΦ≠0\det D\Phi \neq 0 が DD の体積 0 の部分集合(例えば境界 ∂D\partial D)上で破れていても、同じ等式が成り立つ。

∣det⁡DΦ∣\lvert \det D\Phi \rvert をヤコビアン (Jacobian) と呼ぶ(det⁡DΦ\det D\Phi そのものを呼ぶこともある)。記号的には dy=∣det⁡DΦ(x)∣ dxdy = \lvert \det D\Phi(x) \rvert\ dx と覚えればよい。

証明の概略. DD を一辺 hh の小さな立方体 QQ に分割し、各 QQ の中心を cQc_Q とする。Φ\Phi は QQ の上で一次近似 x↦Φ(cQ)+DΦ(cQ)(x−cQ)x \mapsto \Phi(c_Q) + D\Phi(c_Q)(x - c_Q) に近く、その誤差は DΦD\Phi の一様連続性と平均値の不等式(定理 8.4)により o(h)o(h) で一様に抑えられる。したがって、命題 9.12 より

vol⁡(Φ(Q))=∣det⁡DΦ(cQ)∣vol⁡(Q)+o(hn)\operatorname{vol}(\Phi(Q)) = \lvert \det D\Phi(c_Q) \rvert\operatorname{vol}(Q) + o(h^n)

となり、o(hn)o(h^n) は QQ によらず一様である。Φ(D)\Phi(D) は Φ(Q)\Phi(Q) たちに(境界の重なりを除いて)分割されるので

∫Φ(D)f≈∑Qf(Φ(cQ))vol⁡(Φ(Q))≈∑Qf(Φ(cQ))∣det⁡DΦ(cQ)∣vol⁡(Q)≈∫D(f∘Φ)∣det⁡DΦ∣\int_{\Phi(D)}f \approx \sum_Q f(\Phi(c_Q))\operatorname{vol}(\Phi(Q)) \approx \sum_Q f(\Phi(c_Q))\lvert \det D\Phi(c_Q) \rvert\operatorname{vol}(Q) \approx \int_D (f \circ \Phi)\lvert \det D\Phi \rvert

であり、立方体の個数は O(h−n)O(h^{-n}) なので、h→0h \to 0 で誤差の総和は 00 に近づく。Φ(D)\Phi(D) のジョルダン可測性や、DD の境界付近の立方体の扱いなどの細部を詰めるには相当の議論が必要である。完全な証明は、例えば README の参考文献の杉浦『解析入門』や Rudin, Principles of Mathematical Analysis 第10章を参照。□\square

1 変数の置換積分では φ\varphi の単射性は不要だった(向きの付いた積分 ∫αβ\int_\alpha^\beta で打ち消し合う)が、多変数では単射性と絶対値 ∣det⁡∣\lvert \det \rvert が必要である。向きを考えた積分は、第10章と 07-manifolds 第6章 の微分形式の積分で扱う。

9.6 極座標・球座標

極座標 Φ(r,θ)=(rcos⁡θ,rsin⁡θ)\Phi(r, \theta) = (r\cos\theta, r\sin\theta) では

det⁡DΦ=det⁡(cos⁡θ−rsin⁡θsin⁡θrcos⁡θ)=r\det D\Phi = \det\begin{pmatrix} \cos\theta & -r\sin\theta \\ \sin\theta & r\cos\theta \end{pmatrix} = r

なので dx dy=r dr dθdx\ dy = r\ dr\ d\theta。[0,R]×[0,2π][0, R] \times [0, 2\pi] は円板 {x2+y2≤R2}\lbrace x^2 + y^2 \leq R^2 \rbrace に写され、単射性と det⁡≠0\det \neq 0 が破れるのは境界(r=0r = 0 や θ=0,2π\theta = 0, 2\pi)だけなので、定理 9.13 の後半が使える。

例 9.14 (1) 半径 RR の円板の面積は ∫02π∫0Rr dr dθ=πR2\int_0^{2\pi}\int_0^R r\ dr\ d\theta = \pi R^2。 (2) ∫x2+y2≤1(x2+y2) dx dy=∫02π∫01r3 dr dθ=π2\int_{x^2 + y^2 \leq 1}(x^2 + y^2)\ dx\ dy = \int_0^{2\pi}\int_0^1 r^3\ dr\ d\theta = \frac{\pi}{2}。 (3) 放物面 z=x2+y2z = x^2 + y^2 と平面 z=1z = 1 で囲まれた部分の体積は ∫x2+y2≤1(1−x2−y2)=2π∫01(1−r2)r dr=π2\int_{x^2 + y^2 \leq 1}(1 - x^2 - y^2) = 2\pi\int_0^1(1 - r^2)r\ dr = \frac{\pi}{2}。

球座標 Φ(r,ϕ,θ)=(rsin⁡ϕcos⁡θ,rsin⁡ϕsin⁡θ,rcos⁡ϕ)\Phi(r, \phi, \theta) = (r\sin\phi\cos\theta, r\sin\phi\sin\theta, r\cos\phi)(0≤ϕ≤π0 \leq \phi \leq \pi、0≤θ≤2π0 \leq \theta \leq 2\pi)では、第 3 行で余因子展開すると

det⁡DΦ=det⁡(sin⁡ϕcos⁡θrcos⁡ϕcos⁡θ−rsin⁡ϕsin⁡θsin⁡ϕsin⁡θrcos⁡ϕsin⁡θrsin⁡ϕcos⁡θcos⁡ϕ−rsin⁡ϕ0)=r2sin⁡ϕ\det D\Phi = \det\begin{pmatrix} \sin\phi\cos\theta & r\cos\phi\cos\theta & -r\sin\phi\sin\theta \\ \sin\phi\sin\theta & r\cos\phi\sin\theta & r\sin\phi\cos\theta \\ \cos\phi & -r\sin\phi & 0 \end{pmatrix} = r^2\sin\phi

なので dx dy dz=r2sin⁡ϕ dr dϕ dθdx\ dy\ dz = r^2\sin\phi\ dr\ d\phi\ d\theta。例えば半径 RR の球の体積は ∫02π∫0π∫0Rr2sin⁡ϕ dr dϕ dθ=2π⋅2⋅R33=43πR3\int_0^{2\pi}\int_0^\pi\int_0^R r^2\sin\phi\ dr\ d\phi\ d\theta = 2\pi \cdot 2 \cdot \frac{R^3}{3} = \frac{4}{3}\pi R^3。

9.7 広義重積分

非有界な領域や非有界な関数の積分は、1 変数と同様に有界な場合の極限として定める。ただし多変数では「近づき方」がたくさんあるので、それによらないことを確かめる必要がある。

定義 9.15 開集合 D⊂RnD \subset \mathbb{R}^n に対し、コンパクトなジョルダン可測集合の増大列 K1⊂K2⊂⋯⊂DK_1 \subset K_2 \subset \cdots \subset D で ⋃mKm∘=D\bigcup_m K_m^\circ = D を満たすものを、DD の近似列 (exhaustion) という。

近似列 (Km)(K_m) について、DD に含まれる任意のコンパクト集合 CC は、ある KmK_m に含まれる。実際、そうでなければ xm∈C∖Kmx_m \in C \setminus K_m がとれ、収束部分列の極限 x∈C⊂Dx \in C \subset D はある Km0∘K_{m_0}^\circ に属するが、十分先の xmx_m は開集合 Km0∘⊂KmK_{m_0}^\circ \subset K_m に入るので矛盾する。

命題 9.16 f≥0f \geq 0 を開集合 DD 上の連続関数とする。近似列 (Km)(K_m) に対する極限 lim⁡m→∞∫Kmf∈[0,∞]\lim_{m \to \infty}\int_{K_m}f \in [0, \infty] は、近似列のとり方によらない。

証明. ∫Kmf\int_{K_m}f は mm について単調増加なので極限が存在する。別の近似列 (Ll)(L_l) をとると、各 KmK_m はある LlL_l に含まれるので ∫Kmf≤∫Llf≤lim⁡l∫Llf\int_{K_m}f \leq \int_{L_l}f \leq \lim_l\int_{L_l}f。m→∞m \to \infty として lim⁡m∫Kmf≤lim⁡l∫Llf\lim_m\int_{K_m}f \leq \lim_l\int_{L_l}f。逆も同様。□\square

この極限を ∫Df\int_D f と書き、有限のとき広義積分が収束するという。符号の変わる ff については、∫D∣f∣<∞\int_D\lvert f \rvert < \infty(絶対収束)のときに、f±=max⁡{±f,0}f^{\pm} = \max\lbrace \pm f, 0 \rbrace を使って ∫Df=∫Df+−∫Df−\int_D f = \int_D f^+ - \int_D f^- と定める。9.2 節の警告の例が示すように、絶対収束しない場合は近似の仕方で値が変わりうるので、多変数では絶対収束を要求するのが自然である。

9.8 ガウス積分

広義重積分の最も有名な応用は、1 変数の積分 ∫−∞∞e−x2dx\int_{-\infty}^{\infty}e^{-x^2}dx を 2 変数の積分に「持ち上げて」計算する方法である。

定理 9.17(ガウス積分, Gaussian integral)

∫−∞∞e−x2 dx=π\int_{-\infty}^{\infty}e^{-x^2}\,dx = \sqrt{\pi}

証明. f(x,y)=e−(x2+y2)>0f(x, y) = e^{-(x^2 + y^2)} > 0 の R2\mathbb{R}^2 上の広義積分を、二つの近似列で計算する。

(円板) Km={x2+y2≤m2}K_m = \lbrace x^2 + y^2 \leq m^2 \rbrace とすると、極座標で

∫Kmf=∫02π∫0me−r2r dr dθ=2π[−e−r22]0m=π(1−e−m2)→π\int_{K_m}f = \int_0^{2\pi}\int_0^m e^{-r^2}r\,dr\,d\theta = 2\pi\left[-\frac{e^{-r^2}}{2}\right]_0^m = \pi\left(1 - e^{-m^2}\right) \to \pi

(正方形) Lm=[−m,m]2L_m = [-m, m]^2 とすると、フビニの定理より

∫Lmf=∫−mm∫−mme−x2e−y2 dy dx=(∫−mme−x2dx)2→G2\int_{L_m}f = \int_{-m}^{m}\int_{-m}^{m}e^{-x^2}e^{-y^2}\,dy\,dx = \left(\int_{-m}^{m}e^{-x^2}dx\right)^2 \to G^2

ここで G=∫−∞∞e−x2dxG = \int_{-\infty}^{\infty}e^{-x^2}dx(e−x2≤e−∣x∣+1e^{-x^2} \leq e^{-\lvert x \rvert + 1} より収束する)。命題 9.16 より G2=πG^2 = \pi、G>0G > 0 なので G=πG = \sqrt{\pi}。□\square

x=y/2x = y/\sqrt{2} などと置換すると、確率論で基本的な正規分布の密度の積分 ∫−∞∞12πe−y2/2dy=1\int_{-\infty}^{\infty}\frac{1}{\sqrt{2\pi}}e^{-y^2/2}dy = 1 が得られる(11-probability 第4章)。また x=tx = \sqrt{t} と置換すると Γ(12)=∫0∞t−1/2e−tdt=2∫0∞e−x2dx=π\Gamma\left(\frac{1}{2}\right) = \int_0^\infty t^{-1/2}e^{-t}dt = 2\int_0^\infty e^{-x^2}dx = \sqrt{\pi}。

9.9 ガンマ関数とベータ関数、nn 次元球の体積

第5章で予告した関係式を証明する。

定理 9.18 p,q>0p, q > 0 に対して B(p,q)=Γ(p)Γ(q)Γ(p+q)B(p, q) = \dfrac{\Gamma(p)\Gamma(q)}{\Gamma(p + q)}。

証明. D=(0,∞)2D = (0, \infty)^2 上の f(x,y)=xp−1yq−1e−(x+y)>0f(x, y) = x^{p-1}y^{q-1}e^{-(x+y)} > 0 の広義積分を考える。近似列 Lm=[1m,m]2L_m = [\frac{1}{m}, m]^2 についてはフビニの定理より

∫Lmf=∫1/mmxp−1e−xdx⋅∫1/mmyq−1e−ydy→Γ(p)Γ(q)\int_{L_m}f = \int_{1/m}^{m}x^{p-1}e^{-x}dx \cdot \int_{1/m}^{m}y^{q-1}e^{-y}dy \to \Gamma(p)\Gamma(q)

一方、Φ(u,v)=(uv,u(1−v))\Phi(u, v) = (uv, u(1 - v)) は (0,∞)×(0,1)(0, \infty) \times (0, 1) から DD への C1C^1 級全単射(逆写像は u=x+yu = x + y、v=xx+yv = \frac{x}{x+y})で、

det⁡DΦ=det⁡(vu1−v−u)=−u≠0\det D\Phi = \det\begin{pmatrix} v & u \\ 1 - v & -u \end{pmatrix} = -u \neq 0

Km=Φ([1m,m]×[1m,1−1m])K_m = \Phi\left([\frac{1}{m}, m] \times [\frac{1}{m}, 1 - \frac{1}{m}]\right) は DD の近似列であり、変数変換公式とフビニの定理より

∫Kmf=∫1/mm∫1/m1−1/m(uv)p−1(u(1−v))q−1e−uu dv du=∫1/mmup+q−1e−udu⋅∫1/m1−1/mvp−1(1−v)q−1dv\int_{K_m}f = \int_{1/m}^{m}\int_{1/m}^{1 - 1/m}(uv)^{p-1}(u(1-v))^{q-1}e^{-u}u\,dv\,du = \int_{1/m}^{m}u^{p+q-1}e^{-u}du \cdot \int_{1/m}^{1-1/m}v^{p-1}(1-v)^{q-1}dv

これは Γ(p+q)B(p,q)\Gamma(p + q)B(p, q) に収束する。命題 9.16 より Γ(p)Γ(q)=Γ(p+q)B(p,q)\Gamma(p)\Gamma(q) = \Gamma(p+q)B(p,q)。□\square

半径 rr の nn 次元球 Bn(r)={x∈Rn∣∥x∥≤r}B^n(r) = \lbrace x \in \mathbb{R}^n \mid \lVert x \rVert \leq r \rbrace の体積を求めよう。線形変換 x↦rxx \mapsto rx(det⁡=rn\det = r^n)と命題 9.12 より vol⁡(Bn(r))=Vnrn\operatorname{vol}(B^n(r)) = V_nr^n、ただし Vn=vol⁡(Bn(1))V_n = \operatorname{vol}(B^n(1))。

定理 9.19(nn 次元球の体積)

Vn=πn/2Γ(n2+1)V_n = \frac{\pi^{n/2}}{\Gamma\left(\frac{n}{2} + 1\right)}

証明. nn についての帰納法。n=1n = 1 では V1=2V_1 = 2、右辺は πΓ(3/2)=π12π=2\frac{\sqrt{\pi}}{\Gamma(3/2)} = \frac{\sqrt{\pi}}{\frac{1}{2}\sqrt{\pi}} = 2。n≥2n \geq 2 のとき、Bn(1)B^n(1) を最後の座標 tt で切ると、断面は半径 1−t2\sqrt{1 - t^2} の (n−1)(n-1) 次元球なので、フビニの定理より

Vn=∫−11Vn−1(1−t2)n−12dt=Vn−1⋅2∫01(1−t2)n−12dt=Vn−1B(12,n+12)V_n = \int_{-1}^{1}V_{n-1}(1 - t^2)^{\frac{n-1}{2}}dt = V_{n-1} \cdot 2\int_0^1(1 - t^2)^{\frac{n-1}{2}}dt = V_{n-1}B\left(\frac{1}{2}, \frac{n+1}{2}\right)

(最後は t2=st^2 = s と置換した)。定理 9.18 と Γ(12)=π\Gamma\left(\frac{1}{2}\right) = \sqrt{\pi} より

Vn=Vn−1π Γ(n+12)Γ(n2+1)V_n = V_{n-1}\frac{\sqrt{\pi}\,\Gamma\left(\frac{n+1}{2}\right)}{\Gamma\left(\frac{n}{2} + 1\right)}

帰納法の仮定 Vn−1=π(n−1)/2Γ(n+12)V_{n-1} = \frac{\pi^{(n-1)/2}}{\Gamma\left(\frac{n+1}{2}\right)} を代入すれば主張を得る。□\square

具体的には V1=2V_1 = 2、V2=πV_2 = \pi、V3=43πV_3 = \frac{4}{3}\pi、V4=π22V_4 = \frac{\pi^2}{2}、V5=815π2V_5 = \frac{8}{15}\pi^2、V6=π36V_6 = \frac{\pi^3}{6} である。数値では V5≈5.264V_5 \approx 5.264 が最大で、その後 VnV_n は減少して Vn→0V_n \to 0(スターリングの公式より超指数的に)となる。高次元の単位球は、それを囲む一辺 22 の立方体(体積 2n2^n)に比べて極めて小さい。また、半径 rr の球面 Sn−1S^{n-1} の「面積」は ddrVnrn=nVnrn−1\frac{d}{dr}V_nr^n = nV_nr^{n-1} であり、n=3n = 3 で 4πr24\pi r^2 となる。

まとめ

  • 長方形上の積分はダルブー和で定義され、連続関数は可積分である。
  • フビニの定理:可積分な ff の重積分は累次積分に等しい。可積分性がないと積分の順序交換は失敗しうる。
  • 境界が体積 0 の集合がジョルダン可測集合であり、その上の有界連続関数は可積分。連続関数のグラフは体積 0 なので、縦線集合上の積分は累次積分で計算できる。
  • 変数変換公式 ∫Φ(D)f=∫D(f∘Φ)∣det⁡DΦ∣\int_{\Phi(D)}f = \int_D(f \circ \Phi)\lvert \det D\Phi \rvert において、∣det⁡DΦ∣\lvert \det D\Phi \rvert は体積の局所的な拡大率である。線形変換では vol⁡(T(E))=∣det⁡T∣vol⁡(E)\operatorname{vol}(T(E)) = \lvert \det T \rvert\operatorname{vol}(E)。
  • 極座標で dx dy=r dr dθdx\ dy = r\ dr\ d\theta、球座標で dx dy dz=r2sin⁡ϕ dr dϕ dθdx\ dy\ dz = r^2\sin\phi\ dr\ d\phi\ d\theta。
  • 正値関数の広義重積分は近似列によらない。これを使ってガウス積分 ∫e−x2=π\int e^{-x^2} = \sqrt{\pi} と B(p,q)=Γ(p)Γ(q)Γ(p+q)B(p,q) = \frac{\Gamma(p)\Gamma(q)}{\Gamma(p+q)} が示せる。
  • nn 次元単位球の体積は Vn=πn/2Γ(n/2+1)V_n = \frac{\pi^{n/2}}{\Gamma(n/2 + 1)}。

演習問題

問題 9.1 ★ 次を計算せよ。 (1) ∫Dxy dx dy\int_D xy\ dx\ dy、D={0≤y≤x≤1}D = \lbrace 0 \leq y \leq x \leq 1 \rbrace (2) ∫01(∫x1ey3dy)dx\int_0^1\left(\int_{\sqrt{x}}^{1}e^{y^3}dy\right)dx

解答

(1) ∫01∫0xxy dy dx=∫01x32dx=18\int_0^1\int_0^x xy\ dy\ dx = \int_0^1\frac{x^3}{2}dx = \frac{1}{8}。

(2) 領域は {0≤x≤y2, 0≤y≤1}\lbrace 0 \leq x \leq y^2,\ 0 \leq y \leq 1 \rbrace なので、順序を入れ替えて ∫01∫0y2ey3dx dy=∫01y2ey3dy=e−13\int_0^1\int_0^{y^2}e^{y^3}dx\ dy = \int_0^1 y^2e^{y^3}dy = \frac{e - 1}{3}。

問題 9.2 ★ 楕円体 x2a2+y2b2+z2c2≤1\frac{x^2}{a^2} + \frac{y^2}{b^2} + \frac{z^2}{c^2} \leq 1 の体積を求めよ。

解答

線形変換 (u,v,w)↦(au,bv,cw)(u, v, w) \mapsto (au, bv, cw)(行列式 abcabc)は単位球をこの楕円体に写すので、命題 9.12 より体積は abc⋅43π=43πabcabc \cdot \frac{4}{3}\pi = \frac{4}{3}\pi abc。

問題 9.3 ★★ (1) ∫R2dx dy(1+x2+y2)α\int_{\mathbb{R}^2}\frac{dx\ dy}{(1 + x^2 + y^2)^\alpha} が収束する α\alpha の範囲と値を求めよ。(2) ∫0<x2+y2<1dx dy(x2+y2)α/2\int_{0 < x^2 + y^2 < 1}\frac{dx\ dy}{(x^2 + y^2)^{\alpha/2}} が収束する α\alpha の範囲と値を求めよ。

解答

(1) 円板の近似列と極座標で 2π∫0mr dr(1+r2)α2\pi\int_0^m\frac{r\ dr}{(1 + r^2)^\alpha}。α≠1\alpha \neq 1 なら =π1−α((1+m2)1−α−1)= \frac{\pi}{1 - \alpha}\left((1 + m^2)^{1-\alpha} - 1\right) で、α>1\alpha > 1 のとき πα−1\frac{\pi}{\alpha - 1} に収束、α<1\alpha < 1 で発散。α=1\alpha = 1 なら πlog⁡(1+m2)→∞\pi\log(1 + m^2) \to \infty。よって α>1\alpha > 1 で収束し値は πα−1\frac{\pi}{\alpha - 1}。

(2) 開集合 D={0<x2+y2<1}D = \lbrace 0 < x^2 + y^2 < 1 \rbrace の近似列として円環 Km={1m+2≤x2+y2≤1−1m+2}K_m = \lbrace \frac{1}{m+2} \leq \sqrt{x^2 + y^2} \leq 1 - \frac{1}{m+2} \rbrace をとる(KmK_m はコンパクトで増大し、内部の和集合は DD)。極座標で ∫Km=2π∫1/(m+2)1−1/(m+2)r1−αdr→2π∫01r1−αdr\int_{K_m} = 2\pi\int_{1/(m+2)}^{1 - 1/(m+2)}r^{1 - \alpha}dr \to 2\pi\int_0^1 r^{1-\alpha}dr であり、α<2\alpha < 2 で 2π2−α\frac{2\pi}{2 - \alpha} に収束し、α≥2\alpha \geq 2 で発散する。一般に Rn\mathbb{R}^n では、原点の近くで ∥x∥−α\lVert x \rVert^{-\alpha} が可積分 ⇔α<n\Leftrightarrow \alpha < n、無限遠で可積分 ⇔α>n\Leftrightarrow \alpha > n である。

問題 9.4 ★★ a>0a > 0、b∈Rb \in \mathbb{R} とする。次を示せ。 (1) ∫0∞x2e−x2dx=π4\int_0^\infty x^2e^{-x^2}dx = \frac{\sqrt{\pi}}{4} (2) ∫−∞∞e−ax2+bxdx=πa eb2/(4a)\int_{-\infty}^{\infty}e^{-ax^2 + bx}dx = \sqrt{\frac{\pi}{a}}\ e^{b^2/(4a)}

解答

(1) x=tx = \sqrt{t} と置換して 12∫0∞t1/2e−tdt=12Γ(32)=12⋅12π=π4\frac{1}{2}\int_0^\infty t^{1/2}e^{-t}dt = \frac{1}{2}\Gamma\left(\frac{3}{2}\right) = \frac{1}{2}\cdot\frac{1}{2}\sqrt{\pi} = \frac{\sqrt{\pi}}{4}。

(2) −ax2+bx=−a(x−b2a)2+b24a-ax^2 + bx = -a\left(x - \frac{b}{2a}\right)^2 + \frac{b^2}{4a} と平方完成し、y=a(x−b2a)y = \sqrt{a}\left(x - \frac{b}{2a}\right) と置換すると eb2/(4a)1a∫−∞∞e−y2dy=πaeb2/(4a)e^{b^2/(4a)}\frac{1}{\sqrt{a}}\int_{-\infty}^{\infty}e^{-y^2}dy = \sqrt{\frac{\pi}{a}}e^{b^2/(4a)}。

問題 9.5 ★★ (1) n∈Z≥0n \in \mathbb{Z}_{\geq 0} について Γ(n+12)=(2n)!4nn!π\Gamma\left(n + \frac{1}{2}\right) = \frac{(2n)!}{4^nn!}\sqrt{\pi} を示せ。(2) p>−1p > -1 について ∫0π/2sin⁡pθ dθ=π2⋅Γ(p+12)Γ(p2+1)\int_0^{\pi/2}\sin^p\theta\ d\theta = \frac{\sqrt{\pi}}{2}\cdot\frac{\Gamma\left(\frac{p+1}{2}\right)}{\Gamma\left(\frac{p}{2} + 1\right)} を示し、例 5.16 と比較せよ。

解答

(1) Γ(s+1)=sΓ(s)\Gamma(s + 1) = s\Gamma(s) より Γ(n+12)=2n−12⋅2n−32⋯12π=(2n−1)(2n−3)⋯12nπ\Gamma\left(n + \frac{1}{2}\right) = \frac{2n-1}{2}\cdot\frac{2n-3}{2}\cdots\frac{1}{2}\sqrt{\pi} = \frac{(2n-1)(2n-3)\cdots 1}{2^n}\sqrt{\pi}。(2n−1)(2n−3)⋯1=(2n)!2nn!(2n-1)(2n-3)\cdots 1 = \frac{(2n)!}{2^nn!} なので主張を得る。

(2) 定理 5.24 (2) で 2p′−1=p2p' - 1 = p、2q′−1=02q' - 1 = 0 とすると ∫0π/2sin⁡pθ dθ=12B(p+12,12)\int_0^{\pi/2}\sin^p\theta\ d\theta = \frac{1}{2}B\left(\frac{p+1}{2}, \frac{1}{2}\right)。定理 9.18 と Γ(12)=π\Gamma\left(\frac{1}{2}\right) = \sqrt{\pi} から主張を得る。例えば p=2mp = 2m なら (1) より π2⋅(2m)!π/(4mm!)m!=π2⋅(2m)!4m(m!)2\frac{\sqrt{\pi}}{2}\cdot\frac{(2m)!\sqrt{\pi}/(4^mm!)}{m!} = \frac{\pi}{2}\cdot\frac{(2m)!}{4^m(m!)^2} で、例 5.16 の I2mI_{2m} と一致する。

問題 9.6 ★★ 0<a<R0 < a < R とする。xzxz 平面の円板 (x−R)2+z2≤a2(x - R)^2 + z^2 \leq a^2 を zz 軸のまわりに回転してできるトーラスの体積が 2π2Ra22\pi^2Ra^2 であることを示せ。

解答

円柱座標 (x,y,z)=(ρcos⁡θ,ρsin⁡θ,z)(x, y, z) = (\rho\cos\theta, \rho\sin\theta, z)(ヤコビアン ρ\rho)を使うと、トーラスは {(ρ−R)2+z2≤a2}\lbrace (\rho - R)^2 + z^2 \leq a^2 \rbrace なので、体積は ∫02π∫Eρ dρ dz dθ=2π∫Eρ dρ dz\int_0^{2\pi}\int_{E}\rho\ d\rho\ dz\ d\theta = 2\pi\int_E\rho\ d\rho\ dz(EE は ρz\rho z 平面の円板)。ρ=R+s\rho = R + s とすると ∫Eρ=R⋅πa2+∫s2+z2≤a2s=πRa2\int_E\rho = R\cdot\pi a^2 + \int_{s^2 + z^2 \leq a^2}s = \pi Ra^2(第 2 項は対称性より 00)。よって 2π2Ra22\pi^2Ra^2。(体積 = 断面積 × 重心の描く円周の長さ、というパップス–ギュルダンの定理の一例。)

問題 9.7 ★★ 3 次元の正四面体 T={x,y,z≥0, x+y+z≤1}T = \lbrace x, y, z \geq 0,\ x + y + z \leq 1 \rbrace について、vol⁡(T)=16\operatorname{vol}(T) = \frac{1}{6} と ∫Tx dx dy dz=124\int_T x\ dx\ dy\ dz = \frac{1}{24} を示し、TT の重心を求めよ。

解答

vol⁡(T)=∫01∫01−x(1−x−y) dy dx=∫01(1−x)22dx=16\operatorname{vol}(T) = \int_0^1\int_0^{1-x}(1 - x - y)\ dy\ dx = \int_0^1\frac{(1 - x)^2}{2}dx = \frac{1}{6}。∫Tx=∫01x(1−x)22dx=12B(2,3)=12⋅1! 2!4!=124\int_T x = \int_0^1 x\frac{(1-x)^2}{2}dx = \frac{1}{2}B(2, 3) = \frac{1}{2}\cdot\frac{1!\ 2!}{4!} = \frac{1}{24}。対称性から重心は (1/241/6,14,14)=(14,14,14)\left(\frac{1/24}{1/6}, \frac{1}{4}, \frac{1}{4}\right) = \left(\frac{1}{4}, \frac{1}{4}, \frac{1}{4}\right)。

問題 9.8 ★★★ (1) ∫Rne−∥x∥2dx=πn/2\int_{\mathbb{R}^n}e^{-\lVert x \rVert^2}dx = \pi^{n/2} を示せ。(2) 連続関数 g ⁣:[0,∞)→[0,∞)g\colon [0, \infty) \to [0, \infty) について ∫∥x∥≤Rg(∥x∥) dx=nVn∫0Rg(r)rn−1dr\int_{\lVert x \rVert \leq R}g(\lVert x \rVert)\ dx = nV_n\int_0^R g(r)r^{n-1}dr を示せ。(3) (1)、(2) から Vn=πn/2Γ(n/2+1)V_n = \frac{\pi^{n/2}}{\Gamma(n/2 + 1)} を再び導け。

解答

(1) 立方体の近似列 [−m,m]n[-m, m]^n とフビニの定理より (∫−mme−t2dt)n→πn/2\left(\int_{-m}^{m}e^{-t^2}dt\right)^n \to \pi^{n/2}。

(2) [0,R][0, R] の分割 0=r0<⋯<rk=R0 = r_0 < \cdots < r_k = R をとると、球殻 {ri−1≤∥x∥≤ri}\lbrace r_{i-1} \leq \lVert x \rVert \leq r_i \rbrace の体積は Vn(rin−ri−1n)=Vnnξin−1ΔriV_n(r_i^n - r_{i-1}^n) = V_n n\xi_i^{n-1}\Delta r_i(平均値の定理、ξi∈(ri−1,ri)\xi_i \in (r_{i-1}, r_i))。その上で g(∥x∥)g(\lVert x \rVert) は gg の [ri−1,ri][r_{i-1}, r_i] 上の最小値 mim_i と最大値 MiM_i の間にあるので

∑iminVnξin−1Δri≤∫∥x∥≤Rg(∥x∥) dx≤∑iMinVnξin−1Δri\sum_i m_inV_n\xi_i^{n-1}\Delta r_i \leq \int_{\lVert x \rVert \leq R}g(\lVert x \rVert)\,dx \leq \sum_i M_inV_n\xi_i^{n-1}\Delta r_i

分割を細かくすると、両辺は gg の一様連続性によりリーマン和 ∑g(ξi)nVnξin−1Δri\sum g(\xi_i)nV_n\xi_i^{n-1}\Delta r_i との差が 00 に近づき、nVn∫0Rg(r)rn−1drnV_n\int_0^R g(r)r^{n-1}dr に収束する。

(3) (2) で g(r)=e−r2g(r) = e^{-r^2}、R→∞R \to \infty とすると、(1) より πn/2=nVn∫0∞e−r2rn−1dr=nVn⋅12Γ(n2)=VnΓ(n2+1)\pi^{n/2} = nV_n\int_0^\infty e^{-r^2}r^{n-1}dr = nV_n\cdot\frac{1}{2}\Gamma\left(\frac{n}{2}\right) = V_n\Gamma\left(\frac{n}{2} + 1\right)(r=tr = \sqrt{t} と置換した)。

この章を読み終えたら

「読了」にすると学習記録とロードマップに反映されます。演習の自己採点もお忘れなく。

この章の誤りを報告GitHub で見る