% Formule di quadratura - TRAPEZI
% 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

% Metodo dei trapezi composito
h=(b-a)/N ;
xi=a:h:b;
fTR=libreria_funzioni_quad(xi,label) ;
Iapp=0.5*h*(fTR(1)+fTR(end))+h*sum(fTR(2:N)) ;
disp(Iapp)
figure
title(['Approssimazione dell''integrale con metodo dei Trapezi - ' num2str(Iapp)])
plot(x,y,'r--')
grid
hold on
area(xi,fTR)