My question relates to using spherical Bessel functions in Python.
If my ODE contains a spherical Bessel function of the form $$j_\ell(tx)$$ and similarly $$y_\ell(tx)$$
for given values of $t$ and $x$.
I used the following:
from scipy.special import spherical_jn(n, z, derivative)
from scipy.special import spherical_yn(n, z, derivative)
I thought I could maybe define jn as follows:
def Jn(tx,n):
return (sqrt(pi/(2*(t*x))*jn(n+0.5,t*x))
And then just input it into my ODE.
Any help here is appreciated. I can post the ODE if required.