home *** CD-ROM | disk | FTP | other *** search
/ Oakland CPM Archive / oakcpm.iso / cpm / basic / math.ark / LSTBIN.BAS < prev    next >
Encoding:
BASIC Source File  |  1986-12-21  |  896 b   |  41 lines
  1. 10 LPRINT
  2. 20 LPRINT "REM: THIS IS 'LSTBIN', SINGLE TERMS OF THE BINOMIAL DISTRIBUTION"
  3. 30 DEFDBL L
  4. 40 PRINT
  5. 50 LPRINT
  6. 60 INPUT; "N:"; N
  7. 70 PRINT
  8. 80 LPRINT "N=";N
  9. 90 INPUT; "R:";R
  10. 100 PRINT
  11. 110 LPRINT "R=";R
  12. 120 PRINT "RUNTIME =" .0908*R "SECONDS"
  13. 130 PRINT
  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 Y = Y + LOG(X)/LOG(10)
  20. 200 NEXT
  21. 210 M = N - R + 1
  22. 220 FOR Z = M TO N
  23. 230 P = P + LOG(Z)/LOG(10)
  24. 240 NEXT Z
  25. 250 NCR = P - Y
  26. 260 S = 1 - T
  27. 270 Q = R*LOG(T)/LOG(10)
  28. 280 U = (N-R)*LOG(S)/LOG(10)
  29. 290 LOGPROB = NCR + Q + U
  30. 300 LPRINT
  31. 310 LPRINT
  32. 320 LPRINT "LOGPROB =" LOGPROB
  33. 330 LPRINT "PROBABILITY =" 10^(LOGPROB-INT(LOGPROB)), "TEN"INT(LOGPROB)
  34. 340 LPRINT
  35. 350 END
  36. 
  37. 320 LPRINT "LOGPROB =" LOGPROB
  38. 330 LPRINT "PROBABILITY =" 10^(LOGPROB-INT(LOGPROB)), "TEN"INT(LOGPROB)
  39. 340 LPRINT
  40. 350 END
  41.