## 2. 準備：記号と前提この記事では次の記号を使います。 - $R = [a,b]\times[c,d]$（$a < b$, $c < d$）を $\mathbb{R}^2$ の閉長方形とし、その面積を $\lvert R\rvert = (b-a)(d-c)$ と書きます。 - 点は $P=(x,y)$ のように書き、$\lVert P - Q\rVert$ はユークリッド距離です。 - $f$ が集合 $S$ 上**有界**とは、ある $K$ があって $S$ 上つねに $\lvert f\rvert \le K$ となることをいいます。前提として次の 3 つを使います。 1. 1 変数のリーマン積分（定義、区間加法性、単調性、線形性）。[積分の基本定理と定積分](/mathematics/calculus/integration-and-ftc) のを参照してください。 2. 有界閉集合上の連続関数は一様連続であり、最大値・最小値を取ること（ハイネ・カントールの定理、最大値原理）。[極限と連続性 (ε-δ論法)](/mathematics/calculus/limits-and-continuity) のを参照してください。 3. 偏微分とヤコビ行列。[多変数関数の微分と偏微分](/mathematics/calculus/multivariable-differentiation) のを参照してください。 ## 3. 長方形上の 2 重積分まず、いちばん素直な領域である長方形の上で積分を定義します。一般の領域は第 5 節で扱います。 $R=[a,b]\times[c,d]$ とする。$[a,b]$ の分点 $a = x_0 < x_1 < \cdots < x_m = b$ と $[c,d]$ の分点 $c = y_0 < y_1 < \cdots < y_n = d$ の組を $R$ の**分割** $\Delta$ と呼び、$R$ を小長方形 $$ R_{ij} = [x_{i-1}, x_i]\times[y_{j-1}, y_j] \qquad (1\le i\le m,\ 1\le j\le n) $$ に分ける。$\Delta x_i = x_i - x_{i-1}$、$\Delta y_j = y_j - y_{j-1}$、$\lvert R_{ij}\rvert = \Delta x_i \Delta y_j$ と書く。また $$ \lvert \Delta\rvert = \max_{i,j}\sqrt{(\Delta x_i)^2 + (\Delta y_j)^2} $$ （小長方形の対角線の長さの最大値）を分割の**幅**という。幅を対角線で測っておくと、「同じ小長方形に属する 2 点の距離は $\lvert\Delta\rvert$ 以下」がそのまま言えるので、あとの一様連続性の議論が楽になります。 $f\colon R\to\mathbb{R}$ を有界関数とする。分割 $\Delta$ と各小長方形から選んだ**代表点** $\xi_{ij}\in R_{ij}$ に対し $$ \Sigma(f,\Delta,\xi) = \sum_{i=1}^{m}\sum_{j=1}^{n} f(\xi_{ij})\,\Delta x_i\,\Delta y_j $$ を**リーマン和**という。ある実数 $I$ が存在して、任意の $\varepsilon > 0$ に対しある $\delta > 0$ があり、$\lvert\Delta\rvert < \delta$ を満たす**すべての**分割 $\Delta$ と**すべての**代表点の選び方 $\xi$ について $$ \lvert \Sigma(f,\Delta,\xi) - I\rvert < \varepsilon $$ が成り立つとき、$f$ は $R$ 上**（リーマン）可積分**であるといい、$I$ を $f$ の $R$ 上の **2 重積分**と呼んで $$ \iint_R f(x,y)\,dx\,dy = I $$ と書く。 $f \ge 0$ のとき、$f(\xi_{ij})\Delta x_i\Delta y_j$ は底面 $R_{ij}$、高さ $f(\xi_{ij})$ の直方体の体積です。リーマン和はこの柱をすべて足したもの、2 重積分はそれを細かくした極限で、これが「グラフの下の立体の体積」の定義になります。定義はこれでよいのですが、このままでは「すべての分割とすべての代表点」を相手にしなければならず、可積分性の判定に使えません。1 変数のときと同じく、代表点を $\inf$ と $\sup$ に振り切った和を考えると扱いやすくなります。分割 $\Delta$ に対し $$ m_{ij} = \inf_{R_{ij}} f,\quad M_{ij} = \sup_{R_{ij}} f,\qquad s(f,\Delta) = \sum_{i,j} m_{ij}\lvert R_{ij}\rvert,\quad S(f,\Delta) = \sum_{i,j} M_{ij}\lvert R_{ij}\rvert $$ とおき、それぞれ**不足和**（下ダルブー和）、**過剰和**（上ダルブー和）と呼びます。 $f\colon R\to\mathbb{R}$ を有界関数とする。$f$ が $R$ 上可積分であるための必要十分条件は、任意の $\varepsilon > 0$ に対して $$ S(f,\Delta) - s(f,\Delta) < \varepsilon $$ を満たす $R$ の分割 $\Delta$ が存在することである。またこのとき $$ \iint_R f(x,y)\,dx\,dy = \sup_{\Delta} s(f,\Delta) = \inf_{\Delta} S(f,\Delta) $$ が成り立つ。 **準備。** 任意の分割 $\Delta$ と代表点 $\xi$ に対し、$\xi_{ij}\in R_{ij}$ だから $m_{ij}\le f(\xi_{ij})\le M_{ij}$ であり、正の数 $\lvert R_{ij}\rvert$ を掛けて足せば $$ s(f,\Delta) \le \Sigma(f,\Delta,\xi)\le S(f,\Delta) $$ を得ます。さらに $\eta > 0$ を任意に取ると、$\inf$ の定義から各 $(i,j)$ について $f(\xi_{ij}) < m_{ij}+\eta$ となる $\xi_{ij}\in R_{ij}$ が選べ、そのとき $\Sigma < s(f,\Delta) + \eta\lvert R\rvert$ です。$\eta$ は任意なので、代表点をすべて動かしたときのリーマン和の下限はちょうど $s(f,\Delta)$ です。同様に上限は $S(f,\Delta)$ です。 **（必要性）** $f$ が可積分で $I = \iint_R f$ とします。$\varepsilon > 0$ に対しの $\delta$ を $\varepsilon/3$ に対して取り、$\lvert\Delta\rvert < \delta$ なる分割を 1 つ固定します。すべての代表点で $I - \varepsilon/3 < \Sigma < I + \varepsilon/3$ ですから、代表点について上限・下限を取れば、上の準備より $S(f,\Delta)\le I+\varepsilon/3$ かつ $s(f,\Delta)\ge I-\varepsilon/3$。よって $S - s \le 2\varepsilon/3 < \varepsilon$ です。 **（十分性）** まず細分についての単調性を確かめます。分割 $\Delta'$ が $\Delta$ の**細分**（分点を追加して得られる分割）であるとし、$\Delta$ の小長方形 $R_{ij}$ が $\Delta'$ の小長方形 $R'_1,\dots,R'_k$ に分かれたとします。$R'_l\subset R_{ij}$ なので $\inf_{R'_l} f \ge \inf_{R_{ij}} f = m_{ij}$ であり、$\sum_l \lvert R'_l\rvert = \lvert R_{ij}\rvert$ だから $$ \sum_{l}\Bigl(\inf_{R'_l} f\Bigr)\lvert R'_l\rvert \ \ge\ m_{ij}\lvert R_{ij}\rvert . $$ $(i,j)$ について足せば $s(f,\Delta')\ge s(f,\Delta)$、同様に $S(f,\Delta')\le S(f,\Delta)$ です。ゆえに任意の 2 つの分割 $\Delta_1,\Delta_2$ に対し、両方の分点を合わせた共通細分 $\Delta_1\vee\Delta_2$ を挟んで $$ s(f,\Delta_1)\le s(f,\Delta_1\vee\Delta_2)\le S(f,\Delta_1\vee\Delta_2)\le S(f,\Delta_2) $$ が成り立ち、$\underline{I} := \sup_\Delta s(f,\Delta) \le \inf_\Delta S(f,\Delta) =: \overline{I}$ を得ます。仮定よりどんな $\varepsilon>0$ に対しても $\overline{I}-\underline{I}\le S(f,\Delta)-s(f,\Delta)<\varepsilon$ となる $\Delta$ が取れるので、$\underline I = \overline I$。これを $J$ とおきます。残るは、リーマン和が $\lvert\Delta\rvert\to 0$ で $J$ に収束することです。$K = \sup_R \lvert f\rvert$ とおき（$f$ は有界）、$\varepsilon>0$ に対し $S(f,\Delta_0)-s(f,\Delta_0)<\varepsilon/2$ となる $\Delta_0$ を取り、その内部の分点が $x$ 方向に $p$ 本、$y$ 方向に $q$ 本あるとします。$C = p(d-c)+q(b-a)$ とおきます。任意の分割 $\Delta$ に対し共通細分 $\Delta\vee\Delta_0$ を考えると、$\Delta$ の小長方形のうち実際に切られるのは $\Delta_0$ の分点を通る直線と交わるものだけで、それらは幅 $\lvert\Delta\rvert$ 以下の帯に含まれるので、面積の合計は $C\lvert\Delta\rvert$ 以下です。切られた小長方形 1 つあたり過剰和の差は $2K\lvert R_{ij}\rvert$ 以下（$\sup$ の値はどれも絶対値 $K$ 以下）だから $$ S(f,\Delta) - S(f,\Delta\vee\Delta_0) \le 2KC\lvert\Delta\rvert . $$ そこで $\delta = \dfrac{\varepsilon}{8KC+1}$ と取れば、$\lvert\Delta\rvert<\delta$ のとき右辺は $\varepsilon/4$ 未満で、細分の単調性 $S(f,\Delta\vee\Delta_0)\le S(f,\Delta_0)$ と合わせて $$ S(f,\Delta) \le S(f,\Delta_0) + \frac{\varepsilon}{4} < s(f,\Delta_0)+\frac{\varepsilon}{2}+\frac{\varepsilon}{4} \le J + \frac{3\varepsilon}{4} $$ （最後は $s(f,\Delta_0)\le \underline I = J$）。同じ議論で $s(f,\Delta) > J - 3\varepsilon/4$ です。準備の不等式 $s\le\Sigma\le S$ より、$\lvert\Delta\rvert<\delta$ なるすべての分割と代表点で $\lvert \Sigma - J\rvert < \varepsilon$。よって $f$ は可積分で積分値は $J$ です。 $R = [a,b]\times[c,d]$（$a $R$ は有界閉集合なので、$f$ は $R$ 上一様連続です（第 2 節の前提 2）。$\varepsilon>0$ を与え、 $$ \lVert P-Q\rVert < \delta \implies \lvert f(P)-f(Q)\rvert < \frac{\varepsilon}{2\lvert R\rvert} $$ となる $\delta>0$ を取ります。$\lvert\Delta\rvert<\delta$ なる分割を 1 つ取ると、の幅の定義より、同じ小長方形 $R_{ij}$ に属する 2 点の距離は $\lvert\Delta\rvert$ 以下、したがって $\delta$ 未満です。また $R_{ij}$ は有界閉集合なので $f$ はそこで最大値と最小値を取り、$M_{ij}=f(Q_{ij})$、$m_{ij}=f(P_{ij})$ となる点 $P_{ij},Q_{ij}\in R_{ij}$ が存在します。ゆえに $$ M_{ij}-m_{ij} = f(Q_{ij})-f(P_{ij}) < \frac{\varepsilon}{2\lvert R\rvert} $$ であり、$\sum_{i,j}\lvert R_{ij}\rvert = \lvert R\rvert$ を使って $$ S(f,\Delta)-s(f,\Delta) = \sum_{i,j}(M_{ij}-m_{ij})\lvert R_{ij}\rvert \le \frac{\varepsilon}{2\lvert R\rvert}\cdot\lvert R\rvert = \frac{\varepsilon}{2} < \varepsilon . $$ より $f$ は可積分です。 $f(x,y)=xy$、$R=[0,1]\times[0,1]$ とします。$f$ は連続なのでより可積分です。可積分だとわかっていれば、幅が $0$ に近づく**都合のよい**分割・代表点の列を 1 つ選んで極限を計算すれば、それが積分値です。 $x$ 方向も $y$ 方向も $n$ 等分し、代表点を小長方形の右上頂点 $\xi_{ij}=(i/n,\ j/n)$ に取ります。$\Delta x_i = \Delta y_j = 1/n$ なので $$ \Sigma = \sum_{i=1}^{n}\sum_{j=1}^{n}\frac{i}{n}\cdot\frac{j}{n}\cdot\frac{1}{n^2} = \frac{1}{n^4}\Bigl(\sum_{i=1}^{n} i\Bigr)\Bigl(\sum_{j=1}^{n} j\Bigr) = \frac{1}{n^4}\cdot\frac{n(n+1)}{2}\cdot\frac{n(n+1)}{2} = \frac{(n+1)^2}{4n^2}. $$ 分割の幅は $\sqrt{2}/n\to 0$ で、$\dfrac{(n+1)^2}{4n^2}\to\dfrac14$ です。したがって $$ \iint_{[0,1]^2} xy\,dx\,dy = \frac14 . $$ 可積分性を先に確かめた点が大事です。可積分性を知らずに 1 つの列だけを計算しても、他の分割で別の値になる可能性を排除できません。 $f,g$ が $R$ 上可積分、$\alpha,\beta\in\mathbb{R}$ のとき、次が成り立ちます。 - 線形性：$\alpha f+\beta g$ は可積分で $\iint_R(\alpha f+\beta g) = \alpha\iint_R f+\beta\iint_R g$。 - 単調性：$R$ 上つねに $f\le g$ なら $\iint_R f\le \iint_R g$。とくに $f\ge0$ なら $\iint_R f\ge0$。 - 領域についての加法性：$R$ を有限個の小長方形に分けたとき、全体の積分は各部分の積分の和。いずれもリーマン和のレベルで対応する等式・不等式（$\Sigma(\alpha f+\beta g)=\alpha\Sigma(f)+\beta\Sigma(g)$、$f\le g \Rightarrow \Sigma(f)\le\Sigma(g)$）が成り立ち、の極限を取れば従います。詳しい証明は杉浦光夫『解析入門 II』の重積分の章を参照してください。以下ではこれらを自由に使います。 ## 4. フビニの定理と累次積分 $f\colon R=[a,b]\times[c,d]\to\mathbb{R}$ とする。各 $x\in[a,b]$ を固定したとき 1 変数関数 $y\mapsto f(x,y)$ が $[c,d]$ 上可積分で、そうして定まる関数 $$ F(x) = \int_c^d f(x,y)\,dy $$ が $[a,b]$ 上可積分であるとき、 $$ \int_a^b\!\left(\int_c^d f(x,y)\,dy\right)dx = \int_a^b F(x)\,dx $$ を $f$ の**累次積分**（反復積分）という。$x$ と $y$ の役割を入れ替えた $\displaystyle\int_c^d\!\left(\int_a^b f(x,y)\,dx\right)dy$ も同様に定める。 2 つの積分の性格の違いを並べておきます。 | | 重積分 $\displaystyle\iint_R f\,dx\,dy$ | 累次積分 $\displaystyle\int_a^b\!\left(\int_c^d f\,dy\right)dx$ | |---|---|---| | 定義 | 平面の分割によるリーマン和の極限 | 1 変数の定積分を 2 回行う操作 | | 幾何的意味 | グラフの下の符号付き体積 | 断面積 $F(x)$ を $x$ 方向に積み上げたもの | | 存在条件 | 有界性＋ | 各断面の可積分性＋ $F$ の可積分性 | | 計算 | ふつう直接は計算できない | 1 変数の技法がそのまま使える | この 2 つを結ぶのがフビニの定理です。仮定 (ii) は「切り口ごとの面積がそもそも定義できる」ことの要請で、これを落とすと右辺の意味が失われます。 $R=[a,b]\times[c,d]$ とし、$f\colon R\to\mathbb{R}$ が次の 2 条件を満たすとする。 1. $f$ は $R$ 上有界かつ可積分である。 2. 各 $x\in[a,b]$ に対し、1 変数関数 $y\mapsto f(x,y)$ は $[c,d]$ 上可積分である。このとき $F(x)=\displaystyle\int_c^d f(x,y)\,dy$ は $[a,b]$ 上可積分で、 $$ \int_a^b\!\left(\int_c^d f(x,y)\,dy\right)dx = \iint_R f(x,y)\,dx\,dy $$ が成り立つ。$x$ と $y$ の役割を入れ替えた主張（各 $y$ について $x\mapsto f(x,y)$ が可積分なら $\int_c^d\!\left(\int_a^b f\,dx\right)dy = \iint_R f$）も同様に成り立つ。 $R$ の分割 $\Delta$ を任意に取り、その $x$ 方向の分点だけを取り出した $[a,b]$ の分割を $\Delta_x$ と書きます。 $x\in[x_{i-1},x_i]$ を固定します。仮定 2 と 1 変数の積分の区間加法性から $$ F(x) = \int_c^d f(x,y)\,dy = \sum_{j=1}^{n}\int_{y_{j-1}}^{y_j} f(x,y)\,dy . $$ $y\in[y_{j-1},y_j]$ のとき $(x,y)\in R_{ij}$ なので $m_{ij}\le f(x,y)\le M_{ij}$ であり、1 変数の積分の単調性から $$ m_{ij}\,\Delta y_j \le \int_{y_{j-1}}^{y_j} f(x,y)\,dy \le M_{ij}\,\Delta y_j . $$ $j$ について足すと、$[x_{i-1},x_i]$ のすべての $x$ について $$ \sum_{j=1}^{n} m_{ij}\Delta y_j \ \le\ F(x)\ \le\ \sum_{j=1}^{n} M_{ij}\Delta y_j $$ が成り立ちます。左辺・右辺は $x$ によらない定数なので、$[x_{i-1},x_i]$ 上での $F$ の下限・上限について $$ \inf_{[x_{i-1},x_i]} F \ \ge\ \sum_{j} m_{ij}\Delta y_j, \qquad \sup_{[x_{i-1},x_i]} F \ \le\ \sum_{j} M_{ij}\Delta y_j $$ が従います。$\Delta x_i>0$ を掛けて $i$ について足せば $$ s(f,\Delta) = \sum_{i,j} m_{ij}\Delta y_j\Delta x_i \ \le\ s(F,\Delta_x) \ \le\ S(F,\Delta_x)\ \le\ \sum_{i,j} M_{ij}\Delta y_j \Delta x_i = S(f,\Delta) $$ を得ます（真ん中の不等号は不足和 $\le$ 過剰和）。ここで 1 変数の下積分・上積分は $s(F,\Delta_x)\le \underline{\int_a^b} F \le \overline{\int_a^b} F\le S(F,\Delta_x)$ を満たすので、結局すべての分割 $\Delta$ に対して $$ s(f,\Delta)\ \le\ \underline{\int_a^b} F \ \le\ \overline{\int_a^b} F\ \le\ S(f,\Delta) . $$ 仮定 1 とより $\sup_\Delta s(f,\Delta)=\inf_\Delta S(f,\Delta)=\iint_R f$ ですから、$\Delta$ について上限・下限を取ると挟み撃ちで $$ \underline{\int_a^b} F = \overline{\int_a^b} F = \iint_R f(x,y)\,dx\,dy . $$ 上積分と下積分が一致するので $F$ は $[a,b]$ 上可積分であり（1 変数版の）、その積分値は $\iint_R f$ です。$x$ と $y$ を入れ替えても議論はまったく同じです。 $f$ が $R=[a,b]\times[c,d]$ 上連続ならば、2 つの累次積分はともに存在して等しく、重積分に一致する。すなわち $$ \int_a^b\!\left(\int_c^d f(x,y)\,dy\right)dx = \int_c^d\!\left(\int_a^b f(x,y)\,dx\right)dy = \iint_R f(x,y)\,dx\,dy . $$ $f$ は有界閉集合 $R$ 上の連続関数なので有界で、より可積分です。これでの仮定 1 が満たされます。次に $x\in[a,b]$ を固定すると、$y\mapsto f(x,y)$ は $[c,d]$ 上の連続関数なので 1 変数の意味で可積分です。これで仮定 2 も満たされ、から第 1 の式と第 3 の式が等しいことがわかります。同じ議論を $x$ と $y$ を入れ替えて行えば、第 2 の式も第 3 の式に等しくなります。 $$ \iint_{[0,1]\times[0,\pi]} x\cos(xy)\,dx\,dy $$ を計算します。被積分関数は連続なのでよりどちらの順でもよく、値は一致します。 **先に $y$ で積分する場合。** $x$ を固定すると、$y$ の関数として $\cos(xy)$ の原始関数は $\sin(xy)/x$（$x\neq0$）ですから $$ \int_0^{\pi} x\cos(xy)\,dy = x\cdot\left[\frac{\sin(xy)}{x}\right]_{y=0}^{y=\pi} = \sin(\pi x). $$ $x=0$ のときは被積分関数が $0$ なので値は $0$ で、これも $\sin(\pi\cdot 0)=0$ と一致します。よって $$ \int_0^1 \sin(\pi x)\,dx = \left[-\frac{\cos(\pi x)}{\pi}\right]_0^1 = \frac{1+1}{\pi} = \frac{2}{\pi}. $$ **先に $x$ で積分する場合。** $y>0$ を固定して部分積分（$u=x$, $dv=\cos(xy)\,dx$, $v=\sin(xy)/y$）を行うと $$ \int_0^1 x\cos(xy)\,dx = \left[\frac{x\sin(xy)}{y}\right]_0^1 - \int_0^1\frac{\sin(xy)}{y}\,dx = \frac{\sin y}{y} - \frac{1}{y}\left[-\frac{\cos(xy)}{y}\right]_0^1 = \frac{\sin y}{y}+\frac{\cos y - 1}{y^2}. $$ これを $G(y)$ と書きます。$y\to 0$ では $\sin y/y\to1$、$(\cos y-1)/y^2\to-1/2$ なので $G(0^+)=1/2$ となり、$y=0$ での直接計算 $\int_0^1 x\,dx = 1/2$ と一致します。次に $$ \frac{d}{dy}\left(\frac{1-\cos y}{y}\right) = \frac{y\sin y - (1-\cos y)}{y^2} = \frac{\sin y}{y} - \frac{1-\cos y}{y^2} = G(y) $$ なので、$(1-\cos y)/y$ が $G$ の原始関数です（$y\to0$ での値は $0$）。したがって $$ \int_0^{\pi} G(y)\,dy = \left[\frac{1-\cos y}{y}\right]_{0}^{\pi} = \frac{2}{\pi} - 0 = \frac{2}{\pi}. $$ 両者は確かに一致しました。ただし手間はまるで違います。順序を選べること自体がフビニの定理の実用的な価値です。 ## 5. 一般の領域上の重積分と積分順序現実に出てくる領域は長方形ではなく、三角形・円板・2 曲線に挟まれた領域などです。これらを扱う標準的な手口は「長方形まで $0$ で延長する」ことです。 $D\subset\mathbb{R}^2$ を有界集合、$f\colon D\to\mathbb{R}$ を有界関数とする。$D\subset R$ となる閉長方形 $R$ を 1 つ取り、 $$ \tilde f(x,y) = \begin{cases} f(x,y) & ((x,y)\in D)\\ 0 & ((x,y)\in R\setminus D)\end{cases} $$ とおく。$\tilde f$ が $R$ 上可積分であるとき、$f$ は $D$ 上可積分であるといい $$ \iint_D f(x,y)\,dx\,dy := \iint_R \tilde f(x,y)\,dx\,dy $$ と定める。とくに $f\equiv 1$ が $D$ 上可積分であるとき $D$ は**面積確定**であるといい、$\lvert D\rvert = \iint_D 1\,dx\,dy$ を $D$ の**面積**という。この定義が $R$ の取り方によらないことは次のようにわかります。$R\subset R'$ なる長方形を取り、$R'$ の分割で $R$ の 4 辺を分点に含むものを考えると、$R'\setminus R$ の側の小長方形では $\tilde f\equiv 0$ なのでリーマン和に寄与せず、$R$ 側の和とちょうど一致します。一般の 2 つの長方形は、両方を含む長方形を経由して比べればよいです。 $\varphi\colon[a,b]\to\mathbb{R}$ を連続関数とすると、そのグラフ $G=\{(t,\varphi(t)) \mid t\in[a,b]\}$ は面積確定で $\lvert G\rvert = 0$ である。同じことは横向きのグラフ $\{(\psi(t),t)\}$ についても成り立つ。 $\varepsilon>0$ を取ります。$\varphi$ は有界閉区間上の連続関数なので一様連続で、$\lvert t-t'\rvert<\delta \Rightarrow \lvert\varphi(t)-\varphi(t')\rvert<\varepsilon$ となる $\delta>0$ があります。$[a,b]$ を $n$ 等分して $(b-a)/n<\delta$ とし、小区間を $I_1,\dots,I_n$ とします。$\varphi$ は各 $I_k$ 上で最大値 $\beta_k$ と最小値 $\alpha_k$ を取り、$I_k$ 内の 2 点での値の差は $\varepsilon$ 未満なので $\beta_k-\alpha_k\le\varepsilon$ です。 $G$ を含む閉長方形 $R=[a,b]\times[c,d]$ を取り、$G$ 上で $1$、$R\setminus G$ で $0$ となる関数を $g$ とします。$R$ の分割 $\Delta$ として、$x$ 方向は上の $n$ 等分、$y$ 方向は $c,d$ と $2n$ 個の値 $\alpha_1,\beta_1,\dots,\alpha_n,\beta_n$ を分点に持つものを取ります。$I_k$ の上にある $G$ の部分は帯 $I_k\times[\alpha_k,\beta_k]$ に含まれ、$\alpha_k,\beta_k$ は分点なので、この帯は第 $k$ 列の小長方形いくつかの合併になっています。それ以外の小長方形では $g\equiv0$、すなわち $\sup g=0$ ですから $$ S(g,\Delta) \le \sum_{k=1}^{n} \frac{b-a}{n}\,(\beta_k-\alpha_k) \le \sum_{k=1}^{n}\frac{b-a}{n}\,\varepsilon = (b-a)\varepsilon . $$ 一方、どの小長方形も幅と高さが正であり、各 $x$ 座標に対して $G$ の点は高々 1 個しかないので、小長方形が $G$ に含まれることはなく、必ず $G$ の外の点を含みます。そこでは $g=0$ なので $\inf g=0$ となり、$s(g,\Delta)=0$ です。したがって $S-s\le(b-a)\varepsilon$ が任意の $\varepsilon>0$ で成り立ち、より $g$ は可積分で、積分値は $\inf_\Delta S=0$ です。横向きのグラフも $x$ と $y$ を入れ替えれば同じです。以下では、証明を与えずに次の事実を使います。 1. **リーマンの可積分条件（ジョルダン版）**：有界関数 $g\colon R\to\mathbb{R}$ の不連続点全体が面積 $0$ の集合に含まれるならば、$g$ は $R$ 上可積分である。 2. 面積 $0$ の集合の有限個の合併はまた面積 $0$ であり、線分は面積 $0$ である。 2 は定義から直ちに従います（線分は高さをいくらでも小さくした長方形で覆えます）。1 の証明は杉浦光夫『解析入門 II』の重積分の章、あるいは W. Rudin, *Principles of Mathematical Analysis* にあります。1 は「不連続点が少なければ可積分」という直観の正確な形で、境界の上でだけ不連続になる $\tilde f$ を扱うのにちょうど必要な道具です。 $\varphi_1,\varphi_2\colon[a,b]\to\mathbb{R}$ を連続関数で $\varphi_1\le\varphi_2$ を満たすものとし、 $$ D = \{(x,y)\in\mathbb{R}^2 \mid a\le x\le b,\ \varphi_1(x)\le y\le \varphi_2(x)\} $$ とおく（このような $D$ を**縦線集合**という）。$f\colon D\to\mathbb{R}$ が連続ならば、$f$ は $D$ 上可積分で $$ \iint_D f(x,y)\,dx\,dy = \int_a^b\!\left(\int_{\varphi_1(x)}^{\varphi_2(x)} f(x,y)\,dy\right)dx $$ が成り立つ。同様に、$\psi_1\le\psi_2$ を連続関数として $E=\{(x,y)\mid c\le y\le d,\ \psi_1(y)\le x\le\psi_2(y)\}$（**横線集合**）と $E$ 上連続な $f$ に対しては $$ \iint_E f(x,y)\,dx\,dy = \int_c^d\!\left(\int_{\psi_1(y)}^{\psi_2(y)} f(x,y)\,dx\right)dy . $$ $\varphi_1,\varphi_2$ は連続だから有界で、$D$ は有界です。さらに $D$ は有界閉集合なので、その上の連続関数 $f$ は有界です。そこで $D\subset R=[a,b]\times[c,d]$ となる閉長方形を取り、の $\tilde f$ を考えます。 **（仮定 1 の確認）** $\tilde f$ の不連続点を調べます。$D$ の内点では、その近傍で $\tilde f=f$ であり $f$ は連続です。$D$ の外部（補集合の内点）では近傍で $\tilde f\equiv0$ なので連続です。したがって不連続点は境界 $\partial D$ に含まれます。$\partial D$ は 2 つのグラフ $y=\varphi_1(x)$、$y=\varphi_2(x)$ と、2 本の縦線分 $x=a$、$x=b$ の部分集合の合併に含まれます。との 2 より $\partial D$ は面積 $0$、よっての 1 より $\tilde f$ は $R$ 上可積分です。これでの仮定 1 が満たされました。 **（仮定 2 の確認）** $x\in[a,b]$ を固定します。$y\mapsto\tilde f(x,y)$ は、$[c,\varphi_1(x))$ と $(\varphi_2(x),d]$ では $0$、$[\varphi_1(x),\varphi_2(x)]$ では $y\mapsto f(x,y)$ に等しく、$f$ が $D$ 上連続なのでこの区間上連続です。つまり有界で、不連続点は高々 2 点 $y=\varphi_1(x),\varphi_2(x)$ しかありません。1 変数の有界関数は不連続点が有限個なら可積分ですから、仮定 2 も満たされ、さらに区間加法性から $$ \int_c^d \tilde f(x,y)\,dy = \int_{\varphi_1(x)}^{\varphi_2(x)} f(x,y)\,dy $$ です（$\varphi_1(x)=\varphi_2(x)$ のときは両辺 $0$）。以上よりが適用でき、 $$ \iint_D f = \iint_R\tilde f = \int_a^b\!\left(\int_c^d \tilde f(x,y)\,dy\right)dx = \int_a^b\!\left(\int_{\varphi_1(x)}^{\varphi_2(x)} f(x,y)\,dy\right)dx $$ を得ます。横線集合の場合は $x$ と $y$ の役割を入れ替えれば同じ議論です。 $$ \int_0^1\!\left(\int_y^1 e^{x^2}\,dx\right)dy $$ を計算します。内側の $\int_y^1 e^{x^2}dx$ は、$e^{x^2}$ の原始関数が初等関数で書けないため、この順序のままでは手が出ません。積分領域を読み取ります。外側が $0\le y\le 1$、内側が $y\le x\le 1$ ですから、 $$ D = \{(x,y)\mid 0\le y\le 1,\ y\le x\le 1\} $$ で、これは $(0,0),(1,0),(1,1)$ を頂点とする三角形です（$0\le y\le x\le 1$ と書くと見やすいです）。同じ三角形を縦線集合として書き直すと、各 $x\in[0,1]$ に対して $y$ は $0$ から $x$ まで動くので $$ D = \{(x,y)\mid 0\le x\le 1,\ 0\le y\le x\}. $$ 被積分関数 $e^{x^2}$ は $D$ 上連続なのでが両方の表示に使え、どちらの累次積分も $\iint_D e^{x^2}dx\,dy$ に等しくなります。したがって $$ \int_0^1\!\left(\int_y^1 e^{x^2}dx\right)dy = \int_0^1\!\left(\int_0^x e^{x^2}\,dy\right)dx = \int_0^1 e^{x^2}\bigl[y\bigr]_{0}^{x}\,dx = \int_0^1 x e^{x^2}\,dx . $$ 最後は $t=x^2$ と置換して（$dt=2x\,dx$） $$ \int_0^1 xe^{x^2}dx = \left[\frac{e^{x^2}}{2}\right]_0^1 = \frac{e-1}{2}. $$ 内側の積分に現れた $x$ という因子は、三角形の縦の長さ（断面の長さ）です。これが $xe^{x^2}$ の $x$ の正体で、順序交換の御利益はここに集約されています。 $$ f(x,y) = \frac{x^2-y^2}{(x^2+y^2)^2} $$ を $(0,1]\times(0,1]$ 上で考えます。この $f$ は原点に近づけると非有界です。実際 $y$ を固定して $x\to 0^+$ とすると $f(x,y)\to -1/y^2$ ですが、$y=x$ 上では $f=0$、$y\to0^+$ を先にすると $f(x,y)\to 1/x^2$ で、$f(x,x/2)= \frac{x^2-x^2/4}{(5x^2/4)^2}=\frac{12}{25x^2}\to+\infty$ です。したがって $f$ は $[0,1]^2$ 上有界でなく、の意味では**そもそも重積分が定義されません**。それでも 2 つの累次積分は計算できます。$x>0$ を固定すると $y\mapsto f(x,y)$ は $[0,1]$ 上連続です（分母は $x^2$ 以上）。商の微分から $$ \frac{\partial}{\partial y}\left(\frac{y}{x^2+y^2}\right) = \frac{(x^2+y^2)-y\cdot 2y}{(x^2+y^2)^2} = \frac{x^2-y^2}{(x^2+y^2)^2} = f(x,y) $$ なので $$ \int_0^1 f(x,y)\,dy = \left[\frac{y}{x^2+y^2}\right]_{y=0}^{y=1} = \frac{1}{x^2+1}. $$ これを $x$ で積分して $$ \int_0^1\!\left(\int_0^1 f(x,y)\,dy\right)dx = \int_0^1\frac{dx}{1+x^2} = \Bigl[\arctan x\Bigr]_0^1 = \frac{\pi}{4}. $$ 一方 $f(y,x) = \dfrac{y^2-x^2}{(x^2+y^2)^2} = -f(x,y)$ という反対称性があるので、$x$ と $y$ の名前を入れ替えて同じ計算をすれば $$ \int_0^1\!\left(\int_0^1 f(x,y)\,dx\right)dy = -\frac{\pi}{4}. $$ （$x=0$ では内側の積分 $\int_0^1 f(0,y)\,dy=\int_0^1(-y^{-2})dy$ が発散するので、外側の積分は $x\in(0,1]$ 上の広義積分として読みます。$1/(1+x^2)$ は $x=0$ まで連続に延びるので値は $\pi/4$ のままです。） $\pi/4 \ne -\pi/4$ ですから、2 つの累次積分は食い違います。これはと矛盾しません。仮定 1（$R$ 上可積分）が破れているからです。**フビニの定理は「重積分が存在すること」を前提に、それを累次積分に翻訳する定理であって、2 つの累次積分が勝手に一致することを主張するものではありません。** 積分を和に置き換えても同じ現象が現れます。$m,n\ge1$ に対し $$ a_{mn} = \begin{cases} 1 & (m=n)\\ -1 & (m=n+1)\\ 0 & (\text{その他})\end{cases} $$ とおきます。$m$ を固定して $n$ について足すと、$m=1$ のときだけ $a_{11}=1$ で他は $0$ だから和は $1$、$m\ge2$ のときは $a_{mm}=1$ と $a_{m,m-1}=-1$ が打ち消して $0$ です。よって $\sum_m\bigl(\sum_n a_{mn}\bigr)=1$。逆に $n$ を固定して $m$ について足すと、$a_{nn}=1$ と $a_{n+1,n}=-1$ が打ち消して毎回 $0$ なので $\sum_n\bigl(\sum_m a_{mn}\bigr)=0$。和の順序で $1$ と $0$ に分かれます。原因はと同じで、$\sum_{m,n}\lvert a_{mn}\rvert=\infty$（絶対収束しない）ことです。級数の絶対収束と順序交換については [級数と収束判定](/mathematics/calculus/series-and-convergence) のを参照してください。逆に $f\ge0$（または $\lvert f\rvert$ の累次積分が有限）なら順序交換してよい、というのがトネリの定理で、無限区間や非有界関数を扱うにはルベーグ積分の枠組みが必要になります。 ## 6. 変数変換とヤコビアン 1 変数のは、$g$ が $C^1$ 級で単調なら $$ \int_{g(\alpha)}^{g(\beta)} f(x)\,dx = \int_{\alpha}^{\beta} f(g(t))\,g'(t)\,dt $$ でした。向きを気にせず「長さ」だけで読むと、$g$ は微小区間 $[t,t+dt]$ を長さ $\lvert g'(t)\rvert\,dt$ の区間に写すので、$dx = \lvert g'(t)\rvert\,dt$ という置き換えが起きています。2 変数では $dx\,dy$ が面積なので、「微小な面積を何倍にするか」を表す量が必要になります。それがヤコビアンです。 ### 6.1. まず一次変換で考える $A=\begin{pmatrix} a_{11} & a_{12}\\ a_{21} & a_{22}\end{pmatrix}$ を実 2 次正方行列とし、$T(\boldsymbol u)=A\boldsymbol u$ とする。 1. 単位正方形 $[0,1]\times[0,1]$ の像 $T([0,1]^2)$ は $A$ の 2 つの列ベクトル $\boldsymbol a_1=(a_{11},a_{21})$、$\boldsymbol a_2=(a_{12},a_{22})$ の張る平行四辺形であり、その面積は $\lvert\det A\rvert$ に等しい。 2. さらに $\det A\ne0$ とすると、面積確定な有界集合 $E\subset\mathbb{R}^2$ に対し $T(E)$ も面積確定で $\lvert T(E)\rvert = \lvert\det A\rvert\cdot\lvert E\rvert$ が成り立つ。 **(1)** $\boldsymbol u=(u,v)\in[0,1]^2$ に対し $T(\boldsymbol u)=u\boldsymbol a_1+v\boldsymbol a_2$ ですから、像は $\boldsymbol a_1,\boldsymbol a_2$ が張る平行四辺形そのものです。その面積は「底辺 × 高さ」で $$ \lVert\boldsymbol a_1\rVert\cdot\lVert\boldsymbol a_2\rVert\sin\theta $$ （$\theta\in[0,\pi]$ は 2 つのベクトルのなす角、$\sin\theta\ge0$）。$\cos\theta = \dfrac{\langle \boldsymbol a_1,\boldsymbol a_2\rangle}{\lVert\boldsymbol a_1\rVert\lVert\boldsymbol a_2\rVert}$ を使って 2 乗すると $$ \bigl(\lVert\boldsymbol a_1\rVert\lVert\boldsymbol a_2\rVert\sin\theta\bigr)^2 = \lVert\boldsymbol a_1\rVert^2\lVert\boldsymbol a_2\rVert^2 - \langle\boldsymbol a_1,\boldsymbol a_2\rangle^2 . $$ 右辺を成分で書き下すと $$ \begin{aligned} &\ (a_{11}^2+a_{21}^2)(a_{12}^2+a_{22}^2)-(a_{11}a_{12}+a_{21}a_{22})^2\\ =&\ a_{11}^2a_{12}^2+a_{11}^2a_{22}^2+a_{21}^2a_{12}^2+a_{21}^2a_{22}^2 -a_{11}^2a_{12}^2-2a_{11}a_{12}a_{21}a_{22}-a_{21}^2a_{22}^2\\ =&\ a_{11}^2a_{22}^2-2a_{11}a_{22}a_{12}a_{21}+a_{12}^2a_{21}^2 =(a_{11}a_{22}-a_{12}a_{21})^2 = (\det A)^2 . \end{aligned} $$ 両辺の正の平方根を取って、面積は $\lvert\det A\rvert$ です。 **(2)** 辺の長さ $h$ の正方形 $Q$（左下の頂点を $\boldsymbol p$）に対し、$T$ の線形性から $T(Q)=T(\boldsymbol p)+h\,T([0,1]^2)$ です。面積は平行移動で変わらず、$h$ 倍の相似拡大で $h^2$ 倍になるので、(1) より $$ \lvert T(Q)\rvert = h^2\lvert\det A\rvert = \lvert\det A\rvert\cdot\lvert Q\rvert . $$ 次に、平面を一辺 $1/N$ の正方形の格子に分け、$E$ に完全に含まれる正方形の合併を $E_N^{-}$、$E$ と交わる正方形の合併を $E_N^{+}$ とします。$E_N^{-}\subset E\subset E_N^{+}$ であり、$\lvert E_N^{-}\rvert$ と $\lvert E_N^{+}\rvert$ はそれぞれ「$E$ 上で $1$、外で $0$」とした関数の不足和と過剰和にほかならないので、$E$ が面積確定であることから $N\to\infty$ でともに $\lvert E\rvert$ に収束します。$\det A\ne0$ より $T$ は全単射なので $T(E_N^{-})\subset T(E)\subset T(E_N^{+})$ であり、内部が交わらない有限個の平行四辺形の面積の和として $$ \lvert T(E_N^{-})\rvert = \lvert\det A\rvert\,\lvert E_N^{-}\rvert, \qquad \lvert T(E_N^{+})\rvert = \lvert\det A\rvert\,\lvert E_N^{+}\rvert $$ です（面積の有限加法性と、により平行四辺形の辺は面積 $0$ であることを使いました）。面積の単調性で挟み撃ちにすれば、$T(E)$ は面積確定で $\lvert T(E)\rvert=\lvert\det A\rvert\,\lvert E\rvert$ を得ます。 ### 6.2. ヤコビアン：局所的な面積の拡大率 $U\subset\mathbb{R}^2$ を開集合、$\Phi\colon U\to\mathbb{R}^2$、$\Phi(u,v)=\bigl(x(u,v),\,y(u,v)\bigr)$ を $C^1$ 級写像とする。 $$ J_\Phi(u,v) = \begin{pmatrix} \dfrac{\partial x}{\partial u} & \dfrac{\partial x}{\partial v}\\[10pt] \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v} \end{pmatrix} $$ を $\Phi$ の**ヤコビ行列**、その行列式 $$ \det J_\Phi(u,v) = \frac{\partial x}{\partial u}\frac{\partial y}{\partial v}-\frac{\partial x}{\partial v}\frac{\partial y}{\partial u} \ \ \left(=\frac{\partial(x,y)}{\partial(u,v)} \text{ とも書く}\right) $$ を**ヤコビアン**（関数行列式）という。ヤコビ行列は $\Phi$ の全微分（1 次近似）を表す行列でした。[多変数関数の微分と偏微分](/mathematics/calculus/multivariable-differentiation) ので見たとおり、$C^1$ 級写像は各点で $$ \Phi(u+\Delta u,\ v+\Delta v) = \Phi(u,v) + J_\Phi(u,v)\begin{pmatrix}\Delta u\\ \Delta v\end{pmatrix} + o\bigl(\lVert(\Delta u,\Delta v)\rVert\bigr) $$ と 1 次近似できます。そこで $uv$ 平面の微小な長方形 $[u,u+\Delta u]\times[v,v+\Delta v]$ を $\Phi$ で写すと、誤差項を無視すれば、$\Phi(u,v)$ を頂点として 2 つのベクトル $$ \frac{\partial\Phi}{\partial u}\Delta u,\qquad \frac{\partial\Phi}{\partial v}\Delta v $$ が張る平行四辺形になります。この 2 ベクトルを列とする行列は $J_\Phi\,\mathrm{diag}(\Delta u,\Delta v)$ ですから、よりその面積は $$ \lvert\det J_\Phi(u,v)\rvert\,\Delta u\,\Delta v $$ です。「微小面積 $\Delta u\Delta v$ がヤコビアンの絶対値倍に化ける」——これが変数変換公式の正体です。

