urna <- c(1:5)
nest <- 100000
out <- matrix(nrow=nest,ncol=3)
for (i in 1:nest) out[i,] <- sample(urna,size=3,replace=FALSE)
ris <- apply(out,1,min)
frequenze(ris)
mean(ris)
var(ris)