\begin{Schunk}
\begin{Sinput}
> hist(onexset, breaks=20, prob=T, xlim=c(truemean-3*truesd, truemean+3*truesd), ylim=c(0,ymax),   main="",xlab="Observations from one sample")
> lines(x, denx, lty=1)
> lines(onedensity, lty=2)
> legend("topleft", legend=c(paste("Exp. Value=",truemean), paste("Obs. Mean=",round(m1,2))))
> legend("topright", legend=c(paste("Pop. Var.=",round(truesd*truesd,2)), paste("Obs. Var.=",round(sd1*sd1,2))))
\end{Sinput}
\end{Schunk}
