Technique for computing the PDFs and CDFs of non-negative infinitely divisible random variables

We present a method for computing the PDF and CDF of a non-negative infinitely divisible random variable $X$. Our method uses the L'{e}vy-Khintchine representation of the Laplace transform $\mathbb{E}{e}^{-\lambda X}=e^{-\varphi (\lambda )}$, where $\varphi$ is the Laplace exponent. We apply the Post-Widder method for Laplace transform inversion combined with a sequence convergence accelerator to obtain accurate results. We demonstrate this technique on several examples including the stable distribution, mixtures thereof, and integrals with respect to non-negative L'{e}vy processes. Software to implement this method is available from the authors and we illustrate its use at the end of the paper.