LCI = 1000 - 3*4^(1/2)

LCS = 1000 + 3*4/4^(1/2)

k = replicate(10000, (mean(rnorm(4, mean=1000, sd=4))))

k[k < LCI | k > LCS]

n = which(k %in% c(k[k < LCI | k > LCS]))

n[9]