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.