Incomplete-gamma-func
evaluates the incomplete gamma functions in the normalized form
with x ≥ 0 and a > 0, to a user-specified accuracy. With this normalization, P(a, x)+Q(a, x) = 1. The function returns with machine precision as relative accuracy.
For more information please review the s14bac function in the NAG document.