Computing the incomplete gamma function