moneta <-c("TESTA", "CROCE")
nlanci <- 100000
out <- sample(moneta,size=3*nlanci,replace=TRUE)
out <- out=="CROCE"
ris <- matrix(out,ncol=3)
ncr <- apply(ris,1,sum)
(frel <- table(ncr)/sum(table(ncr)))