graphics.off();
cat("\014");
rm(list=ls());


blockLength = 10000000;
nbins = 200;
a = rnorm(blockLength);
hist(a,freq=FALSE, breaks=nbins, col="blue",xlab="x", 
     ylab = expression('f'[X]*'(x)'),
     main="Gaussian PDF")