gumbelrnd
              Random arrays from the Gumbel distribution.
 r = gumbelrnd (mu, beta) returns an array of random
 numbers chosen from the Gumbel distribution (also known as the extreme value
 or the type I generalized extreme value distribution) with location
 parameter mu and scale parameter beta.  The size of r is
 the common size of mu and beta.  A scalar input functions as a
 constant matrix of the same size as the other inputs.
 When called with a single size argument, gumbelrnd 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 maxima.  For modeling minima, use the alternative
 extreme value iCDF, evinv.
Further information about the Gumbel distribution can be found at https://en.wikipedia.org/wiki/Gumbel_distribution
See also: gumbelcdf, gumbelinv, gumbelpdf, gumbelfit, gumbellike, gumbelstat, evrnd
Source Code: gumbelrnd