\begin{Schunk}
\begin{Sinput}
 par(mfcol=c(3,1))
 for (i in 1:3){
      minx <- 20 + 50 * (i-1)
      xvals <- seq(minx,300, length.out=1000)
      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(150, 0.7*max(gam), c("gamma(sh=2,sc=1)", "inv gauss(mu=2,lambda=5)"), lty = 1:2,)
  }
\end{Sinput}
\end{Schunk}
