double sph_neumann(unsigned n, double x); float sph_neumannf(unsigned n, float x); long double sph_neumannl(unsigned n, long double x);

nn(x)=(Ï€/2x)1/2Nn+1/2(x),for xâ‰¥0

where n is n and x is x.