\begin{Schunk}
\begin{Sinput}
 xvals <- seq(0,20, length.out=200)
 gam <- dgamma(xvals, shape=2, scale=1)
 igaus <- dinvgauss(xvals,mu=2,lambda=5)
 plot(xvals,gam,type="l",lty=1,main="")
 lines(xvals,igaus,lty=2)
 legend(6,.2,c("gamma(sh=2,sc=1)","inv gauss(mu=2,lambda=5)"),lty=1:2,)
 
\end{Sinput}
\end{Schunk}
