这个计算器能做什么
本工具可以计算第一类球贝塞尔函数 jv(x)、第二类球贝塞尔函数 yv(x)，以及它们的一阶导数 j'v(x) 和 y'v(x)。这些函数是波动方程和亥姆霍兹方程在球坐标系下的径向解，在物理学中应用十分广泛：散射理论、电磁波与声波辐射，以及量子力学（自由粒子的分波展开）都会用到它们。重建后的版本支持非负整数阶 \(v \ge 0\) 和实数自变量 \(x > 0\)。
使用方法
输入阶数 v（非负整数，如 0、1、2）和自变量 x（正实数），点击计算即可同时得到这四个结果。需要注意的是，当 \(x\) 趋近于 0 时，yv(x) 和 y'v(x) 会发散，因此工具在 \(x = 0\) 处将它们记为无穷大；而 \(j_0(0)\) 作为极限值等于 1。
公式解析
这些函数满足方程 $$x^2 w'' + 2x\,w' + (x^2 - v(v+1))w = 0.$$ 从封闭形式 \(j_0 = \dfrac{\sin x}{x}\) 和 \(y_0 = -\dfrac{\cos x}{x}\) 出发，更高阶的值可通过三项递推公式 $$f_{n+1} = \frac{2n+1}{x}\,f_n - f_{n-1}$$ 求得。对 yv 来说，向上递推是数值稳定的；但对 jv 而言，当 \(n > x\) 时向上递推会不稳定，因此我们采用 Miller 向下递推法：从某个较高阶开始，将 \(f\) 分别设为 0 和 1，逐步向下递推，最后将所有值统一缩放，使 0 阶项与 \(\dfrac{\sin x}{x}\) 相匹配。导数则通过 $$j'_v = j_{v-1} - \frac{v+1}{x}\,j_v$$ 计算。
实例演算（v = 0，x = 2）
\(j_0(2) = \dfrac{\sin 2}{2} = 0.4546487134\)。\(y_0(2) = -\dfrac{\cos 2}{2} = 0.2080734183\)。由于 \(j'_0 = -j_1\)，可得 \(j'_0(2) = -0.4353977750\)；而 \(y'_0 = -y_1\) 则给出 \(0.3506120043\)。
常见问题
支持复数 x 吗？ 不支持。原版页面可以接受复数自变量，而本次重建版本为了简洁和运算速度，仅限于实数 \(x > 0\)。
为什么 yv 在 x = 0 处是无穷大？ 第二类函数在原点处存在极点，因此当 \(x\) 趋近于 0 时，其值会无限增大。
计算精度如何？ 计算采用双精度浮点数，可提供约 15 位有效数字，对于一般的工程和物理计算而言绰绰有余。
