clear
%lähteandmed
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;
ytuletis(8)=(y(8)-y(7))/h;
for i=2:7
    ytuletis(i)=(y(i+1)-y(i-1))/(2*h);
end
ytuletis

%lähendamine
xvork=-1:0.0001:-0.3;
ykuup=interp1(x,y,xvork,'spline');
ytullin=interp1(x,ytuletis,xvork,'linear');

%joonis
plot(x,y,'r*',xvork,ykuup,xvork,ytullin)
grid on
%xmax on ligikaudu -0.7
f=@(z)interp1(x,ytuletis,z,'linear');
xmax=fzero(f,-0.7)
ymax=interp1(x,y,xmax,'spline')