这个计算器能做什么
本工具可计算第一类修正贝塞尔函数 \(I_v(x)\) 与第二类修正贝塞尔函数 \(K_v(x)\)，并同时给出它们的一阶导数 \(I'_v(x)\) 和 \(K'_v(x)\)。这两个函数是修正贝塞尔方程 $$x^2 y'' + x y' - (x^2 + v^2)y = 0$$ 的两个线性无关解，在物理与工程领域应用广泛：例如圆柱体中的热传导、扩散问题、传输线与波导理论，以及统计学等。由于这属于纯数学计算，结果在世界各地完全一致，不受任何地区规则影响。
使用方法
输入阶数 v（可为任意实数）和自变量 x。当 v 为整数时，\(I_v(x)\) 对所有实数 x 均可计算；否则要求 \(x \geq 0\)。\(K_v(x)\) 要求 \(x > 0\)，因为当 \(x \to 0^+\) 时该函数发散；当 \(x \leq 0\) 时，结果会标记为未定义（NaN）。
公式详解
\(I_v(x)\) 通过其幂级数求和得到，其中 Gamma 函数采用 Lanczos 近似计算。\(K_v(x)\) 则通过对积分 $$K_v(x) = \int_0^{\infty} e^{-x\cdot\cosh t}\cdot\cosh(vt)\, dt$$ 进行数值积分求得，该方法对整数阶和非整数阶都同样稳定。导数则使用对称递推公式 $$I'_v(x) = \tfrac{1}{2}\left(I_{v-1}(x) + I_{v+1}(x)\right)$$ 和 $$K'_v(x) = -\tfrac{1}{2}\left(K_{v-1}(x) + K_{v+1}(x)\right)$$，从而避免了除以 x 的运算。
计算示例（v = 0，x = 1）
由级数可得 $$I_0(1) = 1 + 0.25 + 0.015625 + \dots \approx 1.26606588$$ 由积分可得 \(K_0(1) \approx 0.42102444\)。由于 \(I_{-1} = I_1\)，对称公式给出 \(I'_0(1) = I_1(1) \approx 0.56515910\)，以及 \(K'_0(1) = -K_1(1) \approx -0.60190723\)。
常见问题
为什么 K_v(x) 显示为未定义？ \(K_v(x)\) 仅在 \(x > 0\) 时有定义；当 x 等于或小于零时，该函数会发散。
可以使用分数阶吗？ 可以。这两个函数都支持任意实数阶，包括非整数和负数。
计算精度如何？ 结果采用双精度浮点运算（约 12–15 位有效数字），在中等大小的 x 取值范围内与标准参考表完全吻合。
