y_oige=2*exp(3)-1 %vaja y väärtus kohal x=3

x0=0; %algtingimusega määratud x
y0=1; %algtingimusega määratud y
xn=3; %koht, kus tahame leida y väärtuse
n=1000; %võrgupunktide arv, aitab kaasa lahendi täpsusele
h=(xn-x0)/n; %sammupikkus, näitab kui tihedalt x juurde võetakse
f=@(x,y)y+1; %DV parem pool

x(1)=0; %algtingimusest
for i=1:n
    x(i+1)=x(i)+h;
end

y(1)=1; %algtingimusest
for i=1:n
    y(i+1)=y(i)+h*f(x(i),y(i)); %Euleri meetod
end
lahendE=y(n+1)
vigaE=lahendE-y_oige

%lahendame ka prognoosi-korrektsiooni meetodiga
for i=1:n
    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(n+1)
vigaPG=lahendPK-y_oige
