我的问题与在 Python中使用球形贝塞尔函数有关。

如果我的 ODE 包含以下形式的球面贝塞尔函数

对于给定的值$t$和$x$.

我使用了以下内容：

from scipy.special import spherical_jn(n, z, derivative)
from scipy.special import spherical_yn(n, z, derivative)

我想我可以将 jn 定义如下：

def Jn(tx,n):
    return (sqrt(pi/(2*(t*x))*jn(n+0.5,t*x))

然后将其输入到我的 ODE 中。

任何帮助在这里表示赞赏。如果需要，我可以发布 ODE。