% -*- texinfo -*-
% @deftypefn {Function File} {@var{q} =} quadv (@var{f}, @var{a}, @var{b})
% @deftypefnx {Function File} {@var{q} =} quadl (@var{f}, @var{a}, @var{b}, @var{tol})
% @deftypefnx {Function File} {@var{q} =} quadl (@var{f}, @var{a}, @var{b}, @var{tol}, @var{trace})
% @deftypefnx {Function File} {@var{q} =} quadl (@var{f}, @var{a}, @var{b}, @var{tol}, @var{trace}, @var{p1}, @var{p2}, @dots{})
% @deftypefnx {Function File} {[@var{q}, @var{fcnt}] =} quadl (@dots{})
%
% Numerically evaluate integral using adaptive Simpson's rule.
% @code{quadv (@var{f}, @var{a}, @var{b})} approximates the integral of
% @code{@var{f}(@var{x})} to the default absolute tolerance of @code{1e-6}.
% @var{f} is either a function handle, inline function or string
% containing the name of the function to evaluate. The function @var{f}
% must accept a string, and can return a vector representing the
% approximation to @var{n} different sub-functions.
%
% If defined, @var{tol} defines the absolute tolerance to which to
% which to integrate each sub-interval of @code{@var{f}(@var{x})}.
% While if @var{trace} is defined, displays the left end point of the
% current interval, the interval length, and the partial integral.
%
% Additional arguments @var{p1}, etc., are passed directly to @var{f}.
% To use default values for @var{tol} and @var{trace}, one may pass
% empty matrices.
% @seealso{triplequad, dblquad, quad, quadl, quadgk, trapz}
% @end deftypefn