% 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')
% lähendame diskreetsel kujul olevat funktsiooni kuupsplainiga 
y3splain=interp1(x,y,x,'spline')  % S^(3,2)(x) 
hold on 
plot(x,y3splain)
% esitame kuupsplaini funktsioonina
f=@(z) interp1(x,y,z,'spline'); % S^(3,2)(z)
f(-0.87) % S^(3,2)(-0.87)=11.1966 % lihtsalt näide

% diskreetsel kujul olev funktsioon on nüüd funktsioonina kirja pandud
% seega saame leida integraali Simpsoni valemiga
integral(f,-1,-0.3)