home *** CD-ROM | disk | FTP | other *** search
/ Piper's Pit BBS/FTP: ibm 0020 - 0029 / ibm0020-0029 / ibm0028.tar / ibm0028 / MCAD254.ZIP / HISTO.MCD < prev    next >
Encoding:
Text File  |  1989-06-21  |  1.5 KB  |  69 lines
  1. .MCD 25000 0
  2. .CMD SURFACEFORMAT rot=10 tilt=35 vScale=20 size=15,30
  3. .CMD SKETCHFORMAT mag=1.000000,1.000000 center=0.500000,0.500000 size=15,30 box=y
  4. .CMD PLOTFORMAT logs=0,0 subdivs=1,1 size=5,15 type=l
  5. .CMD FORMAT  rd=d ct=10 im=i et=3 zt=15 pr=3 mass length time charge
  6. .CMD SET ORIGIN 0
  7. .CMD SET TOL 0.001000
  8. .CMD MARGIN 0
  9. .CMD LINELENGTH 78
  10. .CMD SET PRNCOLWIDTH 8
  11. .CMD SET PRNPRECISION 4
  12. .TXT 1 0 1 33 
  13. a1,33,31,32
  14. RANDOM COIN FLIPS ... Histogram
  15. .EQN 2 1 1 15 
  16. trials:200
  17. .TXT 0 19 2 28 
  18. a2,27,25,44
  19. 200 trials, each consists 
  20. of 6 coin flips.
  21. .EQN 2 -19 1 12 
  22. flips:6
  23. .EQN 2 0 1 17 
  24. j:1;trials
  25. .EQN 2 0 1 16 
  26. k:1;flips
  27. .EQN 2 0 4 24 
  28. Heads[j:k$flip(1)
  29. .TXT 1 26 1 31 
  30. a1,31,29,30
  31. ... Number of heads per trial
  32. .TXT 2 5 1 48 
  33. a1,48,46,47
  34. (See below for definition of "flip" function.)
  35. .TXT 2 -32 1 29 
  36. a1,29,27,28
  37. The data looks like this...
  38. .EQN 2 2 10 34 
  39. 7&-1&Heads[j{1,1,9,25,+}@trials&0&j
  40. .TXT 11 -1 1 25 
  41. a1,25,23,24
  42. Now create a histogram.
  43. .EQN 2 1 1 21 
  44. k1:0;flips+1
  45. .EQN 2 0 2 19 
  46. intervals[k1:k1
  47. .TXT 0 22 1 22 
  48. a1,21,19,20
  49. ... Interval limits
  50. .EQN 3 -22 1 31 
  51. freq:hist(intervals,Heads)
  52. .EQN 0 33 2 23 
  53. freq[(last(freq)+1):0
  54. .EQN 3 -33 16 59 
  55. 80&0&freq[k1{1,1,15,50,b}@flips+2&-1&k1
  56. .TXT 0 61 2 11 
  57. a2,11,9,21
  58. Frequency 
  59. table ...
  60. .EQN 3 7 10 8 
  61. freq[k1=
  62. .EQN 1 -5 9 5 
  63. k1=
  64. .TXT 13 -65 1 74 
  65. a1,73,71,72
  66. To see another 200 trials, put cursor in equation below and press [F9].
  67. .EQN 2 20 1 29 
  68. flip(coin)~Φ(rnd(1)-.5)
  69.