% Harjutustund 14, ülesanne 7
x=[1:0.5:5];
y=[0,10,18,25,32.2,39,43,48,52.3];
plot(x,y,'o')
trapz(x,y)
% interpoleerime esialgset funktsiooni kuupsplainiga
xx=[1:0.01:5];
ys3=interp1(x,y,xx,'spline');
hold on
plot(xx,ys3)
% esitame kuupsplaini funktsioonina
f=@(z) interp1(x,y,z,'spline');
pindala=integral(f,1,5)