## Number of simulation replications.
nrep = 1e4

## Sample size.
n = 10

## Normally-distributed data.
Z = rnorm(nrep*n)
Z = array(Z, c(nrep,n))

Z_mean = apply(Z, 1, mean)
Z_med = apply(Z, 1, median)

## The first row of the results.
V = c(var(Z_mean), var(Z_med))

for (df in c(2, 4, 6, 8, 10)) {

    ## T2 distributed data.
    T = rt(nrep*n, df)
    T = array(T, c(nrep,n))

    T_mean = apply(T, 1, mean)
    T_med = apply(T, 1, median)

    V = rbind(V, c(var(T_mean), var(T_med)))
}