drunk <- function(n) { t <- 0:n p <- dbinom(t, n, 1/4) x <- 2*t-n plot(x, p, xlim=c(-60,60), ylim=c(0,1)) Sys.sleep(0.2) } ignore <- sapply(0:60, drunk)