evrnd
              Random arrays from the extreme value distribution.
 r = evrnd (mu, sigma) returns an array of random
 numbers chosen from the extreme value distribution (also known as the Gumbel
 or the type I generalized extreme value distribution) with location
 parameter mu and scale parameter sigma.  The size of r is
 the common size of mu and sigma.  A scalar input functions as a
 constant matrix of the same size as the other inputs.
 When called with a single size argument, evrnd 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 Gumbel distribution is used to model the distribution of the maximum (or
 the minimum) of a number of samples of various distributions.  This version
 is suitable for modeling minima.  For modeling maxima, use the alternative
 Gumbel iCDF, gumbelinv.
Further information about the Gumbel distribution can be found at https://en.wikipedia.org/wiki/Gumbel_distribution
See also: evcdf, evinv, evpdf, evfit, evlike, evstat
Source Code: evrnd