Se tienen n cartas numeradas de 1 a n.

Se barajan y se ponen en fila, se tiene un exito si la carta i está en la posición i

X(va) cuenta el número de exitos.

Estimar, usando simulación, la media y la desviación estándar de X

n: número de cartas

nsim: número de simulaciones.

carta <- function(n, nsim) {
    x <- numeric(nsim)
    for (i in 1:nsim) {
        reparte <- sample(1:n)
        pos <- 1:n
        x[i] = sum(reparte == pos)
    }
    return(c(mean(x), sd(x)))
}
carta(100, 10000)

## [1] 1.012 1.014

system.time(carta(100, 10000))

##    user  system elapsed 
##   0.250   0.007   0.256