my_gig function

Draw from generalized inverse Gaussian