普通冪級數失效的地方
在上一篇裡，你充滿信心地用了[[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 被邊界條件丟棄了。這個教訓遠不止於貝塞爾。在同一個模板上換上不同的係數運行弗羅貝尼烏斯，你就一個接一個地造出物理學裡那些有名的函數。這個方法不是只對一個方程管用的把戲；它是一座熔爐，整整一冊特殊函數都在其中被錘打成形，每一個都是它自己那個奇異方程的弗羅貝尼烏斯解。