10  COM X[103,22],M[19,19],U[19,19],Q[20],V[20],P[78]
30  COM M$[60],N$[72]
40  COM N,K,N8,K8,N9,K9,Q9,Q7,Q5,Q4,Q3,Q2,Q1
45  COM I3,I4,U9,X$[20]
50  REM:13OCT73
52  REM:>GAUS
100  DIM A$[12]
220  PRINT 
222  GOTO Q9 OF 230,230,250
230  PRINT "* WHAT IS THE MEAN OF THE NORMAL VARIABLE X ? ";
240  GOTO 260
250  PRINT "* MEAN   = ?";
260  ENTER 30,Q8,G9
280  IF Q8>0 THEN 320
285  PRINT 
290  PRINT "TO CALCULATE PROBABILITIES FOR A NORMAL DISTRIBUTION, YOU"
300  PRINT "MUST SPECIFY THE MEAN AND STANDARD DEVIATION OF THAT DIS-"
310  PRINT "TRIBUTION. YOU WILL BE ASKED FOR THESE SEPARATELY."
315  GOTO 250
320  PRINT 
322  GOTO Q9 OF 330,330,350
330  PRINT "        WHAT IS THE STANDARD DEVIATION OF X  ";
340  GOTO 360
350  PRINT "STD.DEV. = ";
360  INPUT G8
370  IF G8>0 THEN 400
380  PRINT "STANDARD DEVIATION MUST BE POSITIVE";
390  GOTO 320
400  PRINT 
402  GOTO Q9 OF 410,410,430
410  PRINT "* WHAT IS THE LOWER LIMIT L FOR THE CALCULATION ?";
420  GOTO 440
430  PRINT "* LOWER LIMIT L = ?";
440  ENTER 30,Q8,G4
450  IF Q8>0 THEN 500
455  PRINT 
460  PRINT "THE CALCULATION TO BE MADE IS"
470  PRINT "PR(L<X<U) FOR A NORMALLY DISTRIBUTED VARIABLE X."
480  PRINT "YOU WILL BE ASKED SEPARATELY FOR L AND U."
485  PRINT "FOR EXAMPLE, TO OBTAIN A LEFT-TAIL CUMULATIVE PROBABILITY,"
486  PRINT "SET L AT 10 STANDARD DEVIATIONS BELOW THE MEAN."
487  PRINT "TO OBTAIN THE DENSITY AT X, SET L = U = X."
490  GOTO 400
500  PRINT 
502  GOTO Q9 OF 510,510,530
510  PRINT TAB(22);"WHAT IS THE UPPER LIMIT U ";
520  GOTO 540
530  PRINT "  UPPER LIMIT U = ";
540  INPUT G7
545  T=(G7-G9)/G8
550  IF G4=G7 THEN 690
560  IF G4<G7 THEN 595
570  PRINT LIN(1)"UPPER LIMIT CANNOT BE SMALLER THAN LOWER LIMIT!"
580  PRINT "TRY AGAIN."
585  PRINT 
590  GOTO 400
595  Z=0
600  S=(G4-G9)/G8
601  IF ABS(T)>4 AND ABS(T)<5 THEN 604
602  IF ABS(S)>4 AND ABS(S)<5 THEN 604
603  GOTO 609
604  Z=1
609  GOSUB 1000
610  G7=G4
620  G4=G2
630  GOSUB 1000
640  G2=G4-G2
645  PRINT 
650  PRINT  USING 660;G2
660  IMAGE "PR(L<X<U) = ",2D.5D
671  IF Z#1 THEN 730
672  PRINT '10"CAUTION:  IF L OR U IS BETWEEN 4 AND 5 STANDARD DEVIATIONS"
673  PRINT "FROM THE MEAN, PROBABILITIES MAY BE SLIGHTLY IN ERROR IN"
674  PRINT "THE FINAL DECIMAL PLACE."
680  GOTO 730
690  R=EXP(-.5*T^2)
700  R=R/(2.50663*G8)
705  PRINT 
710  PRINT  USING 712;R
712  IMAGE "DENSITY AT X = ",.5DE
730  PRINT 
731  GOTO Q9 OF 732,732,735
732  PRINT "WANT MORE PROBABILITIES ";
734  GOTO 740
735  PRINT "MORE ";
740  INPUT A$
750  IF A$[1,1]="N" THEN 9998
760  GOTO 400
1000  REM: GAUSSIAN PROBABILITIES
1010  G6=(G7-G9)/G8
1020  IF G6>-4.5 THEN 1050
1030  G2=0
1040  GOTO 1230
1050  IF G6<4.5 THEN 1120
1060  G2=1
1070  GOTO 1230
1080  IF G6 >= 0 THEN 1120
1090  G5=-1
1100  G6=-G6
1110  GOTO 1130
1120  G5=1
1130  G2=G6
1140  G3=G6
1150  FOR I=1 TO 100
1160  G3=-(G3*G6^2*(2*I-1))/((2*I+1)*2*I)
1170  G2=G2+G3
1180  IF ABS(G3)<.000005 THEN 1200
1190  NEXT I
1200  G2=.5+G2/SQR(6.28318)
1210  IF G5 >= 0 THEN 1230
1220  G2=1-G2
1230  RETURN 
9998  CHAIN "$IDA",150
9999  END 
