Laplace transform inversions using optimal contours in the complex plane

A numerical method for computing inverse Laplace transforms is proposed. In this method, the complex contour integral defining the inverse transform is computed over an equivalent contour as proposed by Talbot and Evans. Special contours, called optimal contours, are constructed so that the transformed real integrand decreases exponentially to zero as z runs along such a contour to infinity. The efficient Clenshaw-Curtis quadrature is employed for the final evaluation. The presented method is competitive and compares favourably with those of Talbot and Evans.