  We extend the technique of using the Trapezoidal Rule for efficient
evaluation of the Special Functions of Mathematical Physics given by integral
representations. This technique was recently used for Bessel functions, and
here we treat Incomplete Gamma functions and the general Confluent
Hypergeometric Function.
