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