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