fun074.m

来自「matlab讲义」· M 代码 · 共 10 行

function f=fun(x,t1)
t=[0.083 0.167 0.25 0.5 0.75 1.0 1.5 2.25 3 4 6 8 10 12]; 
c=[10.9 21.1 27.3 36.4 35.5 38.4 34.8 24.2 23.6 15.7 8.2 8.3 2.2 1.8];
if nargin==1
   f=c-x(1)*x(3)*(exp(-x(2)*t)-exp(-x(1)*t))/(x(1)-x(2));
else
   plot(t,c,'r+')
   f=x(1)*x(3)*(exp(-x(2)*t1)-exp(-x(1)*t1))/(x(1)-x(2));
end