R from Part 27

urn = c(0,0,1,1,1)
sum(sample(urn,3,replace=FALSE)) #simulates one Xk

n = 2000
sum(replicate(n,sum(sample(urn,3,replace=FALSE))))/n
#simulates average (X1+...+Xn)/n

m = 2000
outcomes = replicate(m, sum(replicate(n,sum(sample(urn,3,replace=FALSE))))/n)
hist(outcomes,breaks=50)