
function S=composta_simpson3ottavi(f,estremi,N)

a=estremi(1);
b=estremi(2);
t=linspace(a,b,N);
S=0;
for k=1:N-1
    aL=t(k); bL=t(k+1);
    estremiL=[aL bL];
    SL=simpson3ottavi(f,estremiL);
    S=S+SL;
end
