pausa <- function() { cat ("\n >> press Enter to continue\n"); scan() }

# n=10000; x <- sample(outcomes, n, replace=TRUE); hist(x, nc=40, col='cyan')

outcomes <- as.vector(outer(1:6,1:6))
N = 100
n = 10000
sx = rep(0, n)
for(ns in 1:N) {
  xi <- sample(outcomes, n, replace=TRUE)
  sx = sx + xi
  hist(sx, nc=100, col='cyan', xlab='Sum X', freq=FALSE, 
       main=sprintf("ns = %d; mean = %.1f, std = %.1f ",
                    ns, mean(sx), sd(sx)))
  if(ns>3) plot(function(x) dnorm(x, mean(sx), sd(sx)), add=TRUE, lwd=3,
                xlim=c(mean(sx)-4*sd(sx), mean(sx)+4*sd(sx)), col='magenta')
  pausa()
}