pub fn I0(x: f64) -> f64
This routine computes the regular modified cylindrical Bessel function of zeroth order, I_0(x)