rgpd function

Generalized Pareto simulation