# il problema del Cavalier de Méré. dado <- c(1:6) prove <- 100000 num.1 <- 0; num.2 <- 0; den <- 0 for (i in 1:prove) { l <- sample(dado,size=4,replace=TRUE) ris <- sum(l==1) if (ris>0) num.1 <- num.1 + 1 tmp.1 <- sample(dado,size=24,replace=TRUE) tmp.2 <- sample(dado,size=24,replace=TRUE) tmp.3 <- tmp.1+tmp.2 ris <- sum(tmp.3==2) if (ris>0) num.2 <- num.2 + 1 den <- den + 1 } cat(num.1,num.2,den,num.1/den,num.2/den,"\n")