% Tund 13, ülesanne 3
x=[-1:0.1:-0.3];
y=[10,11,11.5,11.6,11.4,11.1,10,8.7];
plot(x,y,'o')
hold on
trapz(x,y)
% kui h->0, siis integraali väärtus täpsem
% interpoleerime andmeid kuupsplainiga S^(3,2)(x)
xvork=[-1:0.001:-0.3];
y3splain=interp1(x,y,xvork,'spline');
plot(xvork,y3splain,'b-')
% esitame kuupsplaini funktsioonina
f=@(z) interp1(x,y,z,'spline');
integral(f,-1,-0.3)
hold off