This paper presents a new numerical method for the inversion of the Laplace transform, where only real values of the transform are known. This method, based on the approximation of the transform by truncated series of orthogonal functions, related to the Jacobi or Laguerre polynomials, avoids the loss of significant figures inherent in earlier methods.