clear
x=1:0.5:5;
y=[0 10 18 25 32.2 39 43 48 52.3];

pindala1=trapz(x,y) %trapetsvalemiga

%Simpsoni valemiga
%erinevad võimalused f(x) asendada
f=@(z)interp1(x,y,z,'spline'); %kuupsplain
f1=@(w)interp1(x,y,w,'linear'); %lineaarsplain
f2=@(v)interp1(x,y,v,'nearest'); %lähim täisarvuline lähend

pindala2=quadv(f,1,5)
pindala3=quadv(f1,1,5)
pindala4=quadv(f2,1,5)
