double pareto (double x, void *p) { double * c = (double *)p; return gsl_ran_pareto_pdf (x, c[0], c[1]); }
double test_pareto_pdf (double x) { return gsl_ran_pareto_pdf (x, 1.9, 2.75); }