n <- 100 b <- 5 p <- matrix(0,b,n) plot(c(0,100), c(0.5,0.5),type="l",lty=2,ylim=c(0,1),xlab="Flip Number",ylab="Proportion Heads") for (i in 1:b) lines(seq(1,n,1),cumsum(rbinom(n=100,size=1,prob=0.5))/seq(1,n,1))