контрольные работы 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));