第 1 步：画出区域 $D$ （配方法）—— 【很多考生的死穴】#

看到 $x^2 + y^2 \le x$，千万不要发懵！在解析几何中，只要看到 $x^2$ 和 $x$ 同时出现，第一反应必须是配方（平方完成）：

1. 移项：$x^2 - x + y^2 \le 0$
2. 配凑完全平方：$(x^2 - x + \frac{1}{4}) - \frac{1}{4} + y^2 \le 0$
3. 得到标准圆方程：$(x - \frac{1}{2})^2 + y^2 \le \left(\frac{1}{2}\right)^2$

几何直觉：这是一个圆心在 $(1/2, 0)$，半径为 $1/2$ 的实心圆盘。（如书上 图 5.11 所示，圆正好经过原点，且完全在 $y$ 轴的右侧，所以在这个区域内 $x \ge 0$ 恒成立）。

第 2 步：决定“切片”方向，写出上下限#

根据圆的方程 $y^2 \le x - x^2$，我们可以直接解出 $y$ 的上下边界： $-\sqrt{x - x^2} \le y \le \sqrt{x - x^2}$ 同时从图上可以看出，圆在 $x$ 轴上的投影范围是从 $0$ 到 $1$。

所以，我们选择**“先对 $y$ 积分，再对 $x$ 积分”**（竖着切）： $\text{外层（常数）}: 0 \le x \le 1$ $\text{内层（函数）}: -\sqrt{x - x^2} \le y \le \sqrt{x - x^2}$ 原式转化为累次积分： $I = \int_0^1 \left\{ \int_{-\sqrt{x-x^2}}^{\sqrt{x-x^2}} \sqrt{x} dy \right\} dx$

第 3 步：计算内层积分（对 $y$ 积分）#

我们要计算 $\int \sqrt{x} dy$。因为是对 $y$ 积分，$\sqrt{x}$ 被完全当作常数！ $\int_{-\sqrt{x-x^2}}^{\sqrt{x-x^2}} \sqrt{x} dy = \sqrt{x} \Big[ y \Big]_{-\sqrt{x-x^2}}^{\sqrt{x-x^2}} = \sqrt{x} \left( \sqrt{x-x^2} - (-\sqrt{x-x^2}) \right) = \sqrt{x} \cdot 2\sqrt{x-x^2}$

【神仙化简技巧】： $2\sqrt{x} \cdot \sqrt{x-x^2} = 2\sqrt{x \cdot x(1-x)} = 2\sqrt{x^2(1-x)}$ 因为在积分区域内 $x \ge 0$，所以 $\sqrt{x^2} = x$。 化简得到极度舒适的形式：$2x\sqrt{1-x}$。

第 4 步：计算外层积分（巧妙的一元代换）#

现在我们只需计算一元积分： $I = 2 \int_0^1 x\sqrt{1-x} dx$

怎么积？不要硬拆根号！使用一元微积分中最经典的换元法（置換積分）：

• 令 $t = 1 - x$ （这能把根号里面讨厌的多项式变成单项式）。
• 那么 $x = 1 - t$。
• 求导数：$dx = -dt$。
• 换上下限：当 $x=0$ 时，$t=1$；当 $x=1$ 时，$t=0$。

代入积分式： $I = 2 \int_1^0 (1-t)\sqrt{t} (-dt)$ 负号刚好可以用来把积分的上下限颠倒回来（变成从 0 到 1）： $I = 2 \int_0^1 (1-t)t^{1/2} dt = 2 \int_0^1 (t^{1/2} - t^{3/2}) dt$ 这就是极其简单的幂函数积分了： $= 2 \left[ \frac{2}{3}t^{3/2} - \frac{2}{5}t^{5/2} \right]_0^1$ 把 $1$ 代入： $= 2 \left( \frac{2}{3} - \frac{2}{5} \right) = 2 \left( \frac{10 - 6}{15} \right) = 2 \times \frac{4}{15} = \mathbf{\frac{8}{15}}$

大功告成！

考场解题步骤：#

