\begin{Schunk}
\begin{Sinput}
  hist(x, prob = TRUE, breaks = 50,  xlab="Number of Successes Out of 2000 Experiments", main="481 Draws from Binom(2000, 0.2)")
  lines(density(x), lty = 2)
  xseq <- seq( min(x), max(x), by = 1)
  xprob <- dnorm(xseq, m = 2000*0.2, sd = sqrt(2000*(0.2)*(0.8)))
  lines(xseq, xprob, lwd = 6)
  xbinprob <- dbinom(xseq, size = 2000, prob = 0.2)
  lines(xseq, xbinprob, col = 2, lwd = 2, lty = 1)
  legend("topright", legend = c("observed KDE","Binom Exact", "Normal Approx"), lty = c(2,1,1), lwd = c(1,2,3), col = c(1,2,1))
\end{Sinput}
\end{Schunk}
