\begin{Schunk}
\begin{Sinput}
 hist(ex133$x1, xlab = "x1", prob = TRUE, main = "Sample 133", ylim = c(0, 0.07))
 mtext("Histogram with KDE smooth", 3, -1)
 lines(density(ex133$x1))
 legend("right", legend = paste(c("mean =", "std.dev ="), round(c(mean(ex133$x1), sd(ex133$x1, 2)),2)))
\end{Sinput}
\end{Schunk}
