A method for numerical inversion of the Laplace transform based on Tikhonov regularization is described. The method requires values of u(s) only for real s, and appears relatively reliable in practice, since the ill posedness of this problem has been considered. The precondition to use the algorithm is that the exact solution of the inversion exists and it is Lipschitz continuous with finite (and known) support.