rndWeibull ============================================== Purpose ---------------- Computes Weibull pseudo-random numbers with the choice of underlying random number generator. Format ---------------- .. function:: x = rndWeibull(r, c, shape, scale) { x, newstate } = rndWeibull(r, c, shape, scale, state) :param r: number of rows of resulting matrix. :type r: scalar :param c: number of columns of resulting matrix. :type c: scalar :param shape: shape parameter, ExE conformable matrix with *r* and *c*. :type shape: matrix or vector or scalar :param scale: scale parameter, ExE conformable matrix with *r* and *c*. :type scale: matrix or vector or scalar :param state: Optional argument. **scalar case** *state* = starting seed value only. If -1, GAUSS computes the starting seed based on the system clock. **opaque vector case** :*state* = the state vector returned from a previous call to one of the rnd random number functions. :type state: scalar or opaque vector :return x: Weibull distributed random numbers. :rtype x: RxC matrix :return newstate: the updated state. :rtype newstate: Opaque vector Remarks ------- The properties of the pseudo-random numbers in *x* are: .. math:: E(x) = scale * \Gamma(1 + 1/shape)\\ Var(x) = scale^2*( \Gamma(1 + \frac{2}{shape}) - (\Gamma(1 + \frac{1}{shape}))^2 ) *r* and *c* will be truncated to integers if necessary. .. seealso:: Functions :func:rndCreateState, :func:rndStateSkip