Project 1

From phys660
Jump to navigationJump to search

Part I for both PHYS460 and PHYS660 students

Consider a radiactive decay problem involving two types of nuclei, A and B, with populations and . Suppose that type A nuclei decay to form type B nuclei, which then also decay, according to differential equations:

,

,

where and are the decay time constants for each type of nucleus. Use the Euler method to solve these coupled equations numerically for and as a function of time.

Note that this problem can also be solved analytically either by using ``paper-and-pencil method or Mathematica. Obtain the analytic solutions for and and compare them with your numerical results. Present your results as graphs, with different plot for each of the three cases

  • (a) \tau_A > \tau_B
  • (b) \tau_A = \tau_B
  • (c) \tau_A < \tau_B

Each plot should contain numerical solutions (using different values for the time step in the Euler algorithm) contrasted with the exact analytic solution. To avoid having to assign too many numerical values, use , =0, and as the unit of time. In particular, try to interpret the short and long time behaviors for different value of this ratio.

Part II for PHYS660 students only

Part II for PHYS 660 only