I have a set of differential equations that I would like to fit to experimental data to get chemical kinetics parameters out of. The equations are:
d(Mzb(t))/dt = -[Mzb(t)-Mzb0]/T1b - kb*Mzb(t) + kf*Mzf(t);
d(Mzf(t))/dt = -[Mzf(t)-Mzf0]/T1f - kf*Mzf(t) + kb*Mzb(t);
The data I have is of Mzb(t) and Mzf(t) at different time points. From that data, I can glean Mzb0 and Mzf0 (the equilibrium values Mzb(t) and Mzf(t), which will be constants), and from other data I know the constants T1b and T1f. I also know that this will be a least-squares optimization to get the best fit for the kb and kf parameters. The part I don't understand is how to code all of this-- I'm a novice with the syntax and would appreciate any and all help about it.