A numerical method for computing inverse Laplace transforms is proposed. In
this method, the complex contour integral defining the inverse transform i
s computed over an equivalent contour as proposed by Talbot and Evens. Spec
ial contours, called optimal contours, are constructed so that the transfor
med real integrand decreases exponentially to zero as z runs along such a c
ontour to infinity. The efficient Clenshaw-Curtis quadrature is employed fo
r the final evaluation. The presented method is competitive and compares fa
vourably with those of Talbot and Evans.