
% Suddivido [a,b] in N intervalli
% Calcolo del metodo di CS in ciascun intervallo iesimo
h=(b-a)/N ;
xi=a:h:b ;
PM=a+h/2:h:b ;
%
fi=libreria_funzioni_quad(xi,label) ;
fPM=libreria_funzioni_quad(PM,label) ;
h=h/6 ;
Iapp=0;
for i=1:N
    Iapp=Iapp+h*(fi(i)+4*fPM(i)+fi(i+1)) ;
end
disp(Iapp)

plot(x,y,'r--')
grid
title(['Approssimazione dell''Integrale con Cavalieri Simpson - ' num2str(Iapp)])