double cyl_bessel_j(double nu, double x); float cyl_bessel_jf(float nu, float x); long double cyl_bessel_jl(long double nu, long double x);
Jν(x)=∞∑k=0(−1)k(x/2)ν+2kk!Γ(ν+k+1) , for x≥0,
where ν is nu and x is x.