/* compile and link with
gcc -o testcase testcase.c -lgsl -lgslcblas -lm
*/

#include <stdio.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_randist.h>

int main () {
  const gsl_rng_type * T = gsl_rng_default;
  gsl_rng * r = gsl_rng_alloc (T);
  gsl_rng_set(r , 0); 
  
  int IT;
  
  for (IT = 0; IT < 1; IT++)
    {
      double max = GSL_MAX (0.0, gsl_ran_gaussian (r, 2));
      printf ("%g\n", max);
    }
}