\begin{Schunk}
\begin{Sinput}
 var1 <- rexp(n = 1500, rate = 1/50)
 hist(x = var1, prob = T, breaks = 20, xlim = c(-10, 
     300), ylim = c(0, 0.03), xlab = "An Exponential Random Sample", 
     ylab = "Proportion of Observations", main = "")
 den1 <- density(var1)
 lines(den1, lty = 2, col = "red")
 legend("topleft", legend = c(paste("mean=", 
     round(mean(var1), 3)), paste("sd=", round(sd(var1), 
     3))))
\end{Sinput}
\end{Schunk}
