  This papers discusses the leaky aquifer function considered in a recent paper
by Frank Harris in the Journal of Computational and Applied Mathematics (2008).
We describe properties of an integral representing this function and give
details on how to compute this function with a single algorithm for a wide
range of the parameters.
