摘要

We present a method for computing the probability density PDF) and the cumulative distribution CDF) of a nonnegative infinitely divisible random variable X. Our method uses the Levy-Khintchine representation of the Laplace transform Ee(-lambda X) = e(-phi(lambda)) where phi 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 nonnegative Levy processes.

  • 出版日期2011-3