clear
f=@(x,y)sin(x.^2)./(y+3);
h=0.1;
x=1:h:4;
n=3/h;
y(1)=3; 
%Euleri meetod sammuga 0.1
for i=1:n
    y(i+1)=y(i)+h*f(x(i),y(i));
end
y
%interpoleerime lahendit kuupsplainiga
xvork=1:1e-3:4;
ysplain=interp1(x,y,xvork,'spline');

%kasutame käsku ode45
ode45(f,[1,4],3)
hold on %hoiab joonise alles
plot(xvork,ysplain)