dado <- c(1:6)
nlanci <- 100000
out <- sample(dado,size=3*nlanci,replace=TRUE)
out <- out==3
ris <- matrix(out,ncol=3)
ncr <- apply(ris,1,sum)
(frel <- table(ncr)/sum(table(ncr)))
mean(ncr)
var(ncr)