l>
%
% * ErrBd  REPLACE_WITH_DASH_DASH- estimated integral error | I - Q |
%
% * optParams REPLACE_WITH_DASH_DASH- optional parameters useful to debug and get better
%                  understanding of the algorithm
%
% * optParams.aMLEAll REPLACE_WITH_DASH_DASH- returns the shape parameters computed
%
%
%%  Guarantee
%
% This algorithm attempts to calculate the integral of function f over the
% hyperbox $[0,1]^{\mathrm{dim}}$ to a prescribed error tolerance 
% tolfun:= max(abstol,reltol*| I |) with guarant