home *** CD-ROM | disk | FTP | other *** search
/ Oakland CPM Archive / oakcpm.iso / sigm / vol283 / math.lbr / BINOM-4.BAS < prev    next >
BASIC Source File  |  1986-12-18  |  1KB  |  43 lines
  1.    10 REM: NEVADA BASIC 12 DIGIT PRECISSION PROGRAM
  2.    20 LPRINT "PROGRAM BINPROBS: SINGLE TERMS OF THE BINOMIAL DISTRIBUTION"
  3.    30 REM N IS NUMBER OF TRIES
  4.    40 INPUT "N=",N
  5.    50 PRINT 
  6.    60 LPRINT "N =";N
  7.    70 REM R IS NUMBER OF SUCCESSES
  8.    80 INPUT "R=",R
  9.    90 PRINT 
  10.   100 LPRINT "R =";R
  11.   110 PRINT "RUNTIME =";.23*R;"seconds"
  12.   120 PRINT 
  13.   130 REM THETA IS PROBABILITY OF SUCCESS PER TRY
  14.   140 INPUT "THETA =",T
  15.   150 LPRINT "THETA =";T
  16.   160 PRINT 
  17.   170 PRINT "CALCULATIONS IN PROGRESS"
  18.   180 FOR X=1 TO R
  19.   190 LET Y=Y+LOG10(X)
  20.   200 NEXT X
  21.   210 LET M=N-R+1
  22.   220 FOR Z=M TO N
  23.   230 LET P=P+LOG10(Z)
  24.   240 NEXT Z
  25.   250 LET C=P-Y
  26.   260 LET S=1-T
  27.   270 LET Q=R*LOG10(T)
  28.   280 LET U=(N-R)*LOG10(S)
  29.   290 LET L=C+Q+U
  30.   300 PRINT "LOGPROB =";L,"H = ";-L/N/LOG10(4)
  31.   310 LPRINT "LOGPROB =";L,"H = ";-L/N/LOG10(4)
  32.   320 PRINT "probability =";10^(L-INT(L)+1);"ten";INT(L)-1
  33.   330 LPRINT "PROBABILITY =";10^(L-INT(L)+1);"TEN";INT(L)-1
  34.   340 LET G2=10-INT(LOG10(ABS(L)))-INT(LOG10(1.25*R))
  35.   350 PRINT "CALCULATIONS GOOD TO ";G2;" DIGITS"
  36.   360 LPRINT "CALCULATIONS GOOD TO ";G2;" DIGITS"
  37.   370 PRINT "--oo000oo--"
  38.   380 LPRINT "--oo000oo--"
  39.   390 PRINT 
  40.   400 LPRINT 
  41.   410 END 
  42. ATIONS GOOD TO ";G2;" DIGITS"
  43.   370 PRINT "--oo0