Although the series expansions in §§8.7, 8.19(iv), and 8.21(vi) converge for all finite values of , they are cumbersome to use when is large owing to slowness of convergence and cancellation. For large the corresponding asymptotic expansions (generally divergent) are used instead. See also Luke (1975, pp. 101–102) and Temme (1994b).
DiDonato and Morris (1986) describes an algorithm for computing and for , , and from the uniform expansions in §8.12. The algorithm supplies 14S accuracy. A numerical inversion procedure is also given for calculating the value of (with 10S accuracy), when and are specified, based on Newton’s rule (§3.8(ii)). See also Temme (1987, 1994b).
Expansions involving incomplete gamma functions often require the generation of sequences , , or for fixed and . An efficient procedure, based partly on the recurrence relations (8.8.5) and (8.8.6), is described in Gautschi (1979b, 1999).