pub fn qgamma_raw(a: f64, x: f64) -> f64
Regularized upper incomplete gamma function Q(a, x) = 1 - P(a, x).