第 1 步：分析区域与定序 题目已经非常仁慈地把上下限给我们写好了！这是一个标准的“横向切割（先 $x$ 后 $y$）”区域。 外层积分限是常数：$y$ 从 $0$ 到 $1$。 内层积分限是函数：$x$ 从 $y$ 到 $10y$。 原式转化为： $I = \int_0^1 \left\{ \int_y^{10y} \sqrt{xy - y^2} dx \right\} dy$

第 2 步：计算内层积分（对 $x$ 积分，视 $y$ 为常数） 我们要算 $\int \sqrt{xy - y^2} dx$。为了不看晕，我们把里面提取一个 $y$ 出来： $\int \sqrt{y(x - y)} dx = \sqrt{y} \int \sqrt{x - y} dx$ 因为 $y$ 是常数，$\int \sqrt{x-y} dx$ 直接套用 $\int u^{1/2} du = \frac{2}{3}u^{3/2}$： $= \sqrt{y} \cdot \frac{2}{3}(x - y)^{3/2}$ (注：书后答案写的是 $\frac{2}{3}(xy-y^2)^{3/2} \cdot \frac{1}{y}$，这与我们提公因式的做法是完全等价的，提公因式不容易算错！)

代入上下限 $x = 10y$ 和 $x = y$：

• 上限 $x = 10y$ 代入：$\sqrt{y} \cdot \frac{2}{3}(10y - y)^{3/2} = \sqrt{y} \cdot \frac{2}{3}(9y)^{3/2} = y^{1/2} \cdot \frac{2}{3} \cdot 27 \cdot y^{3/2} = 18y^2$
• 下限 $x = y$ 代入：$(y - y)^{3/2} = 0$ 所以内层积分的结果非常清爽：$18y^2$。

第 3 步：计算外层积分 $I = \int_0^1 18y^2 dy = \left[ 6y^3 \right]_0^1 = \mathbf{6}$ （答案完美命中官方解答）

考场解题步骤：#

第 1 步：极坐标替换法则（⚠️ 千万别漏掉雅可比式 $r$！） 只要看到带有 $x^2+y^2$ 或者极坐标方程 $r = f(\theta)$，马上条件反射：

• $x = r\cos\theta, \ y = r\sin\theta$
• 微元替换：$dxdy = \mathbf{r} \, drd\theta$ （这个多出来的 $r$ 每年都有无数人忘记！）

第 2 步：转化被积函数与上下限

• 被积函数：$y$ 变成 $r\sin\theta$。加上微元的 $r$，总的被积式变成了 $r^2 \sin\theta$。
• 积分区域：因为是上半部分（$y \ge 0$），所以角度 $\theta$ 从 $0$ 转到 $\pi$。 半径 $r$ 从原点 $0$ 伸展到心形线的边界 $a(1+\cos\theta)$。 $I = \int_0^\pi \left\{ \int_0^{a(1+\cos\theta)} r^2 \sin\theta dr \right\} d\theta$

第 3 步：计算内层积分（对 $r$ 积分） $\int_0^{a(1+\cos\theta)} r^2 \sin\theta dr = \sin\theta \left[ \frac{r^3}{3} \right]_0^{a(1+\cos\theta)} = \frac{a^3}{3} (1+\cos\theta)^3 \sin\theta$

第 4 步：计算外层积分（一元换元法） $I = \frac{a^3}{3} \int_0^\pi (1+\cos\theta)^3 \sin\theta d\theta$ 【换元技巧】：令 $t = 1+\cos\theta$。 则 $dt = -\sin\theta d\theta$。 上下限：当 $\theta = 0$ 时 $t = 2$；当 $\theta = \pi$ 时 $t = 0$。 $I = \frac{a^3}{3} \int_2^0 t^3 (-dt) = \frac{a^3}{3} \int_0^2 t^3 dt$ $= \frac{a^3}{3} \left[ \frac{t^4}{4} \right]_0^2 = \frac{a^3}{3} \cdot \frac{16}{4} = \mathbf{\frac{4}{3}a^3}$ （答案完美命中官方解答）

考场解题步骤：#

第 1 步：画图并寻找交点（关键！） 在草稿纸上画出：

1. 原点出发的斜线 $y = x/2$ （即 $x = 2y$）
2. 斜向下的直线 $y = -x + 3$ （即 $x = 3-y$）
3. 底部的 $x$ 轴（$y=0$） 这三条线围成了一个三角形（对应书后第 256 页图 3.1(1) 下方的图 2.1(3)）。

