Returns the upper incomplete gamma function Gamma(a,x) = int(exp(-t)t^(a-1),t=0..x) for real a > 0, x > 0.
double gammaUpperIncomplete(double a, double x) => gammaUpperRegularized(a, x) * gamma(a);