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