% Formule di quadratura - PUNTO MEDIO
% Approssimo l'integrale di f(x) in dx, tra a e b
% a e b inseriti da tastiera
% f(x) definita in una function a parte


% Suddivido [a,b] in N intervalli
% Calcolo del punto medio PM in ciascun intervallo iesimo
h=(b-a)/N ;
PM=a+h/2:h:b ;
% Calcolo di f(PM(i))
%effe=inline('exp(x.^2)') ;
fPM=libreria_funzioni_quad(PM,label) ;
size(fPM)
% Calcolo dell'integrale approssimato
Iapp=sum(fPM*h) ;
disp('Valore approssimato dell''integrale')
disp(Iapp)

plot(x,y,'r--')
grid
hold on
bar(PM,fPM,0.8)
title(['Approssimazione dell''Integrale con Punto Medio - ' num2str(Iapp)])