\begin{Schunk}
\begin{Sinput}
  xrange <- c(m-4*sd, m+4*sd)
  xseq <- seq(xrange[1], xrange[2], length.out = 100)
  dxseq <- dnorm(xseq, m = m, sd = sd)
  onedensity <- density(onexset)
  ymax <- 1.2 * max(dxseq)
  hist(onexset, breaks = 20, prob = TRUE, xlim = c(m-3*sd, m+3*sd), ylim = c(0,ymax),   main="",xlab = "x")
  lines(xseq, dxseq, lty = 1)
  lines(onedensity, lty=2)
  legend("topleft", legend = c(paste("Exp. Value=", m), paste("Obs. Mean=", round(m1,2))))
  legend("topright", legend = c(paste("Pop. Var.=", round(sd*sd,2)), paste("Obs. Var.=",round(sd1*sd1,2))))
  legend("left", legend = c("Theoretical","Observed"), lty = c(1,2))
\end{Sinput}
\end{Schunk}
