Writing a user defined function for inverse-Gaussian distribution

Thank you. Now it works well!!