population.size = 1000000
population.mean = 300
population.std = 50
plot.delay = 0.0
maxsample.size = 100

sample.var = 0

population.data = rnorm(population.size,population.mean,population.std)
population.mean = mean(population.data)

par(mar=c(1,1,1,1))

# plot more samples
for (plot.counter in 2:maxsample.size)
{
  sample.data = sample(population.data,plot.counter,replace = TRUE)
  sample.var[plot.counter-1] = var(sample.data)
}

barplot(sample.var)