Посмотрите пожалуйста тут где может быт ошибка
Код:
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./pi);
h=0.1;
x0=1.7;
y0=5.3;
[X,Y]=Eiler_13(y0,x0,x1,h);
i=1:length(X);
x(i)=y0+1/2*X(i).^2+pi *sin(Y9i)/pi);
plot(X,Z,Y)