• 找顶点：联立 $x=2y$ 和 $x=3-y$，解得 $2y = 3-y \implies 3y = 3 \implies y=1, x=2$。 顶点坐标为 $(2, 1)$。

第 2 步：决定切割方向

• 如果竖着切（先 $y$ 后 $x$）：天花板在 $x=2$ 的左右两侧不一样（左边是 $y=x/2$，右边是 $y=3-x$），必须把积分拆成两半算，非常麻烦！
• 【满分解法】横着切（先 $x$ 后 $y$）：从下往上看，$y$ 的范围是 $0$ 到最高点 $1$。在这个范围内，左墙壁永远是 $x=2y$，右墙壁永远是 $x=3-y$。一气呵成！ $I = \int_0^1 \left\{ \int_{2y}^{3-y} x dx \right\} dy$

第 3 步：计算内层积分（对 $x$ 积分） $\int_{2y}^{3-y} x dx = \left[ \frac{x^2}{2} \right]_{2y}^{3-y} = \frac{1}{2} \left( (3-y)^2 - (2y)^2 \right)$ 展开化简： $\frac{1}{2} (9 - 6y + y^2 - 4y^2) = \frac{1}{2} (9 - 6y - 3y^2)$

第 4 步：计算外层积分 $I = \frac{1}{2} \int_0^1 (9 - 6y - 3y^2) dy$ $= \frac{1}{2} \left[ 9y - 3y^2 - y^3 \right]_0^1$ 代入上限 1： $= \frac{1}{2} (9 - 3 - 1) = \frac{1}{2} \times 5 = \mathbf{\frac{5}{2}}$ （答案完美命中官方解答）

1. 雅可比行列式的计算公式#

$J = \frac{\partial(x, y)}{\partial(u, v)} = \det \begin{pmatrix} \frac{\partial x}{\partial u} & \frac{\partial x}{\partial v} \\ \frac{\partial y}{\partial u} & \frac{\partial y}{\partial v} \end{pmatrix} = \frac{\partial x}{\partial u} \frac{\partial y}{\partial v} - \frac{\partial x}{\partial v} \frac{\partial y}{\partial u}$ (口诀：横着写偏导，主对角线相乘 减去 副对角线相乘)。

2. 二重积分换元终极公式（式 1）#

