MINICERN/mathlib/gen/v/rnmnml.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:56  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       SUBROUTINE RNMNML(N,NSUM,PCUM,NVEC,IERR)
  11
  12       PARAMETER (NRV = 100)
  13       DIMENSION PCUM(*),NVEC(*),RV(NRV)
  14
  15       IERR=2
  16       IF(PCUM(N) .NE. 1) GO TO 9
  17       IERR=1
  18       DO 1 I = 2,N
  19       IF(PCUM(I) .LT. PCUM(I-1)) GO TO 9
  20     1 CONTINUE
  21       IERR=0
  22       DO 2 K = 1,N
  23     2 NVEC(K)=0
  24       J0=NSUM/NRV+1
  25       I0=NRV
  26       DO 3 J = 1,J0
  27       IF(J .EQ. J0) I0=MOD(NSUM,NRV)
  28       CALL RANLUX(RV,I0)
  29       DO 4 I = 1,I0
  30       DO 5 K = 1,N
  31       IF(PCUM(K) .GE. RV(I)) THEN
  32        NVEC(K)=NVEC(K)+1
  33        GO TO 4
  34       ENDIF
  35     5 CONTINUE
  36     4 CONTINUE
  37     3 CONTINUE
  38     9 RETURN
  39       END