Exponential decay problem.

14 views (last 30 days)
Adnan Sardar
Adnan Sardar on 24 Oct 2018
Commented: David Goodmanson on 30 Oct 2018
Radioactive decay is modeled with the exponential function f(t)=f(0)e^(kt), where t is time, f(0) is the amount of material at t=0, f(t) is the amount of material at time t, k is a constant. If 100 mg are present at t=0, determine the amount that is left after 7 days. Material is Gallium-67, which has a half-life of 3.261 days. Write a script file for the problem. The program should first determine the constant k, then calculate f(7).
  8 Comments
Adnan Sardar
Adnan Sardar on 30 Oct 2018
% k = constant
% D7 = amount of Gallium-67 left after 7-days (grams)
format compact
syms k
k=solve(50==100*exp(k*3.261),k);
D7=round(double(100*exp(k*7)),1)
David Goodmanson
David Goodmanson on 30 Oct 2018
For sure. I would round D7 to two decimal places giving you four significant figures, since the half life that was provided has four sig figs. Lots of people are using symbolic calculation now but in this case taking the log of both sides gives
1/2 = exp(k*t_half)
log(1/2) = k*t_half
k = log(1/2)/t_half

Sign in to comment.

Answers (0)

Categories

Find more on Particle & Nuclear Physics in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!