[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / v / binomi.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1996/04/01 15:02:55  mclareni
* Mathlib gen
*
*
#include "gen/pilot.h"
      SUBROUTINE BINOMI(NTOT,P,N,IERROR)
C    BINOMIAL GENERATOR
C    CODED FROM LOS ALAMOS REPORT      LA-5061-MS
C    PROB(N)=FACT(NTOT)*P**N*(1.-P)**(NTOT-N)/(FACT(N)*FACT(NTOT-N))
C        WHERE FACT(N) STANDS FOR FACTORIAL OF N
C    ON RETURN IERROR.EQ.0 NORMALLY
C              IERROR.EQ.1 IF P.LE.0. OR P.GE.1.
C
      IF(P.GT.0. .AND. P.LT.1.) GO TO 100
      IERROR=1
      GO TO 999
  100 CONTINUE
      IERROR=0
      N=0
      DO 200 I=1,NTOT
      IF(RNDM(N).GT.P) GO TO 200
      N=N+1
  200 CONTINUE
  999 CONTINUE
      RETURN
      END
Commit	Line	Data
fe4da5cc	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)PN(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