fn gamma_series(a: f64, x: f64) -> f64
Series expansion for the regularized incomplete gamma function.