clear
x=-1:0.1:-0.3;
y=[10 11 11.5 11.6 11.4 11.1 10 8.7];

vastus1=trapz(x,y)
f=@(t)interp1(x,y,t,'spline');
vastus2=quadv(f,-1,-0.3)