\begin{Schunk}
\begin{Sinput}
 library(statmod)
 mu <- 1
 lambda <- 2
 xrange<- seq(from=0.0,to=2*mu+5/lambda,by=0.02)
 mainlabel<- expression(paste("IG(",mu,",",lambda,")",sep=""))
 xprob <- dinvgauss(xrange, mu = mu, lambda = lambda)
 plot(xrange, xprob, type = "l", main = mainlabel, xlab = "possible values of x", ylab = "probability of x")
 
\end{Sinput}
\end{Schunk}