$\iint_D f(x, y) dxdy = \iint_{D'} f(\varphi(u, v), \psi(u, v)) |J| dudv$

🚨【考场最高危、最致命的陷阱】🚨 仔细看书上公式的最后一个项：是 $|J|$ 而不是 $J$！ 雅可比行列式外面套着一个绝对值符号！ 很多考生在考场上算出 $J = -2$，然后直接带着负号扔进积分里，导致体积算出来是负数，整道大题直接 0 分。 原理：$J$ 的正负只代表坐标系翻转了方向，但 $dxdy$ 代表的是一小块微元的面积，面积永远是正的！所以代入公式时，必须加上绝对值保证面积缩放倍数为正！

3. 【顿悟时刻】：极坐标的那个 $r$ 到底是怎么来的？#

在讲第 136 页时我再三强调极坐标代换 $dxdy = r \, drd\theta$，那个 $r$ 就是雅可比行列式算出来的！ 我们现场算一下极坐标的 Jacobian： 已知 $x = r\cos\theta, \ y = r\sin\theta$ （这里 $u=r, v=\theta$）。 $J = \begin{vmatrix} x_r & x_\theta \\ y_r & y_\theta \end{vmatrix} = \begin{vmatrix} \cos\theta & -r\sin\theta \\ \sin\theta & r\cos\theta \end{vmatrix}$ 算行列式：$(\cos\theta)(r\cos\theta) - (-r\sin\theta)(\sin\theta) = r\cos^2\theta + r\sin^2\theta = \mathbf{r}$ 因为半径 $r \ge 0$，所以绝对值 $|J| = r$。 这就是极坐标代换多出来一个 $r$ 的底层理论来源！

1. 为什么需要“扩展”？#

在第 132 页讲黎曼和定义时，我们有一个大前提：积分区域 $D$ 必须是有界闭区域，且函数 $f(x,y)$ 必须连续（不能有无穷大的点）。 但在考研中，教授经常出这种题：

• 求 $\iint_D \frac{1}{\sqrt{x^2+y^2}} dxdy$，区域是原点附近的圆。 (致命问题：当 $(x,y) \to (0,0)$ 时，分母为 0，函数值爆炸到无穷大！)

2. “近似增加列”的思想（图 5.15）#

为了解决这种含有**不连续点/特异点（集合 $E$）**的积分，书本引入了非常严密的极限思想：

1. 挖坑：既然集合 $E$ 上的点（比如原点）会让函数爆炸，那我们就先拿个保护罩把它“罩起来”并挖掉。
2. 构建近似列 $\{D_n\}$：挖掉坏点后剩下的安全区域记为 $D_n$（图 5.15 中蓝色的部分）。在这个安全区域内，函数是连续的，我们可以放心大胆地算普通二重积分。
3. 逼近：让 $n \to \infty$（即让挖掉的保护罩越来越小，无限逼近那个坏点 $E$）。 如果这个逼近过程算出来的极限存在，那么这个极限值就是广义二重积分的值！

（书上在底部列出了构建 $D_n$ 必须满足的三个数学条件：包容性、避开奇点、面积逼近。这是为了保证数学严密性，考场上主要考的是实战计算。）

解决思想：近似增加列（Expanding Exhaustion）#

和对付奇点的“挖洞逼近”相反，对付无限大区域，我们用的是**“膨胀逼近”**（如 图 5.17 所示）：

1. 先在原点附近画一个半径为 $n$ 的四分之一圆盘，记为 $D_n$（这就是一个有界闭区域了）。
2. 在 $D_n$ 里面老老实实算二重积分。
3. 算完之后，让半径 $n \to \infty$。圆盘就会无限膨胀，最终覆盖整个第一象限！取这个极限，就是无限区域上的积分值。

1. 定理 5.6 到底在说什么？#

假设被积函数 $f(x, y) \ge 0$（比如求体积、概率密度等）： 如果你用某种特定形状（比如正方形）去逼近无限区域，算出极限存在且等于 $I$。 那么，无论你换用其他任何形状（比如圆形、三角形、星形）去膨胀逼近，算出来的极限一定也是 $I$！

2. 💡【考场提分神技：随心所欲选形状】#

在考卷上让你积整个第一象限 $\iint_{x \ge 0, y \ge 0} e^{-(x^2+y^2)} dxdy$。

• 老实的同学：画一个正方形区域 $0 \le x \le n, 0 \le y \le n$。算完后发现极其难算。
• 懂定理 5.6 的同学：既然被积函数 $e^{-(x^2+y^2)} > 0$ 恒成立，那根据定理 5.6，我凭什么非要用正方形逼近？我偏要用四分之一圆盘（$x^2+y^2 \le n^2$）去逼近！
• 降维打击：一旦用了圆盘，瞬间就可以引入极坐标代换 $x^2+y^2 = r^2$ 且 $dxdy = rdrd\theta$。原本极难的指数积分，瞬间化简为极其简单的普通一元积分！

考试金律：只要被积函数是非负的（$f \ge 0$），在无限区域上求积分时，哪种形状的边界容易算，你就大胆地用哪种形状去逼近！ 绝大多数情况下，使用圆形区域膨胀 + 极坐标换元是秒杀无限区域重积分的唯一正解。

第 1 步：祭出“椭圆代换神技”（Route）#

书本 route 栏给出了处理椭圆的千古名招： 令 $x = au, \quad y = bv$

• 几何意义（图 5.18）：把这个代换放进椭圆方程里，$\frac{(au)^2}{a^2} + \frac{(bv)^2}{b^2} = u^2 + v^2 \le 1$。 我们成功地把一个扁扁的椭圆 $D$，拉伸成了一个完美的单位实心圆 $D'$！

第 2 步：计算雅可比行列式（Jacobian）#

微元替换绝不是简单的 $dxdy = dudv$，必须乘上面积缩放比例 $|J|$。 $J = \frac{\partial(x, y)}{\partial(u, v)} = \begin{vmatrix} \frac{\partial x}{\partial u} & \frac{\partial x}{\partial v} \\ \frac{\partial y}{\partial u} & \frac{\partial y}{\partial v} \end{vmatrix} = \begin{vmatrix} a & 0 \\ 0 & b \end{vmatrix} = ab - 0 = ab$ 因为 $a>0, b>0$，所以绝对值 $|J| = ab \neq 0$。 于是，微元替换公式为：$dxdy = ab \, dudv$。

第 3 步：转化原二重积分#

把 $x, y$ 和 $dxdy$ 全部换成 $u, v$ 的世界： $I = \iint_{D'} \left( (au)^2 + (bv)^2 \right) \cdot (ab) \, dudv = ab \iint_{D'} (a^2u^2 + b^2v^2) dudv$ 此时积分区域 $D'$ 已经是单位圆 $u^2+v^2 \le 1$ 了。

第 4 步：化简积分并利用对称性#

由于积分区域是单位圆（上下左右完全对称），且被积函数 $a^2u^2 + b^2v^2$ 全是偶次幂。 【考场省时技巧】：我们不需要积整个圆，只需要积第一象限（右上角的四分之一圆），然后乘以 4 即可！ 第一象限的界限是：$0 \le u \le 1$，且 $0 \le v \le \sqrt{1-u^2}$。 $I = 4ab \int_0^1 \left\{ \int_0^{\sqrt{1-u^2}} (a^2u^2 + b^2v^2) dv \right\} du$

先对 $v$ 积分（把 $u$ 当常数）： $\int_0^{\sqrt{1-u^2}} (a^2u^2 + b^2v^2) dv = \left[ a^2u^2 v + \frac{b^2}{3}v^3 \right]_0^{\sqrt{1-u^2}} = a^2u^2\sqrt{1-u^2} + \frac{b^2}{3}(1-u^2)^{3/2}$

第 5 步：一元三角换元（硬骨头计算）#

把上面的结果代入外层积分： $I = 4ab \int_0^1 \left\{ a^2u^2\sqrt{1-u^2} + \frac{b^2}{3}(1-u^2)^{3/2} \right\} du$ 看到 $\sqrt{1-u^2}$，本能反应是三角换元：令 $u = \sin\theta$。 则 $du = \cos\theta d\theta$；上下限从 $u \in [0, 1]$ 变成 $\theta \in [0, \pi/2]$。 代入后，$\sqrt{1-u^2}$ 变成了 $\cos\theta$。积分式化为了纯三角函数积分： $I = 4ab \int_0^{\pi/2} \left\{ a^2 \sin^2\theta \cos\theta + \frac{b^2}{3} \cos^3\theta \right\} \cdot \cos\theta d\theta$ 展开后得到： $I = 4ab \int_0^{\pi/2} \left( a^2 \sin^2\theta \cos^2\theta + \frac{b^2}{3} \cos^4\theta \right) d\theta$ 利用 $\sin^2\theta \cos^2\theta = \cos^2\theta(1-\cos^2\theta) = \cos^2\theta - \cos^4\theta$，可以整理成书上的形式： $= 4ab \int_0^{\pi/2} \left\{ a^2\cos^2\theta + \left(\frac{b^2}{3}-a^2\right)\cos^4\theta \right\} d\theta$

【考研绝招：华里士公式 Wallis Formula】 到了这一步千万别用倍角公式硬拆！书上右侧旁注提示了 p.82 第3章 例题25。 记住这个能救命的结论：$\int_0^{\pi/2} \cos^2\theta d\theta = \frac{\pi}{4}$，且 $\int_0^{\pi/2} \cos^4\theta d\theta = \frac{3\pi}{16}$。 直接把这两个数值代入上面的式子，展开化简后，瞬间得到最终答案： $I = \mathbf{\frac{\pi ab}{4}(a^2+b^2)}$

第 1 步：转化积分区域 $D \to D'$（图 5.19 核心奥秘）#

• 原来的区域 $D$：是在 $x-y$ 坐标系下的一个实心圆（$x^2 + y^2 \le a^2$）。
• 代换法则：令 $x = r\cos\theta, \ y = r\sin\theta$。这代表平面上任何一个点，都可以用距离原点的半径 $r$ 和旋转角度 $\theta$ 来表示。
• 新区域 $D'$ 是什么？ 要画满一个半径为 $a$ 的实心圆： 角度 $\theta$ 需要转整整一圈：$0 \le \theta \le 2\pi$ 半径 $r$ 需要从圆心一直伸到边缘：$0 \le r \le a$
• 【几何奇迹】：大家看图 5.19 的左边！在新的 $r-\theta$ 坐标系里，这两个条件画出来，居然是一个完美的矩形 $D'$！ 极坐标换元的本质，就是把弯曲的、难切的圆形区域，扯平拉直成了一个上下左右全是常数的矩形！这就是为什么换元后极其好算的原因。

第 2 步：计算雅可比行列式 $J$ （理论的必经之路）#

我们在前一页强调过，微元变换不是直接换，必须乘以面积缩放系数（Jacobian）。书本在这里给出了严格的推导： $J = \frac{\partial(x, y)}{\partial(r, \theta)} = \begin{vmatrix} \frac{\partial x}{\partial r} & \frac{\partial x}{\partial \theta} \\ \frac{\partial y}{\partial r} & \frac{\partial y}{\partial \theta} \end{vmatrix} = \begin{vmatrix} \cos\theta & -r\sin\theta \\ \sin\theta & r\cos\theta \end{vmatrix}$ 计算这个行列式： $J = (\cos\theta)(r\cos\theta) - (-r\sin\theta)(\sin\theta) = r\cos^2\theta + r\sin^2\theta = \mathbf{r}$ 因为半径 $r \ge 0$，所以绝对值 $|J| = r$。 【考场铁律】：考试时绝对不需要重新推导这一步！直接把 $dxdy = r \, drd\theta$ 当成固定公式默写上去！千万千万别漏掉那个 $r$！

（注：书上在 (i) 里提到了一句极其严谨的废话：虽然当 $r = 0$ 时 $\theta$ 可以是任何值，导致 $1$ 对 $1$ 映射在原点失效，但这只是一条线/一个点，它的面积是 $0$，所以对积分结果完全没有影响。）

第 3 步：在新区域上完成“矩形积分”#

把所有东西扔进新的积分号里：

• 被积函数：$x^2 + y^2$ 直接变成了 $r^2$。
• 面积微元：$dxdy$ 变成了 $r \, drd\theta$。
• 积分上下限：变成了全常数。

$\iint_D (x^2 + y^2) dxdy = \iint_{D'} (r^2) \cdot (r) \, drd\theta = \int_0^{2\pi} \left\{ \int_0^a r^3 dr \right\} d\theta$ 【考场秒杀技巧：变量分离法】 只要积分上下限全是常数，且被积函数是 $r$ 和 $\theta$ 乘积的形式，我们完全可以直接把两个积分拆开，分别算完再相乘！（书上就是这么做的）： $= \left( \int_0^{2\pi} 1 d\theta \right) \times \left( \int_0^a r^3 dr \right) = \Big[ \theta \Big]_0^{2\pi} \times \left[ \frac{r^4}{4} \right]_0^a = 2\pi \times \frac{a^4}{4} = \mathbf{\frac{\pi a^4}{2}}$ 计算瞬间完成，清爽无比！

考场满分 4 步曲：#

• 第 1 步：构造“近似增加列（膨胀的圆盘）” 我们取原点为圆心，画一个半径为 $n$ 的四分之一圆盘，记为 $D_n$（如 图 5.21 所示）。 $D_n: x \ge 0, \ y \ge 0, \ x^2+y^2 \le n^2$ 随着 $n \to \infty$，$D_n$ 将覆盖整个第一象限。

• 第 2 步：在 $D_n$ 上引入极坐标代换 既然是圆盘，立刻写出极坐标三连： $x = r\cos\theta, \ y = r\sin\theta$ $dxdy = \mathbf{r} \, drd\theta$ （雅可比行列式千万别忘！） 积分限：$0 \le \theta \le \pi/2$ (第一象限)， $0 \le r \le n$ (半径从 0 到 $n$)。

• 第 3 步：计算 $D_n$ 上的积分（极其巧妙的凑微分） 把式子代入： $\iint_{D_n} e^{-x^2-y^2} dxdy = \int_0^{\pi/2} \int_0^n e^{-r^2} \cdot r \, drd\theta$ 变量分离，各自相乘： $= \left( \int_0^{\pi/2} 1 d\theta \right) \times \left( \int_0^n r e^{-r^2} dr \right)$ 前半部分是 $\frac{\pi}{2}$。后半部分用凑微分（令 $u = -r^2, du = -2rdr$）： $\int r e^{-r^2} dr = -\frac{1}{2} \int e^{-r^2} d(-r^2) = -\frac{1}{2} e^{-r^2}$。 代入上下限 $n$ 和 $0$： $\left[ -\frac{1}{2} e^{-r^2} \right]_0^n = -\frac{1}{2} e^{-n^2} - \left(-\frac{1}{2} e^0\right) = \frac{1}{2}(1 - e^{-n^2})$ 两者相乘，我们在 $D_n$ 上的积分为：$\frac{\pi}{4}(1 - e^{-n^2})$。

• 第 4 步：让圆盘无限膨胀，取极限 $I = \lim_{n\to\infty} \frac{\pi}{4}(1 - e^{-n^2})$ 因为当 $n \to \infty$ 时，$-n^2 \to -\infty$，所以 $e^{-n^2} \to 0$。 最终答案完美得出：$I = \frac{\pi}{4}$。

⚠️【考场危机分析（双重炸弹）】#

这道积分题极其凶险，它面临着两个致命问题：

1. 边界爆炸（奇点）：如果 $p < 1$ 或 $q < 1$，那么 $x^{p-1}$ 和 $y^{q-1}$ 在 $x=0$ 和 $y=0$（坐标轴）处会变成负指数，导致分母为 0，函数值爆炸无穷大。
2. 区域无限大（无界）：积分区域 $D$ 是整个第一象限，无穷无尽。

🛡️【破局法宝：万能的定理 5.6】#

因为被积函数 $e^{-x-y} x^{p-1}y^{q-1} > 0$ 恒为正，我们再次祭出上一页的“免死金牌”定理 5.6：我们可以用任何形状的安全区域去逼近它，最终算出来的极限值一定相等！ 这就为我们“用两种方法算同一个积分，然后画等号”提供了理论合法性。

1. 划定安全区域 $D_n$（图 5.23 解析）#

为了避开“坐标轴”这个奇点，又为了“无限膨胀”覆盖第一象限，书上构建了一个非常精妙的正方形区域 $D_n$： $D_n: \frac{1}{n} \le x \le n, \quad \frac{1}{n} \le y \le n$ 您看 图 5.23 那个蓝色的方块，它故意和 $x$ 轴、$y$ 轴保留了一点点缝隙（$1/n$），同时向外延伸到 $n$。 当 $n \to \infty$ 时，缝隙会闭合，外边界会飞向无穷，它完美覆盖了整个第一象限！

2. 分离变量，拆解积分#

把被积函数拆开：$e^{-x-y} x^{p-1}y^{q-1} = (e^{-x} x^{p-1}) \cdot (e^{-y} y^{q-1})$。 因为在正方形区域上，上下限都是常数，所以二重积分可以直接“劈成两半”，变成两个独立的一元积分相乘： $\iint_{D_n} = \left( \int_{1/n}^n e^{-x} x^{p-1} dx \right) \times \left( \int_{1/n}^n e^{-y} y^{q-1} dy \right)$

3. 取极限，见证 Gamma 函数的诞生#

现在，我们让小缝隙闭合，让正方形无限放大（令 $n \to \infty$）。

• 左边的括号变成了：$\int_0^\infty e^{-x} x^{p-1} dx$。 这不就是第 89 页定义过的 伽马函数 $\Gamma(p)$ 吗！
• 右边的括号变成了：$\int_0^\infty e^{-y} y^{q-1} dy$。 这就是 伽马函数 $\Gamma(q)$！

因此，我们得到了 式 (1)： 通过正方形区域逼近，这个二重积分的最终值为： $\iint_D e^{-x-y} x^{p-1}y^{q-1} dxdy = \Gamma(p) \Gamma(q)$