Given a fixed number $n$ of observations on a variate $x$ which has the Inverse Gaussian probability density function $$\exp\big\{-\frac{\phi^2x}{2\lambda} + \phi ...