普通幂级数失效的地方
在上一篇里，你充满信心地用了[[power-series-method|幂级数法]]：设 y = sum of a_n x^n，把它代入方程，再磨出一个系数的递推关系。这做得很漂亮——但只在[[ordinary-and-singular-points|常点]]处成立，那里方程的系数函数都规规矩矩。麻烦在于：物理学中最重要的那些方程，偏偏在我们最想要展开的那个点（通常是 x = 0）处出问题。在那里尝试普通幂级数，要么它崩溃，要么它根本无法表示真正的解。
看一个最简单的警示例子，柯西–欧拉方程 x^2 y'' - 2y = 0。你或许还记得，柯西–欧拉方程的解是纯幂次：y = x^2 与 y = x^{-1}。其中 x^2 是一个完全普通的幂级数，但 x^{-1} 不是——任何 a_0 + a_1 x + a_2 x^2 + ... 之和都绝不可能生出一个 1/x。普通幂级数从结构上就看不见这第二个解。方程在告诉我们：解想要从一个负数次幂开始，而我们的拟设却禁止了它。
于是修补办法几乎是自己写出来的。如果解坚持要从 x 的某个并非非负整数的次幂开始，那就一开始就把那个次幂内置进拟设里。给一个正常的幂级数乘上一个首因子 x^r，其中 r 是一个我们将去求解的未知指数。这一个想法——让级数从一个可调的次幂 x^r 起步——就是[[frobenius-method|弗罗贝尼乌斯方法]]的全部。
正则奇点：可以修复的损坏
并非每个奇点都能挽救，所以我们必须诚实地划出分界线。把方程写成标准形 y'' + P(x) y' + Q(x) y = 0。一个点 x_0 是[[regular-singular-point|正则奇点]]，如果 P 或 Q 确实在那里发散，但只是温和地发散：(x - x_0) P(x) 与 (x - x_0)^2 Q(x) 都必须在 x_0 处保持有限——实际上是解析的。用文字说，就是 P 发散得不能比 1/(x - x_0) 更猛，Q 发散得不能比 1/(x - x_0)^2 更猛。如果奇性比这更凶，该点就是非正则的，弗罗贝尼乌斯方法不提供任何保证。
有了一个正则奇点，弗罗贝尼乌斯拟设就是 y = x^r 乘 (sum of a_n x^n for n from 0 to infinity)，按惯例 a_0 不为零——首项系数必须非零，否则我们只是把 r 重新贴了个标签。逐项求导，就像你在第一卷对普通幂级数做的那样，只是现在每一项都多带一个指数 r：第 n 项 a_n x^{n+r} 的导数是 (n+r) a_n x^{n+r-1}，二阶导数是 (n+r)(n+r-1) a_n x^{n+r-2}。指数 r 一路搭车穿过每一次求导，而这正是它能从最低次幂的平衡中掉出来的原因。
指标方程定下首项指数
关键的一步来了。把弗罗贝尼乌斯级数代入方程后，按 x 的次幂归并各项，并要求每一个系数都为零——这无非又是级数法那条规则：一个幂级数为零当且仅当它所有系数为零。现在看出现的那个唯一的最低次幂，依你怎么排列，它是 x^{r-2} 项或 x^r 项。它完全由 a_0 搭成，而既然 a_0 不许为零，乘在它前面的那个括号本身就必须为零。那个括号置零，就是[[indicial-equation|指标方程]]。
指标方程永远是关于 r 的一个简单二次方程，形如 r(r-1) + p_0 r + q_0 = 0，其中 p_0 与 q_0 是被驯服的组合 x P(x) 与 x^2 Q(x) 在 x_0 处的取值。它是你为常系数 ODE 解过的特征方程在奇点处的表亲：在那里，一个关于增长率的二次方程挑出指数函数；在这里，一个关于 r 的二次方程挑出首项幂次。它的两个根 r_1 与 r_2（标记使 r_1 大于等于 r_2）称为该奇点的指标，它们是解可以如何起步的仅有候选。
y = x^r ( a_0 + a_1 x + a_2 x^2 + ... ), a_0 =/= 0
Standard form: y'' + P y' + Q y = 0
Regular singular point at 0: p0 = lim_{x->0} x P(x)
q0 = lim_{x->0} x^2 Q(x)
Lowest power of x -> only a_0 survives ->
INDICIAL EQUATION: r(r-1) + p0 r + q0 = 0
Roots r1 >= r2 = the two exponents of the singularity.
Higher powers give the recurrence that fixes a_1, a_2, ...第二个解的三副面孔
较大的根总是干净利落地奏效。把 r = r_1 代入更高次幂的方程，得到一个递推关系，依次定出 a_1、a_2、……，于是你得到一个诚实的弗罗贝尼乌斯解 y_1 = x^{r_1} 乘一个幂级数。所有的戏剧性都集中在第二个解上，而它取三种形式中的哪一种，由差 r_1 - r_2 决定。那一个数就是全部的预报。
- 情形 1——两根之差不是整数（r_1 - r_2 不是整数）。此时两个根都好使，你得到两条干净独立的弗罗贝尼乌斯级数：y_1 = x^{r_1}(...) 与 y_2 = x^{r_2}(...)。通解就是它们的组合；大功告成。（相等的根算作整数 0，因此归入情形 2，而非这里。）
- 情形 2——两根相等（r_1 = r_2）。一条弗罗贝尼乌斯级数根本无法提供第二个独立解——可供起步的指数只有一个。缺失的搭档被迫带上一个对数：y_2 = y_1 ln(x) + x^{r_1}(一条新级数)。这个对数不是装饰；它是重指数不可避免的标志，正如在常系数世界里重特征根逼出一个额外的 x 因子那样。
- 情形 3——两根相差一个正整数（r_1 - r_2 = N，一个整数）。这是最棘手的情形。当你试图从较小的根 r_2 出发构造级数时，递推关系会在第 N 项恰好撞上一个要除以零的步骤——而那正是较大根的级数起步之处。有时这个零是宽容的，你仍得到一条纯级数；有时它像情形 2 那样逼出一个对数：y_2 = c y_1 ln(x) + x^{r_2}(一条级数)，其中常数 c 有可能恰好为零。
贝塞尔方程：方法的现身说法
让整台机器在那个几乎催生了本方法的方程上运转：贝塞尔方程 x^2 y'' + x y' + (x^2 - nu^2) y = 0。化为标准形 P(x) = 1/x、Q(x) = (x^2 - nu^2)/x^2，于是 x P = 1、x^2 Q = x^2 - nu^2，二者在零点都有限：x = 0 是教科书级的正则奇点。指标方程为 r(r-1) + 1 乘 r - nu^2 = 0，化简成 r^2 = nu^2。于是两个指数是 r_1 = nu 与 r_2 = -nu——而它们的差 2 nu，恰恰是选定我们落入三种情形中哪一种的那个旋钮。
看着三种情形仅靠改变 nu 就从同一个方程里全部冒出来。若 nu = 1/2，差 2 nu = 1 是整数（情形 3），可奇性是宽容的，解竟是初等的 sin(x)/sqrt(x) 与 cos(x)/sqrt(x)——纯级数，无对数。若 nu 取比如 1/3，差 2/3 不是整数（情形 1），你得到两条干净独立的级数，即 ±1/3 阶的贝塞尔函数。但若 nu = 0，两根重合（情形 2），不可避免的对数就现身了：第二个解，即第二类贝塞尔函数 Y_0(x)，带着一个 ln(x) 项，并在原点处俯冲向负无穷。
这正是为什么振动的圆形鼓面在中心处有限、而点热源却不有限：物理上允许的解是那个保持有界的解，即第一类贝塞尔函数，而带对数的 Y_0 被边界条件丢弃了。这个教训远不止于贝塞尔。在同一个模板上换上不同的系数运行弗罗贝尼乌斯，你就一个接一个地造出物理学里那些有名的函数。这个方法不是只对一个方程管用的把戏；它是一座熔炉，整整一册特殊函数都在其中被锤打成形，每一个都是它自己那个奇异方程的弗罗贝尼乌斯解。