sim_gig function

Generate Generalized Inverse Gaussian Distribution