Algorithm 654: Fortran subroutines for computing the incomplete gamma function ratios and their inverses - DiDonato, A. R.; Morris, A. H. (ACM Trans. Math. Software, 1987)
Single-precision Fortran, maximum accuracy 14D.