double logistic (double x, void *p) { double * c = (double *)p; return gsl_ran_logistic_pdf (x, c[0]); }
double test_logistic_pdf (double x) { return gsl_ran_logistic_pdf (x, 3.1); }