Курсовая работа: Численные методы интегрирования и оптимизации сложных систем

i=i+1;

A(i)=x;

end

hold on

for n=1:i

plot(n,A(n),'b-o')

end

grid on

xlabel('iteraciya')

ylabel('roots')

disp('ответ');

disp(x);

8. Приложение 2 (Листинг скриптов для решения ДУ)

function Difer

clc

T=4;

a0=638.89;

a1=56.35;

a2=5.60;

b0=595.24;

h=0.0005;

A_X(1,1:3)=[0 0 0];

A=[0 1 0;

0 0 1;

a0 a1 a2];

B=[0 0 b0]';

k=0;

t=0;

while (t < (T-h))