V112 Multinomial Random Numbers

OBSOLETE

Please note that this routine has been obsoleted in CNL 218. Users are advised not to use it any longer and to replace it in older programs. No maintenance for it will take place and it will eventually disappear.
Suggested replacement: RNMNML (V138)

MUNOMI generates a vector of random integers $n$_i (i=1,2,...,N) with probabilities $p$_i

according to a multinomial law:

$Prob(n$₁,n₂,...,n_N) = {(n₁+n₂+...+n_N)!n₁! n₂! ... n_N!} p₁^n₁ p₂^n₂ ... p_N^n_N .

Structure:

SUBROUTINE subprogram
User Entry Names: MUNOMI
External References: RNDM (V104), UZERO (V300)

Usage:

    CALL MUNOMI(NCH,NTOT,P,N,IERR)

NCH

( INTEGER) Number N of random integers $n$_i

requested.

NTOT

( INTEGER) Equals $\sum$_i=1^Nn_i , specified by the user.

P

( REAL) One-dimensional array of length NCH at least. On entry, it contains in P(i) the probability of channel i. On return, it contains the cumulative channel probabilities so that $P(NCH)=1$ . If $P(NCH)=1$ on entry, it is assumed that P(i) contains the cumulative probabilities rather than the individual probabilites, which saves some time.

N

( INTEGER) One-dimensional array of length NCH at least. On return, N(i) contains the generated random integers.

IERR

Error flag.
$=\; 0:$ Normal case,
$=\; 1:$ At least one ,
$=\; 2:$ $\sum P(i)\; >\; 1$ .

V113