
function S=simpson3ottavi(f,estremi)

a=estremi(1);
b=estremi(2);
h=(b-a)/3;
x=linspace(a,b,4);
y=feval(f,x);
S=(3*h/8)*(y(1)+3*y(2)+3*y(3)+y(4));