[Help-gsl] Zipf Random Number Generator

[Michael Duerr]

Hi,

I'm looking for a function that can produce random numbers that are 
Zipf-distributed?
Something like:

int zipfian(double skew, int N),

where skew is the skew factor (e.g. 0.4) and N the upper limit of the range 
(e.g. [0, N) )

Does gsl hold something similar or can you provide me a reference for a sample 
implementation?

Thanks,
Michael