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

    clc, clear
syms x1(t) y1(t) k11 k12 k21 N a x2(t) y2(t) T 
[x1,y1]=dsolve(diff(x1)==-(k11+k12)*x1+N*a/T,diff(y1)==k12*x1-k21*y1,...
    x1(0)==0,y1(0)==0);
x1=simplify(x1), pretty(x1)
y1=simplify(y1), pretty(y1)
[x2,y2]=dsolve(diff(x2)==-(k11+k12)*x2,diff(y2)==k12*x2-k21*y2,...
    x2(0)==subs(x1,t,T)+N*a,y2(0)==subs(y1,t,T));
x2=simplify(x2), pretty(x2),
y2=simplify(y2), pretty(y2)