#include "stdio.h"
#include "stdlib.h"
#include "math.h"
#include "randomlib.h"

#define N 1000
#define M 100000

int main(int argc,char **argv)
{
   int i;
   double r,rmin=1e32,rmax=-1e32;
   double bins[N+1];
   double sum=0,sum2=0;

   for (i=0;i<N;i++)
      bins[i] = 0;

   RandomInitialise(1802,9373);
   for (i=0;i<M;i++) {
      r = RandomGaussian(1.0,0.5);
      sum += r;
      sum2 += (r*r);
      /* printf("%g\n",r); */
      if (r < -4)
         r = -4;
      if (r > 4)
         r = 4;
      bins[(int)(N/2 + r * (N/2) / 4.0)]++;
   }
   fprintf(stderr,"Mean = %g\n",sum / M);
   fprintf(stderr,"Standard deviation = %g\n",sqrt((sum2 - sum*sum/M)/M));
/*
   for (i=0;i<N;i++)
      printf("%g %g\n",i/(double)N,bins[i]);
*/
}

#include "randomlib.c"