Using Stable Distributions in Mathematica
Chambers, Mallows, and Stuck (1976) published an algorithm for the generation of stable random variables. It is implemented in the Stable MathLink Package as the function, StableRandom; similar code is implemented in the supplemental package as the function, SRandom. The example below generates a stable random sequence of 10 numbers in the S(α, β, γ, δ; 1) parameterization where p1 is a parameter set, {α, β, γ, δ}.
p1 = {1,3, 0.5, 1, 0.5}
{-1.30558, 8.32448, 2.26377, -0.14893, 0.669163, -0,25539, -0.216304, -0.146468, 0.0379229, 0.133912}
