% Ylesanded 14, ylesanne 1
syms yy(x);
vorrand=diff(yy,x)==yy-2*x;
% dif. vorrandi lahend algtingimusel y(0)=1
lahend=matlabFunction(dsolve(vorrand,yy(0)==1))
fplot(lahend,[-1,4])
% tuleb leida y(2)=?
tapnelahend=lahend(2)

% lahendame Euleri meetodiga
h=0.4;
x0=0;
xn=2;
n=(xn-x0)/h  % osaloikude arv
x(1)=0; % algtingimuse x-i vaartus
y(1)=1; % algtingimuse y-i vaartus
f=@(x,y) y-2*x;

for i=2:(n+1)
  x(i)=x(i-1)+h;
  y(i)=y(i-1)+h*f(x(i-1),y(i-1));
end
x
y
yeuler=y(n+1) % kui x=2, siis y¬0.6218

% joonisel
hold on
plot(x,y,'o-')
