Project 1

Part I for both PHYS460 and PHYS660

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

$\frac{d{N}_A}{dt}=-\frac{N_A}{{\tau }_A}$,

$\frac{d{N}_B}{dt}=\frac{N_A}{{\tau }_A}-\frac{N_B}{{\tau }_B}$,

where ${\tau }_A$ and ${\tau }_B$ are the decay time constants for each type of nucleus. Use the Euler method to solve these coupled equations numerically for $N_A(t)$ and $N_B(t)$ as a function of time.