clear
x=1:0.5:5;
y=[0 10 18 25 32.2 39 43 48 52.3];
%pindala trapetsvalemiga
pindala1=trapz(x,y)
%pindala Simpsoni valemiga
f=@(z)interp1(x,y,z,'spline');
pindala2=quadv(f,1,5)