clear
x=-1:0.1:-0.3;
y=[10 11 11.5 11.6 11.4 11.1 10 8.7];
%integraal trapetsvalemiga
inttrapets=trapz(x,y)
%integraal Simpsoni valemiga
f=@(z)interp1(x,y,z,'spline');
intSimpson=quadv(f,-1,-0.3)