// gsl routine to calculate random numbers (0,1)

     #include <iostream>
     #include <gsl/gsl_rng.h>
     
     int  main ()
     {
       const gsl_rng_type * T;
       gsl_rng * r;
     
       int i, n = 10; // generate 10 random numbers
     
       gsl_rng_env_setup();
     
       T = gsl_rng_default;
       r = gsl_rng_alloc (T);
     
       for (i = 0; i < n; ++i) 
         {
           double u = gsl_rng_uniform (r);
	   std::cout << u << std::endl; 
         }
     
       gsl_rng_free (r);
     
       return 0;
     }