rndKMnb ============================================== Purpose ---------------- Computes negative binomial pseudo-random numbers. Format ---------------- .. function:: { x, newstate } = rndKMnb(r, c, k, p, state) :param r: number of rows of resulting matrix. :type r: scalar :param c: number of columns of resulting matrix. :type c: scalar :param k: "event" argument for negative binomial distribution, scalar or ExE conformable matrix with *r* and *c*. :type k: matrix or vector or scalar :param p: "probability" argument for negative binomial distribution, scalar or ExE conformable matrix with *r* and *c*. :type p: matrix or vector or scalar :param state: **scalar case** *state* = starting seed value only. If -1, GAUSS computes the starting seed based on the system clock. **500x1 vector case** *state* = the state vector returned from a previous call to one of the ``rndKM`` random number functions. :type state: scalar or 500x1 vector :return x: negative binomial distributed random numbers. :rtype x: RxC matrix :return newstate: the updated state. :rtype newstate: 500x1 vector Remarks ------- The properties of the pseudo-random numbers in *x* are: .. math:: E(x) = \frac{k * p}{1 - p} Var(x) = \frac{k * p}{(1 - p)^2}\\ x = 0, 1,....k > 0\\ 0 < p < 1\\ *r* and *c* will be truncated to integers if necessary. Technical Notes --------------- :func:`rndKMnb` uses the recur-with-carry KISS+Monster algorithm described in the :func:`rndKMi` Technical Notes. Source ------ randkm.src