clear
x=-1:0.1:-0.3;
y=[10 11 11.5 11.6 11.4 11.1 10 8.7];
h=0.1;
%tuletised
ytuletis(1)=(y(2)-y(1))/h;
ytuletis(8)=(y(8)-y(7))/h;
for i=2:7
    ytuletis(i)=(y(i+1)-y(i-1))/(2*h);
end
ytuletis

%interpoleerimine
xvork=-1:0.001:-0.3;
ysplain=interp1(x,y,xvork,'spline');
ytuletislin=interp1(x,ytuletis,xvork,'linear');
%joonis
plot(x,y,'r*',xvork,ysplain,xvork,ytuletislin)
grid on
legend('algandmed','kuupsplain','tuletis')
%leiame tuletise nullkoha
f=@(s)interp1(x,ytuletis,s,'linear'); %saime lineaarsplaini võrrandi
%jooniselt näeme, et võimalik nullkoht on -0,7 juures
xmax=fsolve(f,-0.7)
%leiame ymax
ymax=interp1(x,y,xmax,'spline')