double cyl_bessel_i(double nu, double x); float cyl_bessel_if(float nu, float x); long double cyl_bessel_il(long double nu, long double x);
Iν(x)=i−νJν(ix)=∞∑k=0(x/2)ν+2kk!Γ(ν+k+1) , for x≥0,
where ν is nu and x is x.