曲线/曲面积分

形式上有

实际计算时, 为曲线 `L` (曲面 `Sigma`) 选取合适的参数, 就有 `"d"bm r = bm v dt`, `quad "d"s = |bm v| dt`;
`"d"bm S = bm n "d"u"d"v`, `quad "d"sigma = |bm n|"d"u"d"v`. 其中 切向量 `bm v = ("d"bm r)/dt` `= (dx/dt, dy/dt, dz/dt)`, `quad` 法向量 `bm n = bm r_u xx bm r_v` `= |bm i, bm j, bm k; x_u, y_u, z_u; x_v, y_v, z_v|` 依赖于参数选取.

1. 把曲面的单位法向量记为 `(cos alpha, cos beta, cos gamma)`, 则 `"d"bm S = (cos alpha, cos beta, cos gamma) "d"sigma`.
2. 对于平面曲线, 其法向微元定义为 `"d"bm n = (dy, -dx)`, 从而 `("d"bm r, "d"bm n)` 成左手系, 即 `|dx, dy; dy, -dx| = -"d"s^2 lt 0`.

求 `int_L x y "d"s`, 其中 `L: {x+y+z = 0; x^2+y^2+z^2=1:}`

利用积分区域的对称性, `int_L x y "d"s` `= 1/3 int_L (x y+y z+z x) "d"s` `= 1/6 int_L [(x+y+z)^2 - (x^2+y^2+z^2)] "d"s` `= -1/6 int_L "d"s` `= -pi/3`.

微分算子

偏微分算子 `del_x := del/(del x)`, `del_i := del/(del x_i)`.

Nabla 算子

标量场 (数量场) `f` 具有梯度 `"grad" f = nabla f`. 梯度是矢量, 其与坐标选取无关. 梯度是导数的推广, 它简单地将求导运算推广到各分量, 其加减乘除及函数复合的运算和导数完全一致.

方向导数定义为函数值沿某个方向的变化率: `("d"f)/("d"bm n) = lim_(t to 0) (f(bm x_0 + t bm n) - f(bm x_0))//t` 可以证明 `("d"f)/("d"bm n) = bm n * nabla f`, 因此梯度方向是该点处方向导数取得最大的方向, 大小等于这个最大的方向导数.

矢量场 (向量场) `bm F` 具有散度 `"div" bm F = nabla * bm F` `= (del F_1)/(del x_1) + cdots + (del F_n)/(del x_n)`. 和旋度 (三维空间): `"rot" bm F = nabla xx bm F` `= |bm i, bm j, bm k; del/(del x), del/(del y), del/(del z); F_1, F_2, F_3|`. 散度是标量而旋度是矢量.

散度定义为封闭曲面包围区域趋于无穷小时, 单位体积的通量之极限, 规定向外为正方向: `"div" bm F = lim_("diam"V to 0) (iint_(del V) bm F * "d"bm S)/V`. 散度的符号指出矢量场在一点处有正源/负源, 散度的大小指出源的强度.

方向旋量定义为给定法向量的平面上单位面积的环量之极限, 规定右手螺旋为正方向: `"rot"_(bm v) bm F = lim_("diam"S to 0) (int_(del S) bm F * "d"bm r)/S`. `bm F` 与切向 `"d"bm r` 方向越近似, 环量就越是取正值. 因此环量反映了矢量场沿曲线是否有旋涡, 以及旋涡的强度. 方向旋量反映了矢量场在一点处沿给定方向的旋转情况, 显然只要 `bm v` 取反方向, 方向旋量的符号就取反. 由 Stokes 公式可以证明 `"rot"_(bm v) bm F = bm v * "rot" bm F`. 旋度的方向是该点处方向旋量取得最大的方向, 大小等于这个最大的方向旋量.

设 `f`, `bm F` 二阶连续可微, 因此二阶混合偏导相等. 梯度场的散度是 Laplace 算子 (见下节); 梯度场的旋度为零向量, 旋度场的散度为 0. `nabla * (nabla f) = (nabla * nabla) f := nabla^2 f`,
`nabla xx (nabla f) = bm 0`, `quad nabla * (nabla xx bm F) = 0`.

`bm A * nabla` 算子可以用于标量或矢量. 输入标量, 输出标量; 输入矢量, 输出矢量: `(bm A * nabla) f` `= A_1 (del f)/(del x_1) + cdots + A_n (del f)/(del x_n)`,
`(bm A * nabla) bm F` `= ((bm A * nabla) F_1, cdots, (bm A * nabla) F_n)`.

