為何是這些方程，而非千百個別的
貫穿這個梯級，你學到了一門通用的手藝：餵給一個線性二階方程一個級數，硬磨出一條遞迴關係，再讓係數告訴你解。這門手藝對無窮多個方程都管用，而其中絕大多數既無名亦無聲。那麼，為何每本教科書都要停下來，點名其中寥寥幾個——貝索、勒讓德、厄米特、艾里——並用整整數章去厚待它們？因為這些特定的方程並非隨機而來：它們恰恰是當你試圖在圓形或球形的區域裡求解波動方程、熱方程或拉普拉斯方程時所掉落出來的東西。
其機制如下。當你用分離變數法去解一個偏微分方程時，那個多變數問題會碎裂成好幾個單變數的常微分方程——每個座標一個。在直角座標下，那些常微分方程就是熟悉的 y'' + lambda y = 0，其答案是正弦與餘弦。但改用極座標或球座標，座標因子 1/r 與 1/sin(theta) 便悄悄爬進方程，把常係數變成可變係數，並在原點或極點處安放一個奇點。徑向部分變成貝索方程；極角部分變成勒讓德方程。正弦與餘弦是平直幾何的特殊函數；貝索與勒讓德則是圓形幾何的特殊函數。
勒讓德：當級數客氣地停了下來
先從兩者中較友善的那個開始。勒讓德方程 (1 - x^2) y'' - 2x y' + n(n+1) y = 0 在 x = 0 處有一個常點，所以你可以從本梯級開頭那個直接的冪級數法去進攻它——在那裡並不需要弗羅貝尼烏斯法。代入 y = sum a_k x^k，收集同次冪，那條兩項遞迴便讀作 a_(k+2) = a_k · [k(k+1) - n(n+1)] / [(k+1)(k+2)]。注意那個分子：當 k 達到 n 時，括號 k(k+1) - n(n+1) 恰恰變為零。
那單獨一個零改變了一切。一旦 a_(k+2) = 0，那條鏈裡往後的每一個係數也都是零，於是無窮級數被鍘刀斬成一個 n 次的多項式。對每個非負整數 n，你都得到一個有限、性質良好的答案——這些就是**勒讓德多項式** P_n(x)：P_0 = 1，P_1 = x，P_2 = (3x^2 - 1)/2，如此這般。另一個獨立解（那條不終止的鏈）在 x = +/- 1 處爆掉，那正是球的兩極，所以物理學悄悄丟棄它，留下多項式。「解必須在兩極處保持有限」這項要求，正是一開始就逼著 n 必須是整數的原因。
勒讓德多項式有一份單純的正弦與餘弦不會誇耀的天賦：它們在區間 [-1, 1] 上彼此正交。只要 m 不等於 n，P_m(x) P_n(x) 在該區間上的積分就是零。這種正交性（這裡的權函數為 1）正是把球面上一個任意函數展開成勒讓德各塊之和背後的引擎——恰如傅立葉級數把一個週期函數展開成正弦與餘弦。每個 P_n 攫取函數的一份成分，而其餘各項則積分成零。
貝索：終於來了個道地的弗羅貝尼烏斯情形
貝索方程 x^2 y'' + x y' + (x^2 - p^2) y = 0 在 x = 0 處有一個正則奇點——尋常的冪級數在那裡失效，你必須動用兩篇之前所搭建的完整弗羅貝尼烏斯法。指標方程算出來是 r^2 - p^2 = 0，給出根 r = +p 與 r = -p。這正是指標方程當初為之而建的示範場合：這兩個根要嘛相差一個非整數（容易的情形，兩個乾淨的弗羅貝尼烏斯級數），要嘛相差一個整數（微妙的情形，此時第二個解可能拖進一個對數）。
取較大的根 r = +p，轉動弗羅貝尼烏斯的曲柄。遞迴逼使所有奇數係數消失，並把每個偶數係數綁到前一個上，於是第一類貝索函數 J_p(x) 就出來了。寫開來，J_p(x) = 對 m 求和 [ (-1)^m / (m! · Gamma(m + p + 1)) ] · (x/2)^(2m + p)。分母裡的階乘以兇猛的速度增長，所以收斂半徑是無窮大——J_p 是一個老老實實、近乎整函數的東西，像正弦那樣振盪，但振幅隨 1/sqrt(x) 緩緩衰減，且其零點並非等距分布。
那第二個解呢？當 p 不是整數時，較小的根 r = -p 交出一個完好的獨立級數 J_(-p)(x)，通解不過是這兩者的組合。但當 p 是整數時，這兩個根相差一個整數，那個容易的第二級數便塌陷成第一個的一個倍數——正是指標方程那一篇曾警告你的那個障礙。弗羅貝尼烏斯法於是逼著一個對數進入第二個解，喚作第二類貝索函數 Y_p(x)，它在 x = 0 處無界。完整的通解是 y = c1 J_p(x) + c2 Y_p(x)，而在一個圓心必須保持有限的物理圓盤裡，你留下 J_p，把 Y_p 扔掉。
一座小小的動物園，以及牠們的共通點
貝索與勒讓德也有表親，每一個都是各自那個可變係數方程的有名解，每一個都是從物理學裡拉出來的。以下是這個家族的一覽。
function equation / setting shows up in -------------------------------------------------------------------- Legendre P_n (1-x^2)y'' - 2x y' + n(n+1)y=0 sphere, gravity/electrostatic potential Bessel J_p x^2 y'' + x y' + (x^2-p^2)y = 0 drumhead, cylinder, wave in a pipe Hermite H_n y'' - 2x y' + 2n y = 0 quantum harmonic oscillator Laguerre L_n x y'' + (1-x)y' + n y = 0 hydrogen atom radial part Chebyshev T_n (1-x^2)y'' - x y' + n^2 y = 0 approximation, numerical methods Airy Ai, Bi y'' - x y = 0 turning points, optics caustics
一旦你看出來，這個規律便十分醒目。厄米特、拉蓋爾與切比雪夫都共享勒讓德的好運氣：對整數 n，它們的級數終止成多項式，而且它們都正交——各有自己的權函數，厄米特用 e^(-x^2)，拉蓋爾用 e^(-x)，切比雪夫用 1/sqrt(1-x^2)。相對地，艾里函數從不終止；它們是完整的冪級數，在 x = 0 的一側像衰減的波、另一側像增長的指數，這也正是為何它們主宰著轉折點處的過渡。而這當中許多，都是一個宏大母方程——超幾何方程——的特例，超幾何級數在其參數取特定值時，便涵蓋了勒讓德、切比雪夫以及更多。
正交性才是它們重要的更深原因
人們很容易以為這些函數的重點是寫下閉式。更深一層的回報，是你在勒讓德那裡瞥見的正交性，而這絕非偶然。這些方程每一個都能被揉成一個共同的形狀——一個史特姆－劉維爾問題——而有一條定理保證：屬於不同特徵值的解，關於某個特定的權函數彼此正交。正是這種正交性，讓你能把一個一般的函數展開成特殊函數的級數，並一次一個積分地解出係數，背後是與傅立葉級數同樣的分而治之邏輯。
於是這條觀念之鏈乾淨俐落地閉合了。一個在圓形座標下的偏微分方程分離成有名的常微分方程；弗羅貝尼烏斯法與冪級數法解出那些常微分方程，產出貝索、勒讓德及其表親；有限性的物理要求挑出性質良好的那一支，並把階 n 量子化；而正交性接著讓你把完整的偏微分方程解重新組裝成這些塊的一個級數，逐係數地去配合邊界與初始資料。本梯級的每一件工具，都在那條鏈裡做著真正的活兒。
要帶往前去的一個誠實提醒。為一個函數命名，並不表示你能徒手把它算出來：J_2(3.7) 或 P_5(0.41) 來自數值表、遞迴，或一次函式庫的呼叫，絕非來自一條由初等函數構成的整齊公式。而這些有名的函數縱然光輝，所涵蓋的也只是那一小撮恰好能分離標準幾何的方程——絕大多數可變係數的常微分方程根本沒有特殊函數的答案，必須交給後幾個梯級的數值與定性方法。特殊函數是一座美麗而有限的可解之島，浮在一片大得多的海上。