moneta <- c("TESTA","CROCE") prove <- 100000 num <- 0; den <- 0 for (i in 1:prove) { out <- sample(moneta,size=10,replace=TRUE) den <- den + 1 croci <- sum(out=="CROCE") teste <- sum(out=="TESTA") if (croci>teste) num <- num+1 } cat(num,den,num/den,"\n")