232
-->y=ode(x0,t0,t,syst);
-->//Формирование графического решения
-->plot(t,y);
-->xgrid();
Рис. 25.2. Решение задачи 25.2.
Пример 25.3.
Найти решение задачи Коши для следующей жесткой системы:
;
199
8. 112
75. 181
58. 115
211 .54 182 .63
932 .85
302 .53
577 .8 636 .3
136 .10
395 .10
03. 121
88. 126
38. 185
46. 119
X
dt
dX
−
−
−
−
−
−
−
−
−
−
=
( )
=
1
1
1
1
0
X
.
Решение системы:
-->B=[119.46 185.38 126.88 121.03;-10.395 -10.136 -3.636 8.577;
-->-53.302 -85.932 -63.182 -54.211;-115.58 -181.75 -112.8 -199];
-->function dx=syst1(t,x), dx=B*x,endfunction
-->function J=Jac(t,y),J=B,endfunction
-->x0=[1;1;1;1]; t0=0; t=0:0.01:5;
-->y=ode("stiff",x0,t0,t,syst1,Jac);
-->plot(t,y);
-->xgrid();
Графическое решение показано на рис. 25.3.
Рис. 25.3. Графическое решение задачи 25.3.
