A variation
M = colMeans(replicate(100, rnorm(25, mean = 10, sd = 10))) hist(M, xlab = "sample means")
giving