The practical implementation of a stable Paretian model is a nontrivial tas
k, because-with the exception of a few special cases-its probability densit
y function cannot be expressed analytically. Here, we present an algorithm
for calculating the probability density function of the asymmetric stable P
aretian distribution. Due to the use of the Fast Fourier Transform, the alg
orithm is computationally efficient, easily implemented, and of similar acc
uracy as existing algorithms. (C) 1999 Elsevier Science Ltd. All rights res
erved.