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; %esimeses sõlmes
ytuletis(8)=(y(8)-y(7))/h; %viimases sõlmes
for i=2:7
    ytuletis(i)=(y(i+1)-y(i-1))/(2*h); %keskmistatud diferentsvalem
end
    
%lähendamine
xvork=-1:1e-3:-0.3;
ysplain=interp1(x,y,xvork,'spline');
ytul=interp1(x,ytuletis,xvork,'linear');

%joonis
plot(x,y,'r*',xvork,ysplain,x,ytuletis,'k+',xvork,ytul)
legend('algandmed','splain','tuletis','tuletise lähend')
grid on

%defineerima funktsiooni f
f=@(t)interp1(x,ytuletis,t,'linear');
%jooniselt alglähend=-0.7
x_max=fzero(f,-0.7)
y_max=interp1(x,y,x_max,'spline')

