home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 9 / 09.iso / e / e004 / 1.img / DECAY.MCD < prev    next >
Encoding:
Text File  |  1989-11-02  |  1.6 KB  |  81 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. .CMD FILENAME data1 data1
  11. .CMD FILENAME data2 data2
  12. .TXT 1 0 1 20 
  13. a1,19,17,18
  14. RADIOACTIVE DECAY
  15. .TXT 2 9 4 55 
  16. a4,55,53,162
  17. This document shows the results of radioactive decay. 
  18. Element A decays with a half-life of thalfA into 
  19. element B. Element B decays with a half-life of 
  20. thalfB.
  21. .EQN 5 1 1 9 
  22. gm~1M
  23. .EQN 0 20 1 9 
  24. yr~1T
  25. .TXT 2 -30 1 22 
  26. a1,21,19,20
  27. Initial conditions:
  28. .TXT 2 9 1 30 
  29. a1,30,28,29
  30. Initial amount of element A:
  31. .EQN 0 41 1 14 
  32. A0:100*gm
  33. .TXT 2 -41 1 13 
  34. a1,13,11,12
  35. Half-lives:
  36. .EQN 0 41 1 17 
  37. thalfA:10*yr
  38. .EQN 1 0 1 18 
  39. thalfB:100*yr
  40. .TXT 2 -41 2 31 
  41. a2,31,29,44
  42. Ratio: mass of 1 atom of A to 
  43. 1 atom of B:
  44. .EQN 1 41 1 10 
  45. α:1.1
  46. .TXT 2 -50 1 27 
  47. a1,26,24,25
  48. Compute decay constants:
  49. .EQN 2 10 3 14 
  50. kA:ln(2)/thalfA
  51. .EQN 0 40 3 14 
  52. kB:ln(2)/thalfB
  53. .TXT 4 -50 1 19 
  54. a1,18,16,17
  55. Decay functions:
  56. .EQN 2 10 2 21 
  57. A(t):A0*e^-(kA*t)
  58. .EQN 3 0 2 39 
  59. B(t):-α*A0*e^-(kB*t)*(e^-(kA*t)-1)
  60. .TXT 4 -10 1 39 
  61. a1,38,36,37
  62. Graph the mass of A and B over time:
  63. .EQN 2 10 1 25 
  64. t:0*yr,2*yr;200*yr
  65. .EQN 0 35 1 27 
  66. t2:0*yr,10*yr;200*yr
  67. .EQN 2 -45 16 42 
  68. A0&0*gm&A(t),B(t){1,1,15,30,l}@200*yr&0*yr&t
  69. .EQN 0 46 24 7 
  70. t2/yr=
  71. .EQN 0 10 24 8 
  72. A(t2)/gm{19137}=
  73. .EQN 0 10 24 8 
  74. B(t2)/gm{18993}=
  75. .TXT 10 -24 1 3 
  76. a1,3,1,2
  77. B
  78. .TXT 4 0 1 3 
  79. a1,3,1,2
  80. A
  81.