moneta <- c("TESTA","CROCE") dado <- c(1:6) nrep <- 1000 out <- rep(NA,nrep) for (i in 1:nrep) { x <- sample(moneta,1) if (x == "TESTA") { y <- sample(dado,size=1) } else { y <- sample(dado,size=2,replace=TRUE) y <- sum(y) } out[i] <- y } table(out)