% Harjutustund 13, ylesanne 3
x=[-1:0.1:-0.3];
y=[10,11,11.5,11.6,11.4,11.1,10,8.7];
trapz(x,y)
plot(x,y,'o')

% leiame integraali ka Simpsoni valemiga näiteks
% selleks peab funktsioon olema valemi kujul
% esitame diskreetsel kujul oleva funktsiooni kuupsplainina
y3splain=interp1(x,y,x,'spline')  
% S^(3,2)(x) splaini/funktsiooni väärtused x-de korral
hold on
plot(x,y3splain)
% esitame splaini funktsiooni kujul
f=@(z) interp1(x,y,z,'spline')  % S^(3,2)(z) % argument z
integral(f,-1,-0.3)