pub fn pgamma_raw(a: f64, x: f64) -> f64
Regularized lower incomplete gamma function P(a, x) = γ(a,x) / Γ(a). Uses series expansion for x < a+1, continued fraction otherwise.