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