import scipy.special as sp def bj(m, x): "mth-order Bessel function" return sp.jv(m, x)