home *** CD-ROM | disk | FTP | other *** search
/ Oakland CPM Archive / oakcpm.iso / cpm / basic / math.ark / BINOMIAL.BAS < prev    next >
Encoding:
BASIC Source File  |  1986-12-21  |  1.1 KB  |  55 lines
  1. 10 VAR A,B,LOGPROB,LOG_NCR,N,R,T,X,M,Z,P,Y,Q,U,S =REAL.DOUBLE
  2. 15 VAR C$ = STRING
  3. 20 INPUT "N = ";N
  4. 30 LPRINTER
  5. 40 PRINT "N = ";N
  6. 45 CONSOLE
  7. 50 INPUT "R = ";R
  8. 60 LPRINTER
  9. 70 PRINT "R = ";R
  10. 80 CONSOLE
  11. 90 PRINT "RUNTIME = "; 0.4*R;"SECONDS"
  12. 100 INPUT "THETA =";T
  13. 110 LPRINTER
  14. 120 PRINT "THETA =";T
  15. 130 CONSOLE
  16. 140 PRINT "CALCULATIONS IN PROGRESS"
  17. 150 FOR X = 1 TO R
  18. 160 Y = Y + LOG(X)/LOG(10)
  19. 170 NEXT X
  20. 180 M = N - R + 1
  21. 190 FOR Z = M TO N
  22. 200 P = P + LOG(Z)/LOG(10)   
  23. 210 NEXT Z
  24. 220 LOG_NCR = P - Y
  25. 230 S = 1 - T
  26. 240 Q = R*LOG(T)/LOG(10)
  27. 250 U = (N - R)*LOG(S)/LOG(10)
  28. 260 LOGPROB =  LOG_NCR + Q + U
  29. 270 PRINT
  30. 280 LPRINTER
  31. 290 PRINT
  32. 295 PRINT "LOGPROB =";LOGPROB
  33. 300 CONSOLE
  34. 310 PRINT
  35. 315 PRINT "LOGPROB=";LOGPROB
  36. 320 A = 10^(LOGPROB - INT(LOGPROB))
  37. 330 B = INT(LOGPROB)
  38. 340 C$= "PROBABILITY = #.## TEN +#####"
  39. 350 PRINT USING C$,A,B
  40. 360 LPRINTER
  41. 370 PRINT USING C$,A,B
  42. 380 PRINT "--ooO0Ooo--"
  43. 390 PRINT
  44. 400 PRINT
  45. 410 CONSOLE
  46. 420 PRINT "--ooO0Ooo--"
  47. 430 PRINT
  48. 440 PRINT
  49. 450 END
  50.  
  51.  A,B,LOGPROB,LOG_NCR,N,R,T,X,M,Z,P,Y,Q,U,S =REAL.DOUBLE
  52. 15 VAR C$ = STRING
  53. 20 INPUT "N = ";N
  54. 30 LPRINTER
  55. 40 PRINT "N =