注意 `bm A * nabla` 是一个算子, 而 `nabla * bm A` 是一个标量. 比如在三维情形: `bm A * nabla = A_1 del/(del x) + A_2 del/(del y) + A_3 del/(del z)`,
`nabla * bm A = pp A_1 x + pp A_2 y + pp A_3 z`.

`D` 算子作用于矢量, 得到它的导矩阵: `D bm F` `= [nabla F_1; vdots; nabla F_n]` `= [ (del F_1)/(del x_1), cdots, (del F_1)/(del x_n); vdots, , vdots; (del F_n)/(del x_1), cdots, (del F_n)/(del x_n); ]` 作用于标量时, 就是 `D f = nabla f`.

nabla 算子的综合公式 以下公式都涉及旋度, 因此只限于三维空间:
1. `nabla(bm F * bm G)` `= bm F D bm G + bm G D bm F` `= bm F xx (nabla xx bm G) + bm G xx (nabla xx bm F) + (bm F * nabla) bm G + (bm G * nabla) bm F`;
2. `nabla * (bm F xx bm G) = bm G * nabla xx bm F + bm F * nabla xx bm G`;
3. `nabla xx (bm F xx bm G) = (nabla * bm G) bm F - (nabla * bm F) bm G + (bm G * nabla) bm F - (bm F * nabla) bm G`;
4. `nabla xx (nabla xx bm F) = nabla(nabla * bm F) - (nabla * nabla)bm F`

直接计算验证. 注意到各个分量的地位对称, 因此只验证 `x` 分量相等即可. 以 1. 的证明为例, 左边的 `x` 分量等于 `F_1 (del G_1)/(del x) + F_2 (del G_2)/(del x) + F_3 (del G_3)/(del x)` `+ ** + ** + **`. 右边的 `x` 分量等于 `|F_2, F_3; (del G_1)/(del z) - (del G_3)/(del x), (del G_2)/(del x) - (del G_1)/(del y)|` `+ (F_1 del/(del x) + F_2 del/(del y) + F_3 del/(del z)) G_1` `+ | ** |` `+ **`. 其中 `**` 号省略的部分, 只需将字母 `F, G` 对换就可得到. 因此两式相等.

Laplace 算子

如 `nabla^2(f g) = nabla^2(f g)` `= nabla * (f nabla g + g nabla f)` `= f nabla^2 g + g nabla^2 f + 2 nabla f * nabla g`.

满足 `nabla^2 f -= 0` 的函数称为调和函数.

不同坐标系中的微分算子

部分计算参见多元函数微分学.

其中极坐标下 `bm e_r = (cos theta, sin theta)`, `quad bm e_theta = (-sin theta, cos theta)`. 球坐标下 `bm e_r = (cos varphi cos theta, cos varphi sin theta, sin varphi)`,
`bm e_theta = (-sin theta, cos theta, 0)`,
`bm e_varphi = bm e_r xx bm e_theta` `= (-sin varphi cos theta, -sin varphi sin theta, cos varphi)`. 有些文献中, 球坐标系的 `varphi` 不是指纬度, 而是指与北极方向的夹角. 此时应当把结论中的 `cos varphi`, `sin varphi` 互换, 并把 `del_varphi` 换成 `-del_varphi`.

1. 极坐标的 nabla 算子. 根据一阶微分的形式不变性, 可以通过求解线性方程组完成此类问题, 二阶微分则没有这种做法. 列出 `del_r = cos theta del_x + sin theta del_y`,
   `del_theta = -r sin theta del_x + r cos theta del_y`, Jacobi 行列式 `(del(x,y))/(del(r, theta))` `= |cos theta, sin theta; -r sin theta, r cos theta|` `= r`. 解方程组得 `del_x = (r cos theta del_r - sin theta del_theta) // r`,
   `del_y = (r sin theta del_r + cos theta del_theta) // r`.
