double rayleigh (double x, void *p) { return gsl_ran_rayleigh_pdf (x, *(double*)p); }
double test_rayleigh_pdf (double x) { return gsl_ran_rayleigh_pdf (x, 1.9); }