invgrnd
              Random arrays from the inverse Gaussian distribution.
 r = invgrnd (mu, lambda) returns an array of random
 numbers chosen from the inverse Gaussian distribution with location parameter
 mu and scale parameter lambda.  The size of r is the common
 size of mu and lambda.  A scalar input functions as a constant
 matrix of the same size as the other inputs.
 When called with a single size argument, invgrnd returns a square
 matrix with the dimension specified.  When called with more than one scalar
 argument, the first two arguments are taken as the number of rows and columns
 and any further arguments specify additional matrix dimensions.  The size may
 also be specified with a row vector of dimensions, sz.
 The inverse Gaussian CDF is only defined for mu > 0 and
 lambda > 0.
Further information about the inverse Gaussian distribution can be found at https://en.wikipedia.org/wiki/Inverse_Gaussian_distribution
See also: invgcdf, invginv, invgpdf, invgfit, invglike
Source Code: invgrnd