home *** CD-ROM | disk | FTP | other *** search
/ Oakland CPM Archive / oakcpm.iso / cpm / basic / math.ark / POIS.BAS < prev    next >
Encoding:
BASIC Source File  |  1986-12-21  |  1.4 KB  |  38 lines
  1.    10 REM THIS IS A NEVADA BASIC 12 DIGIT PROGRAM WRITTEM BY MIKE FINERTY
  2.    20 REM IT WILL OPPERATE OUTSIDE OF THE NORMAL RANGE OF THE LANGUAGE AS
  3.    30 REM SET UP BY ELLIS COMPUTING (TM), USING LOG FACORING TECHNIQUES
  4.    40 LPRINT "PROGRAM POISSON:  SINGLE TERMS OF THE POISSON DISTRIBUTION"
  5.    50 PRINT "PROGRAM POISSON:  SINGLE TERMS OF THE POISSON DISTRIBUTION"
  6.    60 LPRINT 
  7.    70 PRINT 
  8.    80 REM X IS NUMBER OF OCCURENCES
  9.    90 INPUT "X = ",X
  10.   100 LPRINT "X = ";X
  11.   110 PRINT "RUNTIME = ";0.12*X;" SECONDS"
  12.   120 PRINT 
  13.   130 REM LAMBDA IS AVERAGE EXPECTATION
  14.   140 INPUT "LAMBDA = ",L
  15.   150 LPRINT "LAMBDA = ";L
  16.   160 PRINT "CALCULATIONS IN PROGRESS"
  17.   170 FOR N=1 TO X
  18.   180 LET Y=Y-LOG10(N)
  19.   190 NEXT N
  20.   200 LET Z=X*LOG10(L)
  21.   210 LET R=-L/LOG(10)
  22.   220 REM P IS THE LOGARITHM OF THE POISSON PROBABILITY
  23.   230 LET P=Y+Z+R
  24.   240 LPRINT 
  25.   250 LPRINT "LOGPROB = ";P
  26.   260 LPRINT "POISSON PROBABILITY = ";10^(P-INT(P)+1);"TEN";INT(P)-1
  27.   270 LET G2=10-INT(LOG10(ABS(P)))-INT(LOG10(1.25*X))
  28.   280 LPRINT "CALCULATIONS GOOD TO ";G2;" DIGITS"
  29.   290 PRINT "LOGPROB = ";P
  30.   300 PRINT "POISSON PROBABILITY = ";10^(P-INT(P)+1);"TEN";INT(P)-1
  31.   310 PRINT "CALCULATIONS GOOD TO ";G2;" DIGITS"
  32.   320 PRINT "--ooOOOoo--"
  33.   330 LPRINT "--oo000oo--"
  34.   340 END 
  35. INT(P)-1
  36.   310 PRINT "CALCULATIONS GOOD TO ";G2;" DIGITS"
  37.   320 PRINT "--ooOOOoo--"
  38.   330 LPRINT "--oo000oo--"