You can use a function such as
ranGenPareto =
Compile[{{loc, _Real}, {scale, _Real}, {shape, _Real}},
Module[
{u},
u = RandomReal[];
loc + scale*(u^(-shape) - 1)/shape
],
CompilationTarget -> "C", RuntimeAttributes -> {Listable}
];
to generate random numbers from the distribution. The above assumes you have a C compiler on your computer to speed up things.
