\begin{Schunk}
\begin{Sinput}
  hist(x, prob= TRUE, xlab="1000 draws from Tweedie(3,50,1)", ylab="Density", main="Tweedie(3,50,1)")
  legnd <- c( paste("mean = ", round(mean(x),2)), paste("var =", round(var(x),2)), paste("sd = ", round(sd(x),2)))
  legend("topright", legend=legnd)
\end{Sinput}
\end{Schunk}
