prob.R <- c(1/6,1/6,1/6,1/6,1/6,1/6) prob.T <- c(1/10,1/10,1/10,1/10,1/10,1/2) dado <- c(rep("R",50),rep("T",50)) prove <- 100000 num <- 0; den <- 0 for (i in 1:prove) { out <- sample(dado,size=1) if (out=="R") lanci <- sample(c(1:6),size=2,replace=TRUE, prob=prob.R) if (out=="T") lanci <- sample(c(1:6),size=2,replace=TRUE, prob=prob.T) if (sum(lanci)==12) { den <- den + 1 if (out=="T") num <- num+1 } } cat(num,den,num/den,"\n")