2. 球坐标的 nabla 算子. 列出 `del_r = cos varphi cos theta del_x + cos varphi sin theta del_y + sin varphi del_z`,
   `del_theta = -r cos varphi sin theta del_x + r cos varphi cos theta del_y`,
   `del_varphi = -r sin varphi cos theta del_x - r sin varphi sin theta del_y + r cos varphi del_z`, Jacobi 行列式 `(del(x,y,z))/(del(r,theta,varphi))` `= |cos varphi cos theta, cos varphi sin theta, sin varphi; -r cos varphi sin theta, r cos varphi cos theta, 0; -r sin varphi cos theta, -r sin varphi sin theta, r cos varphi|` `= r^2 cos varphi`. 解方程组 [附代码] 得 `del_x = cos varphi cos theta del_r - 1/r (sin varphi cos theta del_varphi + (sin theta)/(cos varphi) del_theta)`,
   `del_y = cos varphi sin theta del_r - 1/r (sin varphi sin theta del_varphi - (cos theta)/(cos varphi) del_theta)`,
   `del_z = sin varphi del_r + 1/r cos varphi del_varphi`.

   ---

   In[1]: [simplify(i) for i in solve([
     cos(p)*cos(t)*x+cos(p)*sin(t)*y+sin(p)*z-dr,
     -r*cos(p)*sin(t)*x+r*cos(p)*cos(t)*y-dt,
     -r*sin(p)*cos(t)*x-r*sin(p)*sin(t)*y+r*cos(p)*z-dp
   ], [x, y, z]).values()]
   Out[1]:
   [(-dp*sin(p)*cos(t) + dr*r*cos(p)*cos(t) - dt*sin(t)/cos(p))/r,
   (-dp*sin(p)*sin(t) + dr*r*sin(t)*cos(p) + dt*cos(t)/cos(p))/r,
   dp*cos(p)/r + dr*sin(p)]
       

1. 极坐标的 Laplace 算子. 一个误区是直接计算 `nabla^2` `= (bm e_r del_r + 1/r bm e_theta del_theta)^2` `= del_r^2 + 1/r^2 del_theta^2`. 这种做法的问题在于忽视了 `bm e_r, bm e_theta` 是 `theta` 的函数, 正确做法是: `nabla^2` `= bm e_r del_r * (bm e_r del_r + 1/r bm e_theta del_theta)` `+ 1/r bm e_theta del_theta * (bm e_r del_r + 1/r bm e_theta del_theta)`. 先求 `del_theta bm e_r = bm e_theta`, `quad del_theta bm e_theta = - bm e_r`, 接下来利用 `bm e_r * bm e_theta = 0`, 分别计算 `bm e_r del_r * (bm e_r del_r)` `= bm e_r * del_r (bm e_r del_r)` `= (bm e_r * bm e_r) del_r^2` `= del_r^2`,
   `bm e_r del_r * (1/r bm e_theta del_theta)` `= bm e_r * bm e_theta del_r (1/r del_theta)` `= 0`,
   `1/r bm e_theta del_theta * (bm e_r del_r)` `= 1/r bm e_theta * (del_theta(bm e_r) del_r + bm e_r del_(theta r))` `= 1/r bm e_theta * bm e_theta del_r` `= 1/r del_r`,
   `1/r bm e_theta del_theta * (1/r bm e_theta del_theta)` `= 1/r^2 bm e_theta * (del_theta(bm e_theta) del_theta + bm e_theta del_theta^2)` `= 1/r^2 del_theta^2`. 四式相加即得结论.
2. 事实上, 还有两种计算极坐标下 Laplace 算子的方法. 一是分别计算 `del_x^2` 和 `del_y^2`: `del_x^2 = (cos theta del_r - 1/r sin theta del_theta)^2`
   `quad = cos^2 theta del_r^2` `+ 1/r^2 cos theta sin theta del_theta - 1/r cos theta sin theta del_(r theta)`
   `qquad + 1/r sin^2 theta del_r - 1/r sin theta cos theta del_(theta r)`
   `qquad + 1/r^2 sin theta cos theta del_theta + 1/r^2 sin^2 theta del_theta^2`. `del_y^2 = (sin theta del_r + 1/r cos theta del_theta)^2`
   `quad = sin^2 theta del_r^2` `- 1/r^2 sin theta cos theta del_theta + 1/r sin theta cos theta del_(r theta)`
   `qquad + 1/r cos^2 theta del_r + 1/r cos theta sin theta del_(theta r)`
   `qquad - 1/r^2 cos theta sin theta del_theta + 1/r^2 cos^2 theta del_theta^2`. 两式相加即得结论.
