function ret = rettangoli (fun, n, a, b, x0)
%
% ret = rettangoli (fun, n, a, b, x0)
%
% Formula dei rettangoli a n+1 punti con x0 nel primo intervallo.
h = (b - a) / n;
i = (1:n)';
x = x0 + (i - 1) * (b - a) / n;
ret = (b - a) / n * sum (fun (x));
%!demo
%! fun = @(x) x .^ 2;
%! a = 0;
%! b = 1;
%! nrange = 10:10:50;
%! counter = 0;
%! for n = nrange
%!   counter = counter + 1;
%!   ret = rettangoli (fun, n, a, b, a);
%!   err(counter) = abs (1 / 3 - ret);
%! end
%! loglog (nrange, err, '*', nrange, nrange .^ (-1))
%! axis([10,50,1e-5,1e-1])
%!demo
%! fun = @(x) x .^ 2;
%! a = 0;
%! b = 1;
%! nrange = 10:10:50;
%! counter = 0;
%! for n = nrange
%!   counter = counter + 1;
%!   ret = rettangoli (fun, n, a, b, a + (b - a) / n / 2);
%!   err(counter) = abs (1 / 3 - ret);
%! end
%! loglog (nrange, err, '*', nrange, nrange .^ (-2))
%! axis([10,50,1e-5,1e-1])
