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