3. 还有一种是分别计算 `del_r^2` 和 `del_theta^2`, 注意到 `x, y` 对 `r, theta` 求导比较方便, `r, theta` 对 `x, y` 求导不方便, 因此先展开平方的后一项, 把 `theta` 的函数换到前面来以后, 才展开前一项: `del_r^2 = del_r (cos theta del_x + sin theta del_y)`
   `= cos theta del_r del_x + sin theta del_r del_y`
   `= cos theta (cos theta del_x + sin theta del_y) del_x` `+ sin theta (cos theta del_x + sin theta del_y) del_y`
   `= cos^2 theta del_x^2 + cos theta sin theta del_(y x)` `+ sin theta cos theta del_(x y) + sin^2 theta del_y^2`. `del_theta^2 = del_theta (-r sin theta del_x + r cos theta del_y)`
   `= -r del_theta(sin theta) del_x - r sin theta del_theta del_x` `+ r del_theta(cos theta) del_y + r cos theta del_theta del_y`
   `= -r cos theta del_x - r sin theta (-r sin theta del_x + r cos theta del_y) del_x`
   `quad - r sin theta del_y + r cos theta (-r sin theta del_x + r cos theta del_y) del_y`
   `= -r (cos theta del_x + sin theta del_y)`
   `quad + r^2 sin^2 theta del_x^2 - r^2 sin theta cos theta del_(y x)`
   `quad + r^2 cos^2 theta del_y^2 - r^2 cos theta sin theta del_(x y)`. 现在可以看出 `del_r^2 + 1/r del_r + 1/r^2 del_theta^2 = del_x^2 + del_y^2`.

1. 球坐标的 Laplace 算子. 利用 `bm e_r, bm e_theta, bm e_varphi` 的正交性以及 `del_varphi bm e_r = bm e_varphi`, `quad del_theta bm e_r = cos varphi bm e_theta`,
   `bm e_theta * (del_theta bm e_theta) = 0`,
   `del_varphi bm e_varphi = - bm e_r`, `quad del_theta bm e_varphi = -sin varphi bm e_theta`, 将下式展开: `nabla^2 = (bm e_r del_r + 1/(r cos varphi) bm e_theta del_theta + 1/r bm e_varphi del_varphi)^2`. 分别计算以下九项: `(bm e_r del_r)^2` `= bm e_r * bm e_r del_r^2` `= del_r^2`,
   `bm e_r del_r * (1/(r cos varphi) bm e_theta del_theta)` `= bm e_r * bm e_theta del_r (1/(r cos varphi) del_theta)` `= 0`,
   `bm e_r del_r * (1/r bm e_varphi del_varphi)` `= bm e_r * bm e_varphi del_r (1/r del_varphi)` `= 0`,
   `1/(r cos varphi) bm e_theta del_theta * (bm e_r del_r)` `= 1/(r cos varphi) bm e_theta * (cos varphi bm e_theta del_r + bm e_r del_(theta r))` `= 1/r del_r`,
   `(1/(r cos varphi) bm e_theta del_theta)^2` `= 1/(r^2 cos^2 varphi) bm e_theta * ((del_theta bm e_theta) del_theta + bm e_theta del_theta^2)` `= 1/(r^2 cos^2 varphi) del_theta^2`,
   `1/(r cos varphi) bm e_theta del_theta * (1/r bm e_varphi del_varphi)` `= 1/(r^2 cos varphi) bm e_theta * (-sin varphi bm e_theta del_varphi + bm e_varphi del_(theta varphi))` `= -1/r^2 tan varphi del_varphi`. `1/r bm e_varphi del_varphi * (bm e_r del_r)` `= 1/r bm e_varphi * (bm e_varphi del_r + bm e_r del_(varphi r))` `= 1/r del_r`,
   `1/r bm e_varphi del_varphi * (1/(r cos varphi) bm e_theta del_theta)` `= 1/r bm e_varphi * bm e_theta del_varphi (1/(r cos varphi) del_theta)` `= 0`,
   `(1/r bm e_varphi del_varphi)^2` `= 1/r^2 bm e_varphi * (-bm e_r del_varphi + bm e_varphi del_varphi^2)` `= 1/r^2 del_varphi^2`. 九项相加得到结论.

径向函数

