Computing the probability density function of the stable Paretian distribution

The practical implementation of a stable Paretian model is a nontrivial task, because—with the exception of a few special cases—its probability density function cannot be expressed analytically. Here, we present an algorithm for calculating the probability density function of the asymmetric stable Paretian distribution. Due to the use of the Fast Fourier Transform, the algorithm is computationally efficient, easily implemented, and of similar accuracy as existing algorithms.