clear
x=-1:0.1:-0.3;
h=0.1;
y=[10 11 11.5 11.6 11.4 11.1 10 8.7];

%tuletis
ytuletis(1)=(y(2)-y(1))/h; %diferentsvalem sammuga ette
ytuletis(8)=(y(8)-y(7))/h;
for i=2:7
    ytuletis(i)=(y(i+1)-y(i-1))/(2*h);
end
ytuletis

%interpoleerimine
vork=-1:1e-3:-0.3;
ykuup=interp1(x,y,vork,'spline'); %kuupsplain
ytullin=interp1(x,ytuletis,vork,'linear'); %lineaarsplain

%joonis
plot(x,y,'r*',vork,ykuup,vork,ytullin)
legend('lähteandmed','kuupsplain','tuletise lähend')
grid on
f=@(t)interp1(x,ytuletis,t,'linear');
%alglähend -0,7 jooniselt
xmax=fzero(f,-0.7)
ymax=interp1(x,y,xmax,'spline')