This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / v / binomi.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1996/04/01 15:02:55  mclareni
6 * Mathlib gen
7 *
8 *
9 #include "gen/pilot.h"
10       SUBROUTINE BINOMI(NTOT,P,N,IERROR)
11 C    BINOMIAL GENERATOR
12 C    CODED FROM LOS ALAMOS REPORT      LA-5061-MS
13 C    PROB(N)=FACT(NTOT)*P**N*(1.-P)**(NTOT-N)/(FACT(N)*FACT(NTOT-N))
14 C        WHERE FACT(N) STANDS FOR FACTORIAL OF N
15 C    ON RETURN IERROR.EQ.0 NORMALLY
16 C              IERROR.EQ.1 IF P.LE.0. OR P.GE.1.
17 C
18       IF(P.GT.0. .AND. P.LT.1.) GO TO 100
19       IERROR=1
20       GO TO 999
21   100 CONTINUE
22       IERROR=0
23       N=0
24       DO 200 I=1,NTOT
25       IF(RNDM(N).GT.P) GO TO 200
26       N=N+1
27   200 CONTINUE
28   999 CONTINUE
29       RETURN
30       END