[R] Zipf random number generation

[Weiguang Shi]

Thank you very much again. Your code is really helpful
and it demontrates the power of R to a new-comer like
me.

 --- Ted.Harding at nessie.mcc.ac.uk wrote: 
> Understood! Though I think you mean
> 
>   1/c = sum(1.0 / pow((double) j, alpha)),
> j=1,2,...,N

You are right.

Weiguang