9000  REM  *****  BINOPO  *****  STATISTICS PROGRAM  *****
9001  REM  *****  VERSION 1  *****  7/31/69  *****
9002  REM  COMPARISON OF BINOMIAL, NORMAL, AND POISSON DISTRIBUTIONS
9003  PRINT "PLEASE INPUT THE NUMBER OF TRIALS - N";
9004  INPUT N
9005  PRINT "PLEASE INPUT THE PROBABILITY - P";
9006  INPUT P
9007  PRINT "N="N,"P="P
9008  LET Y=0
9009  PRINT 
9010  PRINT "J","EXACT","NORMAL","POISSON"
9011  PRINT 
9012  LET M=N*P
9013  LET S=SQR(M*(1-P))
9014  LET R=(1-P)^N
9015  LET V=EXP(-M)
9016  FOR J=0 TO N
9017  IF J=0 THEN 9019
9018  LET R=R*P*(N-J+1)/J/(1-P)
9019  LET T=INT(R*10000+.5)/10000
9020  LET H=(J-.5-M)/S
9021  LET K=(J+.5-M)/S
9022  LET C=1/SQR(2)
9023  LET A1=.141128
9024  LET A2=.0886403
9025  LET A3=.0274335
9026  LET A4=-.00039446
9027  LET A5=.00328975
9028  DEF FNO()=1-1/(1+A1*+A2*^2+A3*^3+A4*^4+A5*^5)^8
9029  IF H<0 THEN 9032
9030  LET F=.5+.5*FNO(H*C)
9031  GOTO 9033
9032  LET F=.5-.5*FNO(-H*C)
9033  IF K<0 THEN 9036
9034  LET G=.5+.5*FNO(K*C)
9035  GOTO 9037
9036  LET G=.5-.5*FNO(-K*C)
9037  LET L=.0001*INT(10000*(G-F)+.5)
9038  IF J=0 THEN 9040
9039  LET V=V*M/J
9040  LET Q=.0001*INT(10000*V+.5)
9041  IF T>0 THEN 9047
9042  IF L>0 THEN 9047
9043  IF Q>0 THEN 9047
9044  IF Y=0 THEN 9052
9045  PRINT J,T,L,Q
9046  GOTO 9999
9047  IF J=0 THEN 9050
9048  IF Y=1 THEN 9051
9049  PRINT J-1,0,0,0
9050  LET Y=1
9051  PRINT J,T,L,Q
9052  NEXT J
9053  STOP 
9999  END 
