pub fn gamma_q(a: f64, x: f64) -> f64
Regularized upper incomplete gamma Q(a, x) = Γ(a, x)/Γ(a).
Q(a, x) = Γ(a, x)/Γ(a)