pub fn cyl_bessel_j(nu: f64, x: f64) -> f64
Cylindrical Bessel function of the first kind, J_nu(x)
Uses the series expansion: J_nu(x) = sum_{m=0}^∞ (-1)^m / (m! * Gamma(nu+m+1)) * (x/2)^(2m+nu)