function  I=trapezi_comp(f,a,b,N)
%
% Uso:  I = trapezi_comp(f,a,b,N)
%
% Scopo: Approssimazione di integrali definiti tramite la
%        formula dei trapezi composita
%
% Input:   
%     f    stringa che rappresenta la funzione integranda
%     a,b  intervallo di integrazione
%     N    numero di intervalli della formula composita
%
% Output:
%     I    approssimazione dell'integrale

h=(b-a)/N;
x=linspace(a,b,N+1);

y=eval(f);
I=2*sum(y(2:N));
I=h*(y(1)+I+y(N+1))/2;