I want to show that, when $\mu=\mu_0$, then gamma family $\Gamma(a,b)$ is a conjugate prior to inverse Gaussian with density $f(x,\mu,\lambda)=\sqrt{\frac{\lambda}{2\pi x^2}}exp[-\frac{\lambda(x-\mu)^2}{2\mu^2x}]$.
I was able to verify that $IG$ is an exponential family generated by $\bf{T}$$(X)=-1/2(X,X^{-1})^T$ and $h(x)=\sqrt{\frac{1}{2\pi x^3}}$ but I am pretty confused on how to show the conjugate prior.
Any helps would be deeply appreciated!!