Курсовая работа: Интеграционный метод Эйлера для решения линейных систем алгебраических уравнений

while E>Edop

K1=Fun(t,x);

dx=h*K1;

x0=x;

x1=x0+(h/2)*Fun(t+h,x);

x1=x1+(h/2)*Fun(t+h,x1);

x=x+dx;

E=abs(x1-x);

t=t+h;

t_out=[t_out,t];

y_out=[y_out,x];

end

Текст программы для решения ОДУ методом Эйлера с переменным шагом:

function [t_out,y_out]=RungeKutta1(t0,x0,h,Edop);

% функция решения методом Рунге-Кутта 1

t=t0;

hmax=h;

x=x0;

xmax=max(x0)

t_out=t;

y_out=x0;

E=[1;1];

while E>Edop

K1=Fun(t,x);

dx=h*K1;

x0=x;

x1=x0+(h/2)*Fun(t+h,x);

x1=x1+(h/2)*Fun(t+h,x1);

x=x+dx;