
function I=formula_midpoint_composta(f,a,b,n)

% OGGETTO:
% Formula midpoint composta (n subintervalli) per calcolare 
% un'approssimazione "I" di "integrale(f,a,b)".

I=0;
v=linspace(a,b,n+1);
v0=v(1:end-1);
v1=v(2:end);
x=v0+v1;
x=x/2;
fx=feval(f,x);
I=((b-a)/n)*sum(fx);