R from Part 27

urn = c(0,1)
n = 20
N = 10000

hist(replicate(N, sum(sample(urn, size = n, replace=TRUE))), breaks = seq(-0.5,n+0.5,1), xlim=c(0,n))