這個計算器的用途
本工具可求出實階數 \(v\) 之下，第一類貝塞爾函數 \(J_{v}(x)\) 與第二類貝塞爾函數 \(Y_{v}(x)\)（又稱諾依曼函數或韋伯函數）的第 \(s\) 個正零點。這些零點通常記為 \(j_{v,s}\) 與 \(y_{v,s}\)，廣泛出現在物理與工程領域：圓形振動膜（鼓面）、圓柱體熱傳導、電磁波導，以及傅立葉—貝塞爾級數等。這是一個純數學的特殊函數工具，全球通用，不受地區或單位的限制。
使用方式
輸入階數 \(v\)（實數，常見範圍為 0 到 200）與零點序號 \(s\)（正整數 1、2、3……）。計算器會回傳 \(j_{v,s}\)，也就是 \(J_{v}(x)\) 的第 \(s\) 個正根；以及 \(y_{v,s}\)，即 \(Y_{v}(x)\) 的第 \(s\) 個正根。舉例來說，\(v = 0\)、\(s = 1\) 對應的就是鼓面振動的基本模態。
公式與計算方法
這兩個函數都是貝塞爾方程式 \(x^2 y'' + x y' + (x^2 - v^2) y = 0\) 的解。\(J_{v}\) 由其收斂的冪級數展開求值；\(Y_{v}\) 則採用 $$Y_{v} = \frac{J_{v}\cos(v\pi) - J_{-v}}{\sin(v\pi)}$$ 計算，整數階的情況則以數值極限方式處理。為了定位第 \(s\) 個零點，我們先以 McMahon 漸近估計作為初始值，框出零點所在的區間，再以二分法反覆逼近直到收斂。由於 \(Y_{v}(x)\) 在 \(x = 0\) 處會趨向 \(-\infty\)，搜尋會從一個很小的正數 \(x\) 開始，以避開對數奇異點。
實例演算
當 \(v = 0\)、\(s = 1\) 時：\(J_{0}(x)\) 的第一個正零點為 \(2.4048255577\)，\(Y_{0}(x)\) 的第一個正零點為 \(0.8935769663\)。當 \(v = 1\)、\(s = 1\) 時：$$j_{1,1} = 3.8317059702, \quad y_{1,1} = 2.1971413260$$
零点在應用中的含義
貝塞爾函數零點並非抽象的數學奇想——當在圓形或圓柱形區域上提出波動、熱傳導或勢能問題時，它們是浮現出來的離散特徵值。邊界條件強制解的徑向部分在邊界處消失（或其導數消失），而這一條件僅在零點 \(j_{v,s}\) 處滿足。
振動圓形膜（鼓面）
對於邊緣夾緊的理想圓形鼓，位移分離為由 \((v,s)\) 標記的振動模式，其中 \(v\) 計數角向節點直徑的數量，\(s\) 計數徑向節點圓的數量。允許的特徵頻率為 \(f_{v,s}=\frac{c}{2\pi a}\,j_{v,s}\)，其中 \(c\) 是波速。基本音使用 \(j_{0,1}=2.4048256\)；較高的 \(s\) 和較高的 \(v\) 都會提高音高，由於 \(j_{v,s}\) 不是彼此的整數倍，鼓的泛音是不和諧的。
圓柱形熱傳導
當在表面溫度固定的長圓柱形區域上求解熱方程時，徑向特徵函數為 \(J_0(\lambda_s r/a)\)，其中 \(\lambda_s=j_{0,s}\)。每個模式在時間上衰減為 \(\exp\!\left(-\alpha (j_{0,s}/a)^2 t\right)\)，因此最小零點 \(j_{0,1}\) 控制衰減最慢、壽命最長的溫度分佈。較大的 \(s\) 給出較大的特徵值，因此衰減更快。
波導截止頻率
在半徑為 \(a\) 的空心圓形金屬波導中，橫向磁場（TM）模式在由 \(j_{v,s}\) 設定的頻率處截止，而橫向電場（TE）模式在導數 \(J_v'\) 的零點處截止。對於 TM 模式，截止頻率為 \(f_{c}=\frac{c\,j_{v,s}}{2\pi a}\)；只有超過此頻率時模式才能傳播。最低 TM 模式（TM\(_{01}\)）再次使用 \(j_{0,1}\)。
傅立葉–貝塞爾級數
圓盤上的任何合理函數都可以展開為 \(f(r)=\sum_{s=1}^{\infty} c_s\,J_v\!\left(j_{v,s}\,r/a\right)\)。縮放後的零點 \(j_{v,s}/a\) 的作用完全像普通傅立葉正弦級數的波數 \(n\pi/L\)，而 \(J_v(j_{v,s}r/a)\) 在圓盤上的正交性（以 \(r\) 作為權重）使得係數可以通過積分計算，\(c_s=\frac{2}{a^2 J_{v+1}^2(j_{v,s})}\int_0^a f(r)\,J_v\!\left(j_{v,s}r/a\right) r\,dr\)。
零點如何變化
對於固定的階數，增加索引 \(s\) 會使 \(j_{v,s}\) 以接近 \(\pi\) 的步長增加（遠離原點時振盪變得幾乎週期性）。對於固定的索引，增加階數 \(v\) 會將第一個零點向外推，對於較大的 \(v\)，大致遵循 \(j_{v,1}\approx v + 1.8557\,v^{1/3}+\cdots\)，反映了貝塞爾方程中向心加速度項 \(v^2/x^2\) 如何延遲振盪的開始。實用的結論：較高的角向複雜度（\(v\)）和更多的徑向節點（\(s\)）都對應於較大的特徵值，因此頻率更高或衰減更快。
對於一個具體的鼓例子，其中 \(c=100\ \text{m/s}\) 和 \(a=0.30\ \text{m}\)，基本頻率為 \(f_{0,1}=\frac{100}{2\pi(0.30)}\,(2.4048256)\approx 127.6\ \text{Hz}\)。第一個泛音使用 \(j_{1,1}=3.8317060\)，給出 \(\approx 203\ \text{Hz}\)，因此泛音與基本音的比率為 \(j_{1,1}/j_{0,1}\approx 1.593\)——聽覺上不是八度音程或五度音程，這就是為什麼與弦樂器相比，鼓的聲音不成調。
常見問題
為什麼 \(Y_{v}\) 的第一個零點會小於 1？\(Y_{0}\) 在原點處發散至 \(-\infty\)，並在大約 \(x \approx 0.894\) 處穿越零值，之後才到達第一個極大值，因此它的第一個零點遠小於 \(J_{0}\) 的第一個零點。
\(v\) 可以是非整數嗎？可以。\(Y_{v}\) 的定義公式對任何非整數 \(v\) 都直接成立，而整數階則以連續極限的方式處理。
計算結果的精度如何？計算採用雙精度浮點運算，對於中等範圍的 \(v\) 與 \(s\)，約可提供 10 位有效數字。
