контрольные работы 7 вариант / 13 / Eiler_13

function [X,Y]=Eiler_13(y0,x0,x1,h)
N=(x1-x0)/h;
x(1)=x0;
y(1)=y0;
for i=1:N
x(i+1)=x(1)+h*i;
y(i+1)=y(i)+h*F13(x(i),y(i));
end;
X=x;
Y=y;

function z=F13(x,y)
z=x+cos(y/exp(1));