这些函数从何而来
上一篇里，伽马函数为阶乘在每一个实数与复数自变量上都安了家。这正是本级阶梯的统一套路：物理或统计中的某个问题反复写出同一个定积分，而这个积分没有初等闭形式，于是我们给它起个名字，一劳永逸地研究它。本篇的四个函数正是如此——它们是一批你会反复遇到的积分的、有名字、可制表、可画图的答案。
贝塔函数：伽马函数的双臂兄弟
贝塔函数由一个在单位区间上的积分定义：B(x, y) = 从 0 到 1 对 t^{x-1} (1-t)^{y-1} dt 的积分，在 x > 0 且 y > 0 时有效。可以把它想成两个幂之间的一次加权握手——一个因子偏向左端点 t = 0，另一个偏向右端点 t = 1。它是对称的，B(x, y) = B(y, x)，原因很简单：把 t 换成 1 - t 就把两个端点对调了。
最醒目的事实是：贝塔其实根本不是一个新函数：B(x, y) = Gamma(x) Gamma(y) / Gamma(x+y)。于是一张伽马函数表就为每一个贝塔值定了价。这个恒等式正是贝塔无处不在的原因——只要出现形如 t^{a}(1-t)^{b} 的积分（在概率论和组合学里它到处都是），你就能把答案读作一组伽马值之比。做个验算：B(1, 1) = 从 0 到 1 对 1 dt 的积分 = 1，而 Gamma(1)Gamma(1)/Gamma(2) = 1·1/1 = 1，吻合。
误差函数：驯服钟形曲线
在第一卷你已认识高斯积分：从负无穷到正无穷对 e^{-x^2} dx 的积分 = sqrt(pi)，从 0 到无穷对 e^{-x^2} dx 的积分 = sqrt(pi)/2。这个干净的总值藏着一桩烦恼：部分面积，即从 0 到 z 的积分，没有初等原函数。对于‘钟形曲线在 0 与 z 之间占了多少面积’，并不存在用指数和幂写出的公式。于是我们给它起名。误差函数为 erf(z) = (2/sqrt(pi)) · 从 0 到 z 对 e^{-t^2} dt 的积分，常数 2/sqrt(pi) 是这样选的：使 erf(无穷) = 1。
想象它的图象：erf 是一条奇对称的 S 形曲线，erf(0) = 0，先迅速上升，再逐渐压平、贴向直线 y = 1（左侧则贴向 y = -1）。它的搭档，互补误差函数 erfc(z) = 1 - erf(z)，度量 erf 留下的那条尾巴——z 之外的面积，也就是一个正态样本落在远离均值处的概率。统计学家就住在这条尾巴里，所以 erfc 也配得上自己的名字与自己的按钮。
我们究竟怎样取出数值？在零附近，把被积函数展成泰勒级数——e^{-t^2} = 1 - t^2 + t^4/2 - ...——再逐项积分，它处处收敛，给出 erf(z) = (2/sqrt(pi))(z - z^3/3 + z^5/10 - ...)。对大的 z 这条级数没用（需要巨量相消），此时改用尾部的渐近展开：erfc(z) ~ e^{-z^2}/(z sqrt(pi)) · (1 - 1/(2z^2) + ...)。这条渐近级数若把所有项都保留其实是发散的——可一旦在大 z 处只取前几项截断，它精确得令人咋舌。这种‘小 z / 大 z’的二分法，是特殊函数计算反复出现的节奏。
指数积分、正弦积分与余弦积分
还有三个‘弄丢了闭形式’的积分函数补全这一家族，每当一个看似无害的被积函数在分母上带个 1/t 或 t，你就会遇到它们。指数积分为 Ei(x) = 从负无穷到 x 对 e^{t}/t dt 的积分（在 t = 0 处的极点取主值），以及与之密切相关的 E_1(x) = 从 x 到无穷对 e^{-t}/t dt 的积分。正弦积分与余弦积分为 Si(x) = 从 0 到 x 对 sin(t)/t dt 的积分，以及 Ci(x) = 负的、从 x 到无穷对 cos(t)/t dt 的积分。
Si 最容易画出来。它的被积函数 sin(t)/t 就是那著名的鼓包——在 t = 0 处等于 1，再以越来越小的波纹荡下去。当 x 向外延伸，Si(x) 上升、过冲，随后以衰减振荡安顿到极限 Si(无穷) = pi/2——这个数恰好就是狄利克雷积分，从 0 到无穷对 sin(t)/t dt 的积分 = pi/2。Si 的过冲不是数值故障；它正是本卷稍后会成为傅里叶级数跳跃处吉布斯现象的那种振铃。相比之下，Ci 在 x 趋于 0 时像 ln(x) 一样跌向负无穷，再在大 x 处向 0 振荡。
它们都在哪儿出现？Ei 和 E_1 支配辐射输运、线源散热，以及那著名的、估计某数以下有多少素数的对数积分。Si 和 Ci 刻画狭缝边缘的衍射图样，以及滤波器和天线的响应。它们一点也不冷僻——它们正是物理与工程里那种日常的‘本该是闭形式却偏偏不是’，这恰恰就是它们被命名的原因。
它们如何相连——以及如何使用
这些函数并非五件互不相干的奇珍；它们是同一棵树上的叶子。贝塔是若干伽马函数的闭形式组合。误差函数则是上限可移动的伽马——事实上 erf 正是同一个高斯积分的不完全版本，而那个积分一旦补全，就给出了上一篇的 Gamma(1/2) = sqrt(pi)。而 erf、Ei、Si、Ci 用的全是同一招：取一个原函数逃出初等函数范围的被积函数，把流动的面积定义成一个新函数，再去研究它。认出这一招，才是本级真正要练的本领。
- 认出标志性被积函数：[0,1] 上的 t^{a}(1-t)^{b} 意味着贝塔；e^{-t^2} 意味着误差函数；e^{t}/t 意味着指数积分；sin(t)/t 或 cos(t)/t 意味着正弦或余弦积分。
- 把你的积分揉成标准形：用换元做缩放与平移，让积分限和被积函数与定义严丝合缝——大部分功夫不过是常数的记账。
- 读出数值：对贝塔，转成伽马值；其余则查表或调用 erf、erfc、Ei、Si、Ci，在 0 附近用小自变量级数、在大自变量处用渐近尾。
书架上紧邻的是椭圆积分——下一篇的主题——它由椭圆的弧长和真实单摆的摆动逼出来。教训一字不变地延续：当真实问题变硬时，答案往往不是‘积不出来’的失败，而是‘还没人告诉你它名字’的一个函数的登场。