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.