scatola <- c("ROSSA","ROSSA","ROSSA","ROSSA","BLU","BLU","BLU") prove <- 100000 num1 <- 0; num2 <- 0; den <- 0 for (i in 1:prove) { out <- sample(scatola,size=2,replace=FALSE) den <- den + 1 if (out[1]=="ROSSA") num1 <- num1+1 if (out[2]=="ROSSA") num2 <- num2+1 } cat(num1,num2,den,num1/den,num2/den,"\n")