
function demo

for n=10:10:100
    f=@(x) sin(x); a=0; b=pi/2;
    I(n)=formula_midpoint_composta(f,a,b,n);
    fprintf('\n \t n: %3.0f I(n): %1.15e',n,I(n));
end

fprintf('\n');

% >> demo  
%    n:  10 I(n): 1.001028824142709e+00
%  	 n:  20 I(n): 1.000257067197303e+00
%  	 n:  30 I(n): 1.000114240667270e+00
%  	 n:  40 I(n): 1.000064258127219e+00
%  	 n:  50 I(n): 1.000041124535493e+00
%  	 n:  60 I(n): 1.000028558454002e+00
%  	 n:  70 I(n): 1.000020981610028e+00
%  	 n:  80 I(n): 1.000016063989881e+00
%  	 n:  90 I(n): 1.000012692505260e+00
%  	 n: 100 I(n): 1.000010280911905e+00
% 
% >> format long e
%  -cos(pi/2)-(-cos(0)) % valore esatto 
%     ans =
%       9.999999999999999e-01
% >> % Il valore esatto e' 1.