clear
x=0:6;
F=[0 5.1 13.2 16.1 25.4 40.2 61.9];
for i=2:7
    F(i)=F(i)+(-1)^(i+5)*0.1*5;
end
vastus1=trapz(x,F)
f=@(t)interp1(x,F,t,'spline');
vastus2=quadv(f,0,6)

    