# Compare binomial with normal n <- 16 x <- 0:n plot(x, xlim=c(-1, n+1), ylim=c(0,0.2), type="n", xlab="x", ylab="dens(X=x)", main="Comparison of S16 with N(8,4)") points(x, dbinom(x, n, 0.5), pch=16) x <- (-1000:(1000*(n+1)))/1000 points(x, dnorm(x, mean=n/2, sd=sqrt(n/4)), type="l")