\begin{Schunk}
\begin{Sinput}
  hist(m2000, prob = TRUE, xlab = "x", main = "Means with N=2000")
  lines(density(m2000), lty = 2)
  mp <- mean(m2000)
  sdp <- sd(m2000)
  xp <- seq(range(m2000)[1], range(m2000)[2], length = 200)
  npd <- dnorm(xp, m = mp, s = sdp)
  lines(xp, npd, lty = 1)
  legend("topright",legend=c("Observed","Theoretical"), lty=c(2,1))
\end{Sinput}
\end{Schunk}
