\begin{Schunk}
\begin{Sinput}
 myx1 <- rnorm(1000, 40, 10)
 myx2 <- rnorm(1000, 40, 20)
 mx1 <- round(mean(myx1, na.rm=T), 2)
 sdx1 <- round(sd(myx1, na.rm=T), 2)
 mx2 <- round(mean(myx2, na.rm=T), 2)
 sdx2 <- round(sd(myx2, na.rm=T), 2)
 par(mfcol=c(1,2))
 hist(myx1 , breaks=30, prob=T, xlab="Small Variance", ylab="Density", main="", xlim=c(-20,100),ylim=c(0,0.035))
 text(5, 0.026, label=paste("Mean = ", round(mx1,2))) 
 text(5, 0.024, label=paste("Std. Dev. = ", round(sdx1,2)))
 hist(myx2 , breaks=30, prob=T, xlab="Big Variance", ylab="Density", main="", xlim=c(-20,100), ylim=c(0,0.035))
 text(70, 0.026, label=paste("Mean = ", round(mx2,2))) 
 text(70, 0.024, label=paste("Std. Dev. = ", sdx2))
\end{Sinput}
\end{Schunk}
