程序：Mathematica 求解模型Ⅰ程序：
Matlab 拟合曲线程序： Curvefun.m文件 function f=curvefun(x,t)
f=103.333*x(1)*(exp(-x(2)*t)-exp(-x(1)*t))/(x(1)-x(2));
Run.m文件
t=[0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5 6 7 8 9 10 11 12 13 14 15 16]; c=[30 68 75 82 82 77 68 68 58 51 50 41 38 35 28 25 18 15 12 10 7 7 4]; x0=[2.4,0.16];
x=lsqcurvefit('curvefun',x0,t,c)
在命令窗口中运行run即可求得拟合结果。
Mathematica 作图程序：
Mathematica 求解模型Ⅱ程序：
Mathematica 求解问题一程序：
16
Mathematica 求解问题二程序：
Mathematica 求解问题三程序：
Mathematica 求解问题四程序：
17