ПРИЛОЖЕНИЕ 1. РЕАЛИЗАЦИЯ ПОИСКА АНАЛИТИЧЕСКОГО РЕШЕНИЯ СИСТЕМЫ (3.1) В СРЕДЕ MAPLE.

restart;

>Digits:=50 ; >with(linalg);

>dll:= Vm*diff(X(1),1) =-(kl*X(1)+k3*Tg(1)); >

dl2: = Vg*diff(Y<1),1) = (kl+X(1)+k3*Tg(1))*>dl4: = Vm*diff(Tm(l) ,1) *Cm=av*Vv/Gm*

dl3 :=Vg*diff

sys : = (dll,dX2 # dl3,dl4) ; >

sol:=dsolve({sys},(X(l),Y(1),Tg(l),Tm(l)>); >lhs41:=lhs(sol[4]); >lhs31;-lhs(sol[3]);

>lhs21:=lhs(sol[2]); >lhsll:=lhs(sol[1]); >lhs4r:=rhs(sol[4]); >lhs3r:=rhs(sol[3]); >lhs2r:=rhs(sol[2]); >lhslr:=rhs(sol[1]);

>X_l: = 4if% )); >

if4 (lhsll=Y(l) ,lhslr/ifs (lhs21=Y(l) ,lhs2r, 4 if" (lhs31=Y(l) ,lhs3r,lhs4r))>; >Tm_l: = 4f * (lhsll=Tm(l) , lhslr, *if* (lhs21=Tm(l) , lhs2r, 4 if * (lha31=Tm(l) ,lhs3r,lhs4r>)) ; >Tg_l:-*if% (lhsll-Tg(1) , lhslr, *if4 (lhs21=Tg(l) ,lhs2r,*if * (lhs31=Tg(l) ,lhs3r,lhs4r> >>; >

Tg_l:=collect(Tg_l,_C1); >Tg_l:-collect(Tg_l,_C2); >Tg_l:=collect(Tg_l,_C3) ; >Tg__l:=collect(Tg_l,_C4) ; >

Tg_l_Cl: =coeff (Tg_l ,_C1) ; >

Tg_l_C2: =coeff

Tg_l_C3: =co©ff (TgJL,_C3) ; >

Tg_l_C4 : =co©f f (Tg,_C4) ; >

Tm_l:«collect (TmJL»_C1) ; >TmJL: «collect (Tm_l,jC2) ; >

Tm_l: —collect (Tm__l, _C3) ; >TmJL: «collect (Tm_l,_C4) ; >

Tm_l_Cl: =coeff (Tm_l ,_C1) ; >

Tm_l_C2 : =coeff (TmJL ,_C2) ; >

Tm_l_C3: =coeff (Tm_l ,_C3) ; >

Tm_l_C4 : =coef f (Tm_l ,_C4) ; >X_l:=collectX_l:=collect(X_1,_C2> ; >

X_1: =collect (X_l ,_C3 >; >

X_l: «collect (X_1,_C4) ; >

X_1_C1: =coeff (X_l ,_C1) ; >

X_1_C2 : =coeff (X_l ,_C2) ; >X_l_C3:=coeff (X_1,_C3) ; >

X_1_C4 : «coeff (X_l ,_C4) ; >Y_1: «collect (Y_1,_C1) ; >

Y_1: =collect (Y_1 ,_C2) ; >Y_jL:«collect (Y_1,_C3); >Y_l:=collect(Y_l,_C4>; >

Y_1_C1: =coeff (Y_1 , JU) ; >

Y__1_C2 : =coeff ( Y_1 ,_C2 > ; >

Y_1_C3: «coeff (Y_1 ,_C3 > ; >

Y_1_C4 : =coeff (Y_l , _C4) ;

>A := matrix([[all, *12, аІЗ, al4],[a21, a22, a23, a24],[a31, a32, аЗЗ, a34],[a4l, а42, ad3, a44JJ);

> Ь: =vector ( [Ы, Ь2, ЬЗ, Ь4]); >sollin:=linsolve(А, Ь); > all =Х_1_С1 ; >al2 =Х_1_С2; >al3 =Х_1_СЗ; > al4 =Х_1_С4; >a2l =Y_1_C1; >a22 =Y_1_C2; >a23 =Y_l"c3; >a24 =Y_1_C4; >a31 =Tg_l_Cl ; >a32 =Tg_l_C2 ; > аЗЗ =Tg_l_C3; > a34 =Tg_l_C4 ; >а41 =Tm_l_Cl; >а42 =Tm_l__C2; >а43 =Tm_l_C3 ; >а44 =Tm 1 C4; > Ы : -X__in ; Ь2 := Y_in; b3:=Tg_in; Ь4:=Tm_in; >

CI:=simplifуC2:=siraplify(sollin[2]): >C3:^simplify(sollin[3]): >

C4:«simplify(sollin[4]):

>Cm:»840; Cg:«1040; Vv:=0.2; lambda:=2261000; lambdal:=2261000; av:=70; Gm:=8.77; Gg:«0.125/ >kl:**0.001; >k2:«0; k3:=4e-8;

>X_in:=0.025; Tm_in:=293; Y_in:=0.4; >

Vg:»0.7; >

Via: =0.0037 8;

>bl:»X_in; b2:= Y_in; ЬЗ: =Tg_in; b4:=TmJ.n;

>Tg_in;=472; Vg:=0.7;

>_Cl:=eval>_C2:=eval(C2, 1=0):

>_C3:=eval>_C4:«eval(C4, 1=0): >

X_1; >Y_1' >

Tm_l ; >Tg 1;