www.gusucode.com > 司守奎《Lingo软件及应用习题解答》的源程序 > 司守奎《Lingo软件及应用习题解答》的源程序/《Lingo软件及应用习题解答》程序和数据/12第12章/Lt12_10_9.m

    clc, clear, close all
gamma=57.22; alpha=2.008; beta=0.1855;
t0=(log(alpha)-log(beta))/(alpha-beta)
ft=@(t,N)gamma*N*(exp(-beta*t)-exp(-alpha*t));
hold on
for i=1:5
    fplot(@(t)ft(t,i),[0,16],'k')
end