Scribd
Upload
19 views2 pages

MATLAB

The document models radioactive decay over time using a system of differential equations to calculate the decay of nuclei N1, N2, and N3. It plots the decay curves of the three nuclei on a graph with time on the x-axis and number of nuclei on the y-axis. The graph shows an exponential decay in the number of nuclei N1 over time, with numbers transferring to N2 and then N3, peaking at a set time point.

Uploaded by

akram1978
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
19 views2 pages

MATLAB

The document models radioactive decay over time using a system of differential equations to calculate the decay of nuclei N1, N2, and N3. It plots the decay curves of the three nuclei on a graph with time on the x-axis and number of nuclei on the y-axis. The graph shows an exponential decay in the number of nuclei N1 over time, with numbers transferring to N2 and then N3, peaking at a set time point.

Uploaded by

akram1978
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

%Decay Problem

L1=0.1;L2=0.2;t=0:0.1:50;
[N1, N2, N3] = dsolve('DN1=-L1*N1','DN2=L1*N1-
L2*N2','DN3=L2*N2','N1(0)=100','N2(0)=0','N3(0)=0')
N1=eval(N1);
N2=eval(N2);
N3=eval(N3);
plot(t,N1,t,N2,t,N3)
hold
plot(t(N2==max(N2)),max(N2),'x')
title('Radio N1 to N2 & N3')
xlabel('Time')
ylabel('Number of Nuclei')
legend('N1','N2','N3')

>> DR

N1 =

100*exp(-L1*t)

N2 =

(exp(-L1*t)*exp(-L2*t)*((100*L1*L2*exp(L1*t))/(L1 - L2) - (100*L1*L2*exp(L2*t))/(L1 - L2)))/L2

N3 =

exp(-L1*t)*exp(-L2*t)*(100*exp(L1*t)*exp(L2*t) - (100*L1*exp(L1*t))/(L1 - L2) +


(100*L2*exp(L2*t))/(L1 - L2))

Current plot held

>>
Radio N1 to N2 & N3
100
N1
N2
90 N3

80

70

60
Number of Nuclei

50

40

30

20

10

0
0 5 10 15 20 25 30 35 40 45 50
Time

You might also like