double cyl_neumann(double nu, double x); float cyl_neumannf(float nu, float x); long double cyl_neumannl(long double nu, long double x);
Nν(x)=⎧⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎨⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎩Jν(x)cosνπ−J−ν(x)sinνπ,for x≥0 and non-integral νlimμ→νJμ(x)cosμπ−J−μ(x)sinμπ,for x≥0 and integral ν
where ν is nu and x is x.