设 `f(x_1, cdots, x_n)` 是 `n` 元数量函数, 如果 `f` 是球对称的, 即函数值只和到原点的距离有关, 而与方向无关, 就称 `f` 是径向函数. 记 `bm r = (x_1, cdots, x_n)`, `r = |bm r|`, 则 `f` 可以记为 `r` 的一元函数 `f(r)`.

设 `f(r)` 是 `RR^n` 上的径向函数, 注意到 `r` 本身也是径向函数, 有
1. `nabla * bm r = n`, `nabla r = bm r // r`, `nabla^2 r = (n-1)/r`;
2. `nabla f(r) = dd f r nabla r`, `nabla^2 f(r) = dd^2 f r + (n-1)/r dd f r`.

1. 只证最后一式. `nabla^2 r = nabla * nabla r` `= nabla * (bm r)/r`
   `= 1/r nabla * bm r + bm r * nabla 1/r`
   `= n/r - 1/r^2 bm r * nabla r`
   `= n/r - 1/r^3 bm r * bm r`
   `= (n-1)/r`.
2. 只证最后一式. `nabla^2 f(r)` `= nabla ("d"f)/("d"r) * nabla r + ("d"f)/("d"r) nabla^2 r`
   `= ("d"^2 f)/("d"r^2) nabla r * nabla r + ("d"f)/("d"r) nabla^2 r`
   `= ("d"^2 f)/("d"r^2) + (n-1)/r ("d"f)/("d"r)`. 这个结果与极坐标系、球坐标系下的 Laplace 算子相吻合.

微分方程 `f'' + (n-1)/r f' = 0` 的通解 `f(r) = { c_1 ln r + c_2, if m = 2; c_1 r^(2-n) + c_2, otherwise :}` 给出 `n` 维空间的径向调和函数. 特别 `ln r`, `1//r` 分别是二维和三维空间的径向调和函数. 这个例子告诉我们, 一个函数是否调和, 与它所在空间的维数有关.

多元积分公式

在二维空间中, 约定 `"d"V = dx dy`, 三维空间中则 `"d"V = dx dy dz`.

Green, Gauss, Stokes

Leibniz: `"d"/dt int_(alpha(t))^(beta(t)) f(x, t) dx` `= int_(alpha(t))^(beta(t)) del/(del t) f(x, t) dx` `+ f(beta(t), t) beta'(t)` `- f(alpha(t), t) alpha'(t).` Green (记忆: `del y` 者符号相反): `iint_V |del/(del x), del/(del y); P, Q| "d"V = oint_(del V) P dx + Q dy` Gauss: `iiint_V nabla * bm F "d"V` `= oiint_(del V) bm F * "d" bm S` Stokes: `iint_Sigma nabla xx bm F * "d" bm S = oint_(del Sigma) bm F * "d" bm r`,

用混合积的定义, Stokes 公式可以写为 `iint_S | del/(del x), del/(del y), del/(del z); P, Q, R; dy dz, dz dx, dx dy; | = oint_(del S) P dx + Q dy + R dz`, 取 `dz = 0` 就得到 Green 公式.

[来自 群Scalar] 求 `(t - t^3, 1 - t^4)`, `t in [-1, 1]` 围成图形的面积.

由 Green 公式 (注意 `t` 从 `-1` 到 `1` 曲线为顺时针) `iint_V dx dy` `= 1/2 oint_(del V) x dy - y dx` `= 1/2 int_1^(-1) [(t - t^3)(-4 t^3) - (1-t^4)(1-3t^2)] dt` `= 16/35`.

Green 三大公式

回忆方向导数的定义: `(del f)/(del bm n) = nabla f * bm n`, 在 Gauss 公式中取 `bm F = nabla f`, 则 `bm F * "d"bm S` `= nabla f * bm n |"d"bm S|` `= (del f)/(del bm n) "d"sigma`, 得到有用的方向导数形式: `iiint_V Delta f "d"V = oiint_(del V) (del f)/(del bm n) "d" sigma`, `quad bm n` 是单位外法向量. 在 Green 公式中取 `(Q, -P) = nabla f`, 也得到方向导数形式: `iint_V Delta f "d"V = oint_(del V) (del f)/(del bm n) "d"s`, `quad bm n` 是单位外法向量 `(dy/("d"s), -dx/("d"s))`. 这个公式实际是下面 Green 第一公式的特殊情形.

