clear
x=1:0.1:2;
y=[4 7 2.1 15 3 5.2 6 8.5 14 15.1 17.8];
h=0.1; %sammupikkus

%tuletis
ytuletis(1)=(y(2)-y(1))/h; %diferentsvalem sammuga ette
ytuletis(11)=(y(11)-y(10))/h; %diferentsvalem sammuga taha
%ülejäänud tuletised leiame keskmistatud diferentsvalemiga
for i=2:10
    ytuletis(i)=(y(i+1)-y(i-1))/(2*h);
    yIItuletis(i)=(y(i-1)-2*y(i)+y(i+1))/h^2;
end
ytuletis
yIItuletis

%uurime funktsiooni tuletist kasutades
xvork=1:0.001:2;
%lähendame tuletist
ytullineaar=interp1(x,ytuletis,xvork,'linear');
%joonis
plot(x,y,'k*',xvork,ytullineaar)
grid on