caselle <- 5
prove <- 100000
fass <- rep(NA,prove)
for (i in 1:prove)
{
    out <- rep(0,caselle)
    nfig <- 0
    repeat {
       ris <- sample(caselle,size=1,replace=FALSE)
       out[ris] <- 1
       nfig <- nfig + 1
       if (sum(out)==caselle) break
    }
    fass[i] <- nfig
}
table(fass)
table(fass)/prove
plot(table(fass)/prove)