### 6.3. 変数変換公式 $U,V\subset\mathbb{R}^2$ を開集合、$\Phi\colon U\to V$ を全単射な $C^1$ 級写像で、すべての $(u,v)\in U$ において $\det J_\Phi(u,v)\ne0$ を満たすものとする。$E$ を有界閉かつ面積確定な集合で $E\subset U$ を満たすものとすると、$\Phi(E)$ も有界閉かつ面積確定であり、$\Phi(E)$ 上の任意の連続関数 $f$ に対して $$ \iint_{\Phi(E)} f(x,y)\,dx\,dy = \iint_{E} f\bigl(\Phi(u,v)\bigr)\,\bigl\lvert\det J_\Phi(u,v)\bigr\rvert\,du\,dv $$ が成り立つ。実用上は、$\Phi$ の単射性や $\det J_\Phi\ne0$ が $E$ の中の**面積 $0$ の集合の上で崩れていても**公式はそのまま成り立ちます。次のの極座標変換では、$r=0$ の辺（$\theta$ 軸上の線分）が 1 点 $(0,0)$ に潰れて単射性が崩れ、そこでは $\det J_\Phi = r = 0$ ですが、この辺は面積 $0$ なので積分に寄与せず、結論は変わりません。直観的には、$E$ から面積 $0$ の閉集合を除いた開集合に定理を適用し、除いた部分の寄与が $0$ であることを確かめればよい、ということです。厳密な形と証明は杉浦光夫『解析入門 II』の変数変換の節にあります。の完全な証明はこの記事には収めきれないので、標準的な文献に譲ります（杉浦光夫『解析入門 II』、W. Rudin, *Principles of Mathematical Analysis* 第 10 章、M. Spivak, *Calculus on Manifolds* 第 3 章）。方針だけ述べておきます。 1. 一次変換に対しては、公式はそのものです（$f\equiv1$ の場合が面積の関係式、一般の $f$ はリーマン和を移し替えれば従います）。 2. 一般の $\Phi$ は、各点の近くで「$u$ だけを動かす写像」と「$v$ だけを動かす写像」の合成に分解できます（Rudin の primitive mapping）。1 変数の置換積分とを組み合わせると、この種の写像については公式が示せます。 3. 微分可能性から、小さな正方形の上では $\Phi$ は一次変換 $J_\Phi$ で近似でき、誤差は $\lvert\Delta\rvert$ について高位の無限小です。逆関数定理でこの近似を制御します。 4. コンパクト性（有限被覆）で局所的な結果を貼り合わせます。 $$ \Phi(r,\theta) = (r\cos\theta,\ r\sin\theta) $$ とします。$x=r\cos\theta$、$y=r\sin\theta$ の偏微分を並べると $$ J_\Phi = \begin{pmatrix}\cos\theta & -r\sin\theta\\ \sin\theta & r\cos\theta\end{pmatrix}, \qquad \det J_\Phi = r\cos^2\theta + r\sin^2\theta = r . $$ $\Phi$ は開集合 $U=\{(r,\theta)\mid r>0,\ 0<\theta<2\pi\}$ の上で単射で $\det J_\Phi=r>0$、像は平面から半直線 $\{(x,0)\mid x\ge0\}$ を除いた開集合です。半径 $a$ の円板 $D_a=\{(x,y)\mid x^2+y^2\le a^2\}$ に対しては $E=[0,a]\times[0,2\pi]$ と取れば $\Phi(E)=D_a$ で、単射性が崩れるのは $r=0$ の辺と $\theta=0,2\pi$ の 2 辺だけ、いずれも面積 $0$ です。より $$ \iint_{D_a} f(x,y)\,dx\,dy = \int_0^{2\pi}\!\!\int_0^{a} f(r\cos\theta,\ r\sin\theta)\;r\,dr\,d\theta $$ が使えます（右辺は長方形上の累次積分なのでで順序も自由です）。 **確認 1（円の面積）。** $f\equiv1$ とすると $$ \int_0^{2\pi}\!\!\int_0^a r\,dr\,d\theta = \int_0^{2\pi}\frac{a^2}{2}\,d\theta = \pi a^2 . $$ **確認 2（半球の体積）。** $f(x,y)=\sqrt{a^2-x^2-y^2}$ とすると $f(\Phi(r,\theta))=\sqrt{a^2-r^2}$ で $$ \int_0^{2\pi}\!\!\int_0^{a}\sqrt{a^2-r^2}\;r\,dr\,d\theta = 2\pi\left[-\frac{(a^2-r^2)^{3/2}}{3}\right]_0^{a} = 2\pi\cdot\frac{a^3}{3} = \frac{2}{3}\pi a^3 . $$ これは半径 $a$ の球の体積 $\frac43\pi a^3$ のちょうど半分で、正しい値です。因子 $r$ を忘れるとどうなるかも見ておきます。$\int_0^{2\pi}\int_0^a\sqrt{a^2-r^2}\,dr\,d\theta = 2\pi\cdot\frac{\pi a^2}{4}$ となり、半球の体積とは似ても似つきません。$r$ は「外側の $\theta$ 方向の弧が半径に比例して長くなる」ことを表す、なくてはならない因子です。 $$ I = \int_{-\infty}^{\infty} e^{-x^2}\,dx $$ を求めます。まず収束を確認します。$\lvert x\rvert\ge1$ では $x^2\ge \lvert x\rvert$ なので $e^{-x^2}\le e^{-\lvert x\rvert}$ であり、$\int_1^{\infty}e^{-x}dx = e^{-1}<\infty$ ですから広義積分は収束します。 $I_R=\int_{-R}^{R}e^{-x^2}dx$ とおきます。$e^{-(x^2+y^2)}=e^{-x^2}e^{-y^2}$ は連続なのでが使え、正方形 $S_R=[-R,R]^2$ 上で $$ \iint_{S_R} e^{-(x^2+y^2)}dx\,dy = \int_{-R}^{R}\!\left(\int_{-R}^{R} e^{-x^2}e^{-y^2}\,dx\right)dy = \int_{-R}^{R} e^{-y^2}I_R\,dy = I_R^{\,2}. $$ 次に円板で挟みます。$D_R\subset S_R\subset D_{\sqrt2 R}$（正方形の外接円の半径は $\sqrt2R$）であり、被積分関数は正なので、の加法性と単調性から $$ \iint_{D_R} e^{-(x^2+y^2)}\,dx\,dy \ \le\ I_R^{\,2}\ \le\ \iint_{D_{\sqrt2 R}} e^{-(x^2+y^2)}\,dx\,dy . $$ 円板上の積分はの極座標公式で計算できます。$t=r^2$ と置換して（$dt=2r\,dr$） $$ \iint_{D_a} e^{-(x^2+y^2)}dx\,dy = \int_0^{2\pi}\!\!\int_0^{a} e^{-r^2}\,r\,dr\,d\theta = 2\pi\left[-\frac{e^{-r^2}}{2}\right]_0^{a} = \pi\bigl(1-e^{-a^2}\bigr). $$ したがって $$ \pi\bigl(1-e^{-R^2}\bigr)\ \le\ I_R^{\,2}\ \le\ \pi\bigl(1-e^{-2R^2}\bigr). $$ $R\to\infty$ とすると両端はともに $\pi$ に収束し、$I_R\to I$ なので、はさみうちの原理から $I^2=\pi$。$e^{-x^2}>0$ より $I>0$ ですから $$ \int_{-\infty}^{\infty}e^{-x^2}\,dx = \sqrt{\pi}. $$ 1 変数のままでは $e^{-x^2}$ の原始関数が初等関数で書けないため手が出ませんが、2 乗して 2 次元にし、極座標に移ると $r\,dr$ の $r$ が $e^{-r^2}$ の微分の因子をちょうど供給してくれます。ヤコビアンが計算を可能にした、いちばん有名な例です。 ## 7. 演習 $\displaystyle\iint_{[0,1]\times[0,2]}(x+y^2)\,dx\,dy$ を、2 通りの累次積分で計算し、値が一致することを確かめてください。被積分関数 $x+y^2$ は長方形 $R=[0,1]\times[0,2]$ 上で連続なので、よりどちらの順序でもよく、値は重積分に一致します。先に $y$ で積分すると $$ \int_0^1\!\left(\int_0^2 (x+y^2)\,dy\right)dx = \int_0^1\left[xy+\frac{y^3}{3}\right]_{y=0}^{y=2}dx = \int_0^1\left(2x+\frac83\right)dx = \bigl[x^2\bigr]_0^1+\frac83 = 1+\frac83 = \frac{11}{3}. $$ 先に $x$ で積分すると $$ \int_0^2\!\left(\int_0^1 (x+y^2)\,dx\right)dy = \int_0^2\left[\frac{x^2}{2}+xy^2\right]_{x=0}^{x=1}dy = \int_0^2\left(\frac12+y^2\right)dy = 1+\frac83 = \frac{11}{3}. $$ どちらも $\dfrac{11}{3}$ です。 $$ \int_0^1\!\left(\int_y^1 \frac{\sin x}{x}\,dx\right)dy $$ の値を求めてください。ただし $x=0$ では被積分関数の値を $1$ と定めます（$\lim_{x\to0}\sin x/x=1$ なので、こうすると $[0,1]$ 上連続になります）。 $g(x)=\dfrac{\sin x}{x}$（$x\ne0$）、$g(0)=1$ とおくと $g$ は $[0,1]$ 上連続です。$\int_y^1 g(x)\,dx$ の原始関数は初等関数で書けないので、この順序では計算できません。積分順序を交換します。積分領域は、外側が $0\le y\le 1$、内側が $y\le x\le 1$ ですから $$ D = \{(x,y)\mid 0\le y\le 1,\ y\le x\le 1\} = \{(x,y)\mid 0\le y\le x\le 1\} $$ で、$(0,0),(1,0),(1,1)$ を頂点とする三角形です。これを縦線集合として書き直すと、各 $x\in[0,1]$ に対し $y$ は $0$ から $x$ まで動きます。$g(x)$ は $D$ 上連続なのでを両方の表示に適用でき、 $$ \int_0^1\!\left(\int_y^1 g(x)\,dx\right)dy = \iint_D g(x)\,dx\,dy = \int_0^1\!\left(\int_0^x g(x)\,dy\right)dx = \int_0^1 x\,g(x)\,dx . $$ $x\,g(x)=\sin x$（$x=0$ でも両辺 $0$）なので $$ \int_0^1 \sin x\,dx = \bigl[-\cos x\bigr]_0^1 = 1-\cos 1 . $$ 数値では $1-\cos 1 = 0.4597\ldots$ です。 $a,b,c>0$ とする。楕円体 $$ \left\{(x,y,z)\ \middle|\ \frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}\le 1\right\} $$ の体積が $\dfrac43\pi abc$ であることを、変数変換 $\Phi(r,\theta)=(ar\cos\theta,\ br\sin\theta)$ を使って示してください。楕円体は、楕円板 $D=\left\{(x,y)\ \middle|\ \frac{x^2}{a^2}+\frac{y^2}{b^2}\le1\right\}$ の上で、グラフ $z=h(x,y)$ と $z=-h(x,y)$ に挟まれた立体です。ここで $$ h(x,y)=c\sqrt{1-\frac{x^2}{a^2}-\frac{y^2}{b^2}} . $$ したがって体積は $V=\displaystyle\iint_D 2h(x,y)\,dx\,dy$ です。 $\Phi(r,\theta)=(ar\cos\theta,\ br\sin\theta)$ のヤコビ行列と行列式は $$ J_\Phi=\begin{pmatrix} a\cos\theta & -ar\sin\theta\\ b\sin\theta & br\cos\theta\end{pmatrix}, \qquad \det J_\Phi = abr\cos^2\theta+abr\sin^2\theta = abr . $$ $E=[0,1]\times[0,2\pi]$ とすると $\Phi(E)=D$ で、単射性が崩れるのは $r=0$ と $\theta=0,2\pi$ の辺（面積 $0$）だけなのでより変数変換が使えます。また $$ h(\Phi(r,\theta)) = c\sqrt{1-r^2\cos^2\theta-r^2\sin^2\theta}=c\sqrt{1-r^2} $$ なので、とより $$ V = \int_0^{2\pi}\!\!\int_0^1 2c\sqrt{1-r^2}\cdot abr\,dr\,d\theta = 2abc\cdot 2\pi\int_0^1 r\sqrt{1-r^2}\,dr . $$ 最後の 1 変数積分は $t=1-r^2$（$dt=-2r\,dr$）と置換して $$ \int_0^1 r\sqrt{1-r^2}\,dr = \left[-\frac{(1-r^2)^{3/2}}{3}\right]_0^1 = \frac13 . $$ よって $V = 4\pi abc\cdot\dfrac13 = \dfrac43\pi abc$ です。$a=b=c$ とすれば球の体積 $\frac43\pi a^3$ に一致します。 $D$ を $(0,0),(1,0),(0,1)$ を頂点とする三角形、すなわち $D=\{(x,y)\mid x\ge0,\ y\ge0,\ x+y\le1\}$ とする。 $$ \iint_D \exp\left(\frac{y-x}{y+x}\right)dx\,dy $$ を求めてください（原点では被積分関数が定義されませんが、値を任意に $1$ と定めます）。 **下準備。** $(x,y)\in D$ かつ $(x,y)\ne(0,0)$ のとき $x+y>0$ で、$\lvert y-x\rvert\le x+y$ より指数の値は $[-1,1]$ に入ります。よって被積分関数は $[e^{-1},e]$ の範囲に収まり有界で、不連続点は原点と $\partial D$ だけ、すなわち面積 $0$ の集合に含まれます。の 1 より $D$ 上可積分です。 **変数変換。** $u=y-x$、$v=y+x$ と置きたいので、その逆写像 $$ \Phi(u,v)=\left(\frac{v-u}{2},\ \frac{v+u}{2}\right) $$ を考えます。$\Phi$ は $\mathbb{R}^2$ 全体の線形同型で $$ J_\Phi = \begin{pmatrix} -\tfrac12 & \tfrac12\\[2pt] \tfrac12 & \tfrac12\end{pmatrix}, \qquad \det J_\Phi = -\frac14-\frac14 = -\frac12, \qquad \lvert\det J_\Phi\rvert = \frac12 . $$ $\Phi(u,v)=(x,y)$ とすると、$x\ge0\iff u\le v$、$y\ge0\iff u\ge -v$、$x+y=v\le1$ なので $$ E=\{(u,v)\mid 0\le v\le 1,\ -v\le u\le v\} $$ が $\Phi$ によって $D$ に全単射に写ります（$E$ も $(0,0),(1,1),(-1,1)$ を頂点とする三角形です）。 **適用。** 被積分関数 $f$ は原点で連続でないので、$\eta\in(0,1)$ を取り $E_\eta=\{(u,v)\in E\mid \eta\le v\le1\}$ にを適用します。$E_\eta$ 上で $f(\Phi(u,v))=e^{u/v}$ は連続です。$E_\eta$ は縦線集合（$v$ を「横軸」と見て、各 $v$ に対し $u$ が $-v$ から $v$ まで）なのでが使え、 $$ \iint_{\Phi(E_\eta)} f\,dx\,dy = \frac12\int_\eta^1\!\left(\int_{-v}^{v} e^{u/v}\,du\right)dv = \frac12\int_\eta^1 v\Bigl[e^{u/v}\Bigr]_{u=-v}^{u=v}dv = \frac{e-e^{-1}}{2}\int_\eta^1 v\,dv . $$ （内側で $\int_{-v}^{v}e^{u/v}du = v\,(e^{1}-e^{-1})$ を使いました。$v>0$ なので $u\mapsto e^{u/v}$ の原始関数は $v\,e^{u/v}$ です。）一方 $D\setminus\Phi(E_\eta)$ は $\Phi(E\cap\{v\le\eta\})$ で、その面積は $\lvert\det J_\Phi\rvert\cdot\eta^2=\eta^2/2$（、$E\cap\{v\le\eta\}$ は面積 $\eta^2$ の三角形）。$\lvert f\rvert\le e$ なので、この部分の積分の寄与は $e\eta^2/2$ 以下で、$\eta\to0$ とすれば $0$ に収束します。したがって $$ \iint_D \exp\left(\frac{y-x}{y+x}\right)dx\,dy = \frac{e-e^{-1}}{2}\int_0^1 v\,dv = \frac{e-e^{-1}}{2}\cdot\frac12 = \frac{e-e^{-1}}{4}. $$ 数値では $0.5876\ldots$ です。$y-x$ と $y+x$ という「斜めの座標」に乗り換えることで、絡み合っていた変数が分離した点がこの問題の要点です。 ## 参考文献 - 杉浦光夫『解析入門 II』東京大学出版会、1985 — 重積分、フビニの定理、変数変換公式を扱う章。この記事で証明を省いた事実（リーマンの可積分条件、変数変換公式）の完全な証明があります。 - 高木貞治『解析概論』改訂第三版、岩波書店、1961 — 重積分と面積・体積を扱う章。古典的な記述で、累次積分の扱いが丁寧です。 - W. Rudin, *Principles of Mathematical Analysis*, 3rd ed., McGraw-Hill, 1976 — Chapter 10。primitive mapping への分解による変数変換定理の証明。 - M. Spivak, *Calculus on Manifolds*, W. A. Benjamin, 1965 — Chapter 3 (Integration)。フビニの定理と変数変換定理を、ルベーグ測度 $0$ の言葉で簡潔に扱っています。 - G. Fubini, "Sugli integrali multipli", *Rendiconti della Reale Accademia dei Lincei* 16 (1907) — 累次積分と重積分の関係を一般的に論じた原論文。 - T. Tao, *An Introduction to Measure Theory*, American Mathematical Society, 2011 — ルベーグ積分の枠組みでのフビニ・トネリの定理。非有界関数や無限区間を扱いたいときはこちらへ。 ## Appendix: 3 重積分と球座標ここまでの議論は次元を上げてもそのまま通用します。$\mathbb{R}^3$ の直方体を小直方体に分割してリーマン和を作れば 3 重積分が定義でき、と同じ論法で「3 重積分＝ 1 変数の積分 3 回」が示せます。変数変換公式も、の $\lvert\det A\rvert$ が「平行六面体の体積」を表すことに置き換えれば、同じ形 $$ \iiint_{\Phi(E)} f\,dx\,dy\,dz = \iiint_{E} f(\Phi(u,v,w))\,\lvert\det J_\Phi\rvert\,du\,dv\,dw $$ で成り立ちます。もっともよく使うのが球座標 $$ \Phi(r,\theta,\varphi) = (r\sin\theta\cos\varphi,\ r\sin\theta\sin\varphi,\ r\cos\theta) \qquad (r>0,\ 0<\theta<\pi,\ 0<\varphi<2\pi) $$ です（$\theta$ が天頂角、$\varphi$ が方位角）。ヤコビ行列は $$ J_\Phi = \begin{pmatrix} \sin\theta\cos\varphi & r\cos\theta\cos\varphi & -r\sin\theta\sin\varphi\\ \sin\theta\sin\varphi & r\cos\theta\sin\varphi & r\sin\theta\cos\varphi\\ \cos\theta & -r\sin\theta & 0 \end{pmatrix} $$ で、第 3 行に沿って余因子展開します。第 $(3,1)$ 余因子は $$ +\det\begin{pmatrix} r\cos\theta\cos\varphi & -r\sin\theta\sin\varphi\\ r\cos\theta\sin\varphi & r\sin\theta\cos\varphi\end{pmatrix} = r^2\sin\theta\cos\theta(\cos^2\varphi+\sin^2\varphi) = r^2\sin\theta\cos\theta, $$ 第 $(3,2)$ 余因子は $$ -\det\begin{pmatrix}\sin\theta\cos\varphi & -r\sin\theta\sin\varphi\\ \sin\theta\sin\varphi & r\sin\theta\cos\varphi\end{pmatrix} = -r\sin^2\theta(\cos^2\varphi+\sin^2\varphi) = -r\sin^2\theta $$ なので $$ \det J_\Phi = \cos\theta\cdot r^2\sin\theta\cos\theta + (-r\sin\theta)\cdot(-r\sin^2\theta) = r^2\sin\theta\,(\cos^2\theta+\sin^2\theta) = r^2\sin\theta . $$ $0<\theta<\pi$ では $\sin\theta>0$ なので絶対値も $r^2\sin\theta$ です。これを使えば半径 $a$ の球の体積は $$ \int_0^{2\pi}\!\!\int_0^{\pi}\!\!\int_0^{a} r^2\sin\theta\,dr\,d\theta\,d\varphi = 2\pi\cdot\bigl[-\cos\theta\bigr]_0^{\pi}\cdot\frac{a^3}{3} = 2\pi\cdot2\cdot\frac{a^3}{3} = \frac43\pi a^3 $$ と一息で求まります。の $r$ が「弧長が半径に比例する」ことの反映だったのと同じく、$r^2\sin\theta$ は「球面の面積要素が $r^2\sin\theta\,d\theta\,d\varphi$ である」ことの反映です。