double sph_bessel(unsigned n, double x); float sph_besself(unsigned n, float x); long double sph_bessell(unsigned n, long double x);

jn(x)=(Ï€/2x)1/2Jn+1/2(x) ,\quad {}for xâ‰¥0,

where n is n and x is x.