Rand() doesn't work for NormalInverseGaussian

Please quote your code, and include an MWE.