Description: In the present paper, we propose and investigate a numerical method of computing the probability distribution of S. Inversion of the Laplace transform moment generating function is shown to be a better approach than the existing recursive method for find the distribution of sum of several independent random variables. In the continuous case, in many situations, numerical inversion needs to be used to compute the probability distribution function of the sum.