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),for x≥0

where n is n and x is x.