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.