#include <stdio.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>
#include <time.h>     

     int
     main (void)
     {
       const gsl_rng_type * T;
       gsl_rng * r;
     
       int i, n = 600;
       int run;
       double mu = 1.0;
       double k;
       double meanv;
       long seed;
     
       /* create a generator chosen by the 
          environment variable GSL_RNG_TYPE */
     
       gsl_rng_env_setup();
     
       T = gsl_rng_default;
       r = gsl_rng_alloc (T);
       seed = gsl_rng_default_seed;

       /* print n random variates chosen from 
          the gaussian distribution with standard deviation 
          parameter mu */

       for(run=1;run<11;run++)
	 {
	   seed = seed + 0.456;
	   meanv = 0;
	   for (i = 0; i < n; i++) 
	     {
	       k = gsl_ran_gaussian (r, mu);
	       meanv=k+meanv;
	     }
	    printf("%lf\n",meanv); /*Print mean value after generating 600 samples*/
	 }
       printf ("\n");
       gsl_rng_free (r);
       return 0;
     }