The Laplace transform is often applied to linear partial differential equations to eliminate the time dimension. The analytical solutions thus obtained need to be inverted to the time domain. This paper provides Mathematica functions for approximate inversion of the Laplace transform.