
function I=formula_rettangolo_composta(f,a,b,m)

% OGGETTO:
% Formula composta del rettangolo ("m" subintervalli, valutata nell'estremo
% a sinistra) per calcolare un'approssimazione "I" di "integrale(f,a,b)".

if m <= 0
    I=NaN;
    fprintf('\n \t m non positivo');
else
    I=0;
    v=linspace(a,b,m+1);
    v0=v(1:m);
    fv0=feval(f,v0);
    I=((b-a)/m)*sum(fv0);
end

 
