The authors discuss in detail the usefulness of continued fraction expansions for the numerical computation of the incomplete gamma function
(a,z) for complex values of a and z. Using functional relations for
(a,z), the complementary function
(a,z) and the (related) special Kummer function M(l,b,z), they introduce what they call "additional (or extended)" continued fraction expansions, which they show to be stable when evaluated by the backward recurrence algorithm. A thorough error analysis is presented and the special cases of
(a,z) which correspond to the error functions erf(z), erfc(z) and the exponential integral
are discussed in detail. In connections with the complex error function, the "anomalous convergence" of general T-fractions (the two point Padé approximants) is discussed.