/* mc_pi.c */
double mc_pi(long m, gsl_rng *r);