function I = trapezi(fun,a,b,m)

% I = trapezi(fun,a,b,m)
%
% Applica il metodo dei trapezi  
% per calcolare una approssimazione dell'integrale 
% tra a e b di fun suddividendo [a,b] in 
% m sottointervalli di uguale ampiezza.
% Si assume (b>=a).

h = (b-a)/m; % ampiezza dei sottointervalli

x = a + h*[0:m];
f = feval(fun,x);

I = (h/2)*(f(1) + 2*sum(f(2:end-1)) + f(end));

return

