# include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_complex.h>
#include <gsl/gsl_complex_math.h>
#include <gsl/gsl_randist.h>

main (void)
{
  double  x[650]; /*Random numbers with a gaussian distribution*/
  int N; /*N is the number of samples = 600*/
  int i,k,count;/*Counters*/
  int run; /*run the code 1000 times*/
  N = 600;
  for(run = 1;run<1001;run++)
     {	
       /*Generate x[] with a Gaussian distribution*/
       gsl_rng_env_setup();
       T = gsl_rng_default;
       r = gsl_rng_alloc (T);
       seed = gsl_rng_default_seed;
       seed = seed;/*Update seed in every run run*/
       for(k=0;k<N;k++)
	 {
	   x[k]=gsl_ran_gaussian(r,1);
	 }
       gsl_rng_free(r);
     }
}