#!/usr/bin/Rscript
n <- 10000
par(mfrow=c(3,2))
for(a in c(1, 0.75, 0.5))
{
values <- sapply(1:100, function(x) 1/n^a*sum(rnorm(n)))
plot(values, ylab=paste("a =", a))
hist(values)
}