Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
13
Добавлен:
23.11.2018
Размер:
389 б
Скачать
function [X,Y]= RungeKutt_13_3(y0,x0,h,N)
x(1)=x0;
y(1)=y0;
for i=2:N
x(i)=x(1)+h*(i-1);
k1=h*F(x(i-1),y(i-1));
k2=h*F(x(i-1)+h/2,y(i-1)+k1/2);
k3=h*F(x(i-1)+h/2,y(i-1)+k2/2);
k4=h*F(x(i-1)+h,y(i-1)+k3);
y(i)=y(i-1)+1/6*(k1+2*k2+2*k3+k4);
end;
X=x;
Y=y;

function z=F(x,y)
z=x+cos(y./exp(1));
Соседние файлы в папке 13