Green 第一公式 类比于分部积分公式将导数从一个因子转移到另一个因子上, 此公式将 nabla 算子从 `g` 转移到 `f` 上:
1. (2d) `iint_V nabla f * nabla g "d"V = int_(del V) g{::} (del f)/(del bm n) "d"s - iint_V g nabla^2 f "d"V`;
2. (3d) `iiint_V nabla f * nabla g "d"V = iint_(del V) g{::} (del f)/(del bm n) "d"sigma - iiint_V g nabla^2 f "d"V`.

只证 2d 情形 (3d 情形其实更简单). 利用方向导数的定义和 `bm n = (dy/("d"s), -dx/("d"s))`, `int_(del V) g{::} (del f)/(del bm n) "d"s` `= int_(del V) g nabla f * bm n "d"s` `= int_(del V) g ((del f)/(del x) dy - (del f)/(del y) dx)` `= iint_V [del/(del x)(g{::}(del f)/(del x)) + del/(del y)(g{::}(del f)/(del y))] "d"V` `= iint_V nabla f * nabla g "d"V + iint_V g nabla^2 f "d"V`.

Green 第二公式
1. (2d) `iint_V |nabla^2 f, nabla^2 g; f, g| "d"V = int_(del V) |(del f)/(del bm n), (del g)/(del bm n); f, g| "d"s`;
2. (3d) `iiint_V |nabla^2 f, nabla^2 g; f, g| "d"V = iint_(del V) |(del f)/(del bm n), (del g)/(del bm n); f, g| "d"sigma`.

分别对函数 `f, g` 应用 Green 第一公式即可. `nabla f * nabla g` 的积分刚好抵消.

Green 第三公式 设 `u` 为调和函数, `r` 是点 `(x,y)` 或 `(x,y,z)` 到 `del V` 上积分变动点的距离.
1. (2d) `u(x,y) = 1/(2pi) int_(del V) (u (del ln r)/(del bm n) - ln r (del u)/(del bm n)) "d"s`;
2. (3d) `u(x,y,z) = 1/(4pi) iint_(del V) (1/r (del u)/(del bm n) - u (del(1//r))/(del bm n)) "d"sigma`.

只证 2d 情形. 取 `C` 是以 `(x,y)` 为心, `rho` 为半径的圆周. 对函数 `u` 和 `ln r`, 在 `C` 和 `del V` 所夹的区域上应用 Green 第二公式, 注意 `ln r` 是二维空间的调和函数, 有 `int_(del V) - int_C = 0`, 其中被积函数是 `|(del ln r)/(del bm n), (del u)/(del bm n); ln r, u|`. 于是 `int_(del V)` `= int_C (u (del ln r)/(del bm n) - ln r (del u)/(del bm n)) "d"s` `= int_C (u (del ln r)/(del r) - ln rho (del u)/(del bm n)) "d"s` `= 1/rho int_C u "d"s`. 令 `rho to 0`, 则 `u` 在 `C` 上的平均值 `1/(2pi rho) int_C u "d"s` 趋于 `u(x,y)`, 从而 `u(x,y) = 1/(2pi) int_(del V)`.

Dirichlet 原理 在区域边界上取给定值的连续可微函数, 其 Dirichlet 积分 (函数的梯度的模的平方在区域上的积分) 取最小值当且仅当该函数为调和函数.

设 `nabla^2 f = 0`, 且 `f|{::}_(del V) = g|_(del V)`. 在 Green 第一公式中令 `g = f`, 利用边界条件得 `iint_V |nabla f|^2 "d"V` `= int_(del V) f{::} (del f)/(del bm n) "d"s` `= int_(del V) g{::} (del f)/(del bm n) "d"s` `= iint_V nabla f * nabla g "d"V`. 从而 `iint_V |nabla g|^2 "d"V - iint_V |nabla f|^2 "d"V` `= iint_V |nabla g|^2 "d"V + iint_V |nabla f|^2 "d"V - 2 iint_V nabla f * nabla g "d"V` `= iint_V |nabla g - nabla f|^2 "d"V ge 0`.

微分形式与外微分运算

