/*-------------------------------------------------------------

This program demonstrate the use GSL for generating  
Gassian Random Numbers in range [0-1] 

           - Jayanti Prasad, Sept 13, 2011    
             Comments & Feedback : prasad.jayanti@gmail.com 
--------------------------------------------------------------*/
#include <stdio.h>
#include <gsl/gsl_rng.h>
#include<gsl/gsl_randist.h>

int  main (int argc, char *argv[]){
  const gsl_rng_type * T;
  gsl_rng * r;
  int i, n;
  double sigma,g;

  if (argc < 3){
     fprintf(stderr,"./gran <#  points> <sigma> \n");
      return(-1);
  }
  
  n = atoi(argv[1]);
  sigma=atoi(argv[2]); 

  gsl_rng_env_setup();
  
  T = gsl_rng_default;
  r = gsl_rng_alloc (T);
  
  for (i = 0; i < n; i++) {
     g = gsl_ran_gaussian (r,sigma); 
     fprintf(stdout,"%.6f\n", g);
    }
  
  gsl_rng_free (r);
  
  return 0;
}