Abstract
A numerical method for determining the cumulative probability distribution of a nonnegative random variable is based on the steepest descent approximation of the inverse Laplace transform of its moment-generating function. Good numerical agreement with the cumulative exponential and Poisson distributions is demonstrated.