binom <- function(n) { plot(0:n, dbinom(0:n, n, 1/2), xlim=c((n/2)-2*sqrt(n),(n/2)+2*sqrt(n))) Sys.sleep(0.05) } ignore <- sapply(1:200, binom)