Idea - not using differential equation model, but direct curve fitting. We turn to this laters.


We could solve just the differential equation for a(t) - exponential decay - and then use this with antidifferentiation (e.g., integrating factor) to get b(t) and then easily get c(t).


We set up differential or rate equations for the chemical reactions assuming rate of reaction is proportional to amount of chemical present.

We assume the reactions are first order.

Input := 


eqa = a'[t] == - k1 a[t] 

a'[t] == -(k1 a[t])

Output =


a'[t] == -(k1 a[t])

Output =


a'[t] == -(k1 a[t])

Input := 


eqb = b'[t] == - k2 b[t] + k1 a[t]

b'[t] == k1 a[t] - k2 b[t]

Output =


b'[t] == k1 a[t] - k2 b[t]

Output =


b'[t] == k1 a[t] - k2 b[t]

Input := 


eqc = c'[t] == k2 b[t]

c'[t] == k2 b[t]

Output =


c'[t] == k2 b[t]

Output =


c'[t] == k2 b[t]