\begin{Schunk}
\begin{Sinput}
> onexset <- rnorm(500, m=truemean, sd=truesd)
> m1 <- round( mean(onexset), 2)
> sd1 <- round(sd(onexset),2)
> hist(onexset, breaks=20, prob=T, xlim=c(truemean-4*truesd, truemean+4*truesd), ylim=c(0,ymax), main="",xlab="Observations from one sample")
> legend("topleft", c(paste("Mean =", m1), paste("StdDev =", sd1)))
\end{Sinput}
\end{Schunk}
