clear
%algandmed ülesande tekstist
x0=2
xn=10;
n=1000; %meie valitud väärtus
h=(xn-x0)/n;
f=@(x,y)cos(x*y)
%teeme vektorid x ja y
x(1)=2; %vektori x esimene element
y(1)=2; %vektori y esimene element

%Euleri meetod
for i=1:n
    y(i+1)=y(i)+h*f(x(i),y(i));
    x(i+1)=x(i)+h;
end
lahendEuler=y(i+1)
%graafik
plot(x,y)
hold on

%prognoosi-korrektsiooni meetod
for i=1:n
    x(i+1)=x(i)+h;
    y(i+1)=y(i)+h/2*f(x(i),y(i))+h/2*f(x(i+1),y(i)+h*f(x(i),y(i)));
end
lahendPK=y(i+1)
plot(x,y,'r')
legend('Euleri meetod','Prognoosi-korrektsiooni meetod')