外积 观察多元积分换元公式 `"d"u"d"v = (del(u, v))/(del(x, y)) dx dy`, 如果将 `x, y` 的次序对调, 得到 `dy dx = (del(y, x))/(del(x, y)) dx dy` `= |y_x, y_y; x_x, x_y| dx dy` `= |0, 1; 1, 0| dx dy` `= -dx dy`. 这个性质类似于向量的外积, 我们也把它称为外积, 写作 `dx ^^ dy`. 外积运算满足
1. 线性性. `(P dx + Q dy) ^^ "d"u = P dx ^^ "d"u + Q dy ^^ "d"u`;
2. 反对称性. `dx ^^ dy = - dy ^^ dx`.
由反对称性立即得到 `dx ^^ dx = 0`.
一般地, 在 `n` 维空间中, 规定每交换一对变量, 外积就改变一次符号: `dx_1 ^^ cdots ^^ color(red)(dx_i) ^^ cdots ^^ color(blue)(dx_j) ^^ cdots ^^ dx_n`
`= -dx_1 ^^ cdots ^^ color(blue)(dx_j) ^^ cdots ^^ color(red)(dx_i) ^^ cdots ^^ dx_n`.

1. 零阶微分形式 就是普通的多元函数 `f`.
2. 一阶微分形式 是 `dx, dy, dz` 以及它们的线性组合 (系数是普通函数, 如 `P dx + Q dy`).
3. 二阶微分形式 是 `dx ^^ dy`, `dy ^^ dz`, `dz ^^ dx` 以及它们的线性组合.
4. 三阶微分形式 是 `dx ^^ dy ^^ dz` 以及它的线性组合 (其实在三维空间中只有一种组合).
由 `dx ^^ dx = 0` 知道, 在 `n` 维空间中, 最高只有 `n` 阶微分形式.

设 `omega` 是微分形式, 定义外微分运算 `"d"omega` 如下:
1. 零阶微分形式 `f` 的外微分定义为它的全微分: `"d"f = nabla f * "d"bm r` `= sum (del f)/(del x_i) dx_i`.
2. 若 `omega` 形如 `omega = dx_1 cdots dx_n`, 则单项式 `f omega` 的外微分定义为 `"d"(f omega) = "d"f ^^ omega`. 利用外积运算的线性性, 可将 `"d"f` 展开各项与 `"d"omega` 相乘, 再合并同类项.
3. 最后, 规定外微分具有线性性 `"d"(omega_1 + omega_2) = "d"omega_1 + "d"omega_2`, 欲求整个微分形式的外微分, 只需将各个单项式的外微分相加.

1. 零阶微分形式 (普通函数) 求两次外微分结果为零: `"d"^2 f = 0`.

1. `"d"^2 f` `= "d"(sum_i pp f x_i "d"x_i)` `= sum_i "d"(pp f x_i) ^^ "d"x_i` `= sum_i sum_j pp^2 f (x_i x_j) "d"x_j ^^ "d"x_i`, 从而 `"d"^2 f` 展开为 `("d"x_1, cdots, "d"x_n)` 的二次型, 其系数矩阵反对称, 因此结果为零.

外微分下的三大公式
1. Green: `"d"(P dx + Q dy)` `= ((del Q)/(del x) - (del P)/(del y)) dx ^^ dy`.
2. Stokes: `"d"(P dx+Q dy+R dz)` `= ((del R)/(del y)-(del Q)/(del z)) dy ^^ dz` `+ ((del P)/(del z)-(del R)/(del x)) dz ^^ dx` `+ ((del Q)/(del x)-(del P)/(del y)) dx ^^ dy`.
3. Gauss: `"d"(P dy^^dz + Q dz^^dx + R dx^^dy)` `= ((del P)/(del x) +(del Q)/(del y) + (del R)/(del z)) dx ^^ dy ^^ dz`.
至此三大公式, 包括 Newton-Leibniz 公式都统一写成: `int_(del D) omega = int_D "d" omega`.

1. 左边等于 `"d"P ^^ dx + "d"Q ^^ dy` `= ((del P)/(del x) dx + (del P)/(del y) dy) ^^ dx` `+ ((del Q)/(del x) dx + (del Q)/(del y) dy) ^^ dy` 等于右边.
2. 和 Green 公式完全类似.
3. 第一项 `"d"P ^^ dy ^^ dz` `= ((del P)/(del x)dx+(del P)/(del y)dy+(del P)/(del z)dz) ^^ dy ^^ dz` `= (del P)/(del x) dx ^^ dy ^^ dz`. 其余两项是类似的 (注意 `dx ^^ dy ^^ dz = dy ^^ dz ^^ dx = dz ^^ dx ^^ dy`).