% Harjutustund 13, ülesanne 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 valemi abil
% siis peab olema funktsioon valemi kujul, y=f(x)
% lähendame näiteks diskreetsel kujul olevat funktsiooni kuupsplainiga
% S^(3,2)(x)
y3splain=interp1(x,y,x,'spline'); % S^(3,2)(x)
hold on
plot(x,y3splain)
% esitame splaini funktsioonina
f=@(z) interp1(x,y,z,'spline') % z-argument
% integraal:
integral(f,-1,-0.3)