本指南要解開的謎
在本階梯裡，你已見過三大線性偏微分方程——熱方程、波方程，以及定常的拉普拉斯方程 nabla^2 u = 0——你也在矩形上用分離變量法馴服了它們，得到的是普通的正弦、餘弦模式和一個傅里葉級數。但世界很少是矩形的。鼓面是一張圓盤，受熱的導線是一根圓柱，行星的引力場活在一個球面之外。誠實的問題是：當區域變圓時，同一套方法會發生什麼？
這裡有一句可以貫穿整篇指南的點睛之語。分離變量法依然有效——但這只因為我們先把拉普拉斯算子改寫到貼合邊界的座標裡；一旦這樣做，剩下的單變量方程就不再是那個溫和的 y'' + lambda y = 0 了。它們正是你在級數解階梯和特殊函數階梯裡學過的貝塞爾方程與勒讓德方程。圓形的幾何並非僅僅把代數弄複雜；它伸進你的工具箱，精確地挑出哪些特殊函數獲准住在那個形狀上。讀到最後你會明白，為什麼圓盤逼出貝塞爾函數、球面逼出勒讓德多項式，而這種選擇會讓你覺得是必然的，而非死記硬背的。
當空間彎曲，拉普拉斯算子變成什麼
首先，從直覺上看，拉普拉斯算子究竟是什麼？在笛卡爾座標裡 nabla^2 u = d^2u/dx^2 + d^2u/dy^2 + d^2u/dz^2，是若干二階偏導之和。但那個公式只是一個無關座標的想法的局部面孔：nabla^2 u 度量的是某點處的 u 與它周圍一個小球面上 u 的平均值相差多少。若某點比四周更涼，那裡的 nabla^2 u 就為正，熱量會流入；而拉普拉斯方程 nabla^2 u = 0 恰恰是說「每一點都已經等於它鄰居的平均」——一個徹底鬆弛下來的調和場，再沒有熱點或冷點需要抹平。這幅取平均的圖景在任何座標系裡都一樣；改變的只是記賬方式。
為什麼記賬方式非改不可？因為在曲線座標裡，網格線並非均勻分佈。在圓盤中心附近，把角度 theta 撥過一格，你幾乎沒動；在邊緣附近撥過同樣一格，你卻掃出一道長長的弧。這些把每個座標撥動一格對應多少真實距離記錄下來的換算因子，就是度規因子（scale factors），而它們正是全部關鍵。對極座標 (r, theta)，度規因子是 h_r = 1、h_theta = r：沿 r 移動是誠實的距離，但沿 theta 撥一格掃過的弧長是 r·d(theta)，r 越大弧越長。拉普拉斯算子必須除以這些因子，以抵消網格的拉伸。
圓盤：那些 1/r 項如何孕育出貝塞爾方程
取最平、卻仍有趣的情形：寫在極座標裡的拉普拉斯方程（或亥姆霍茲方程），定義在一張圓盤上。把度規因子 h_r = 1、h_theta = r 代入配方，得到 nabla^2 u = d^2u/dr^2 +（1/r）du/dr +（1/r^2）d^2u/d(theta)^2。前面帶 1/r 和 1/r^2 的那兩塊，正是曲率的指紋；它們在平直的笛卡爾形式裡並不存在，而它們即將包辦全部的活。它們不是錯誤也不是近似——它們是角度 theta 越往外覆蓋越多地盤所付出的、精確的代價。
現在分離：猜一個乘積 u(r, theta) = R(r) Theta(theta)，代入，再整體除以 R Theta。角度與半徑乾淨地分到等號兩側，這只有在兩側都等於同一個常數時才可能——把它叫做 m^2。角度那一半給出 Theta'' + m^2 Theta = 0，跟矩形上一樣友好的振子，解為 cos(m theta) 與 sin(m theta)。但這裡幾何提出了矩形從不曾提的要求：繞圓盤一圈，從 theta 到 theta + 2pi，會回到同一個物理點，所以 Theta 必須以 2pi 為週期重複。這種單值性強迫 m 取整數 0, 1, 2, 3, ...——這個整數是圓的拓撲遞給你的，而非手工選定的。
徑向那一半，正是那個有名字的函數誕生之處。把角度導數換成 -m^2，並（對振動的鼓面）帶上一個頻率常數 k^2 之後，半徑滿足 R'' +（1/r）R' +（k^2 - m^2/r^2）R = 0。令 x = k r，這就是一字不差的 m 階貝塞爾方程：x^2 R'' + x R' +（x^2 - m^2）R = 0。你並沒有從外面搬來貝塞爾方程；是曲線拉普拉斯算子留下的 1/r 與 1/r^2，加上圓所供給的整數 m，當著你的面把它拼了出來。它的有界解是貝塞爾函數 J_m(k r)，而要求鼓面在邊緣被夾緊——R(a) = 0——便逼得 k a 必須是 J_m 的一個零點，這正是鼓的允許頻率（它的泛音）被量子化的方式。
Laplacian on a disk (polar): nabla^2 u = u_rr + (1/r) u_r + (1/r^2) u_thetatheta
Separate u(r,theta) = R(r) Theta(theta) , divide by R Theta:
R''/R + (1/r) R'/R + (1/r^2) Theta''/Theta = -k^2
|__________ depends on r only _________| |__ theta only __|
=> Theta'' + m^2 Theta = 0 (period-2pi => m = 0,1,2,...)
=> R'' + (1/r) R' + (k^2 - m^2/r^2) R = 0
Let x = k r : x^2 R'' + x R' + (x^2 - m^2) R = 0 <-- BESSEL of order m
bounded solution R = J_m(k r) ; clamp R(a)=0 => k a = zero of J_m球面：cos(theta) 如何孕育出勒讓德方程
登到三維，選球座標 (r, theta, phi)，其中 theta 是從北極量起的極角，phi 是經度。度規因子是 h_r = 1、h_theta = r、h_phi = r sin(theta)，把它們送進配方，得到一個含三個嵌套部分的拉普拉斯算子：一個徑向部分、一個帶 1/sin(theta) 因子的極角部分，以及一個方位項 1/(r^2 sin^2(theta)) d^2u/d(phi)^2。那個孤零零的 sin(theta) 是曲率的新指紋——它是經線向兩極擠攏時的那種壓縮，正是這種壓縮讓平面世界地圖把格陵蘭拉得老大。
分三步分離：u = R(r) Theta(theta) Phi(phi)。經度因子最先分出，給出 Phi'' + m^2 Phi = 0，又一次，2pi 週期性逼著 m 取整數——和圓盤相同的圓形拓撲，如今作用在經度上。半徑因子接著分出，給出一個等維方程 r^2 R'' + 2r R' - l(l+1) R = 0，其解是乾淨的冪 r^l 與 r^{-(l+1)}——那些隨距離的冪次增長或衰減的熟悉場，後者在 l = 0 時恰是點電荷或質點的 1/r 勢。
極角，正是第二個有名字的函數登場之處。Theta 方程仍帶著那個彆扭的 1/sin(theta)，但一個靈光的替換就把它清掉：令 x = cos(theta)，於是從極到極、theta 從 0 到 pi 的掃掠，變成 x 走過整潔的區間 -1 到 1。在這一變換下，那個雜亂的三角算子坍縮成 d/dx[(1 - x^2) dTheta/dx] + [l(l+1) - m^2/(1 - x^2)] Theta = 0——即連帶勒讓德方程。當不依賴經度（m = 0）時，它就是普通的勒讓德方程 (1 - x^2) y'' - 2x y' + l(l+1) y = 0。幾何再一次替你選定了函數：要求解在兩極 x = +1 與 x = -1 處保持有限，便逼得 l 取非負整數，而能存活下來的唯一解，正是勒讓德多項式 P_l(cos theta)。
為什麼這非如此不可：更深的規律
退後一步，每個座標系裡的規律都是相同的。分離把多變量偏微分方程一層層剝成每個座標一個常微分方程，而這些單變量方程中的每一個都是一個斯圖姆–劉維爾問題——一個區間上的自伴本徵值問題，正是佔據了本階梯整整一級的那種結構。在矩形上，區間平平無奇，斯圖姆–劉維爾問題是平凡的 y'' + lambda y = 0，所以本徵函數是正弦與餘弦。在圓盤和球面上，區間自帶曲線座標的權重，同一台斯圖姆–劉維爾機器吐出的便改成貝塞爾函數與勒讓德多項式。幾何不同，機器相同，吐出的有名字的本徵函數也不同。
這條斯圖姆–劉維爾血統並非腳註——它正是讓整套方法可用的關鍵。因為每個徑向或角向問題都是自伴的，它的本徵函數正交且完備，於是你能把任何合理的邊界數據展開成 J_m 的貝塞爾–傅里葉級數，或 P_l 的勒讓德級數，正如你在矩形上用普通傅里葉級數展開那樣。配方在精神上是完全一樣的：分離、找出本徵函數，再藉助正交性、用匹配邊界條件來定下剩餘的係數。特殊函數不過是適配圓形邊界的正確字母表，而正交性就是讓你能用那套字母表讀出任意函數的詞典。
- 選與邊界相符的座標——圓而直的形狀用極座標/柱座標，球狀的用球座標——使邊界變成某個座標取定值（比如 r = a）。
- 用度規因子把 nabla^2 改寫到那些座標裡，再代入乘積 u = R(r) Theta(theta) Phi(phi) 並整體相除，把變量分開。
- 先解容易的週期性角向部分；它們的單值性把分離常數 (m, l) 量子化為整數，再餵入徑向/極角方程。
- 認出剩下的方程是貝塞爾或勒讓德方程，只保留在該有限處保持有限的那個解，再把答案拼成一個正交級數，其係數由邊界數據定下。
誠實的告誡，以及魔法在何處止步
要看清楚是什麼讓這一切奏效，因為它很脆弱。分離變量法不是通用求解器；只有當幾何與方程彼此配合時它才成功。邊界必須是某個座標的等值面（一個圓 r = a，一個球面 r = a），而座標系必須是那少數幾個走運的之一——三維拉普拉斯算子能分離的座標系恰好只有十一個，不能再多。把區域歪成一個橢圓，或一個偏心的圓環，變量便拒絕分開；乾淨的乘積 u = R(r)Theta(theta) 那時乾脆就是錯的，你必須訴諸數值方法，或一種全然不同的攻法。那些有名字的函數是對稱性的饋贈，且只是對稱性的饋贈。
再說兩點誠實話。第一，每個徑向方程都有兩個獨立解，而我們只保留一個——保留在區域所及之處不會炸掉的那個。在實心圓盤上我們捨棄第二類貝塞爾函數 Y_m，因為它在中心 r = 0 處發散；在球上當區域包含原點時我們丟掉 r^{-(l+1)}。這種捨棄是一個物理判斷（真實的鼓面在正中並非無窮大），而非數學上的條件反射——而對外部問題它會翻轉：那時你反過來保留衰減的解，扔掉在無窮遠處增長的那個。永遠要問：你實際的區域能容忍哪一個解。