function [s]=cavsimp(a,b,n,f)
% Cavalieri Simpson - si prendono 2n punti
delta=(b-a)/(2*n);
s=f(a)+f(b);
s1=0;
s2=0;
x=a;
for k=1:n-1
    x=x+delta;
    s1=s1+f(x);
    x=x+delta;
    s2=s2+f(x);
end
 x=x+delta;
 s1=s1+f(x);
s=(s+4*s1+2*s2)*delta/3;
end
