Computer simulation of geometric stable distributions

We present a new method for computer simulation of strictly geometric stable random variables. The method is based on their representation as a product of two independent random variables with explicit distribution functions, coupled with the inversion method. We also extend the method to the multivariate case, by deriving new representation and simulation algorithm for multivariate Linnik distribution.