# il paradosso di Monty Hall.
nrep <- 100000
num <- 0; den <- 0
nsca <- 3 # si puņ mettere quante scatole si vuole
x <- c(1:nsca)

for (i in 1:nrep)
{
    scatola <- rep(0,nsca)
    euro <- sample(x,1); scatola[euro] <- 1
    punta <- sample(x,1)
    if (scatola[punta]==1) chiusa <- sample(x[-euro],1)
    if (scatola[punta]!=1) chiusa <- euro
    if (scatola[chiusa]==1) num <- num + 1
    den <- den + 1
}

cat(num,den,num/den,"\n")