clear
f=@(x,y)x.^2-y.^2;
h=0.02;
x=-1:h:2;
n=(2-(-1))/h;
y(1)=0;
%prognoosi-korrektsiooni meetod sammuga 0.02
for i=1:n
    y(i+1)=y(i)+0.5*h*f(x(i),y(i))+0.5*h*f(x(i+1),y(i)+h*f(x(i),y(i)));
end
vastus=y(n)