8999 REM ******* REVISED APRIL, 1972  BY ED BAKER *******
9000 REM ***** BINOPO ***** STAISTICS PROGRAM *****
9001 REM ***** VERSION1 ***** 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:PRINT'N='N,'P='P
9008 Y=0
9009 PRINT
9010 PRINT "J","EXACT","NORMAL","POISSON"
9011 PRINT
9012 M=N*P
9013 S=SQR(M*(1-P))
9014 R=(1-P)^N
9015 V=EXP(-M)
9016 FOR J=0 TO N
9017 IF J=0 THEN 9019
9018 R=R*P*(N-J+1)/J/(1-P)
9019 T=INT(R*10000+.5)/10000
9020 H=(J-.5-M)/S
9021 K=(J+.5-M)/S
9022 C=1/SQR(2)
9023 A1=.141128
9024 A2=8.86403E-02
9025 A3=2.74335E-02
9026 A4=-3.9446E-04
9027 A5=3.28975E-03
9028 DEF FNO(S)=1-1/(1+A1*X+A2*X^2+A3*X^3+A4*X^4+A5*X^5)^8
9029 IF H<0 THEN 9032
9030 F=.5+.5*FNO(-H*C)
9031 GOTO 9033
9032 F=.5-.5*FNO(-H*C)
9033 IF K<0 THEN 9036
9034 G=.5+.5*FNO(K*C)
9035 GOTO 9037
9036 G=.5-.5*FNO(-K*C)
9037 L=.0001*INT(10000*(G-F)+.5)
9038 IF J=0 THEN 9040
9039 V=V*M/J
9040 Q=.0001*INT(10000*V+.5)
9041 IF T>0 THEN 9047
9042 IF L>0 THEN 9047
9043 IF Q>0 THEN 9046
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 Y=1
9051 PRINT J,T,L,Q
9052 NEXT J
9999 END