generate gaussians using GSL

randn.h

double randn(double mean, double sigma)
{ 
  gsl_rng *gnr = gsl_rng_alloc(gsl_rng_default);
  return mean + gsl_ran_gaussian(gnr, sigma);
  gsl_rng_free(gnr);
}