B = NULL

## Number of simulation replications.
nrep = 1000

## Sample size.
n = 20

## Trim k values off each end.
for (k in c(0,1,2,3,4,5))
{
  X = array(rnorm(nrep*n), c(nrep,n))
  X = apply(X, 1, sort)
  X = t(X)  ## R peculiarity.
  TX = X[,(k+1):(20-k)]
  TM = apply(TX, 1, mean)
  B = c(B, mean(TM))
}