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

vastus_trapets=trapz(x,y) %vastus trapetsvalemiga

%leiame f(x) avaldise
f=@(z)interp1(x,y,z,'spline'); %f(x) on kuupsplain
f1=@(v)interp1(x,y,v,'linear'); %f(x) on lineaarsplain
f2=@(w)interp1(x,y,w,'nearest'); %f(x) on lähim täisarvuline lähend
vastus_simpson=quadv(f,1,5)
vastus_simpson1=quadv(f1,1,5)
vastus_simpson2=quadv(f2,1,5)