moneta <- c("T","C") nlanci <- 200000 campione <- sample(moneta,nlanci,replace=TRUE,prob=c(1/2,1/2)) testa <- campione=="T" nteste <- cumsum(testa) x <- c(1:nlanci) y <- nteste - x/2 k<-c(seq(10,90,10),seq(100,900,100),seq(1000,9000,1000),seq(10000,90000,10000)) plot(x[k],y[k],type="l",log="x") abline(h=0,lty=2) plot(x[k],y[k]/x[k],type="l",log="x") abline(h=0,lty=2)