Function Reference
— Function File: q = quadl (f, a, b)
— Function File: q = quadl (f, a, b, tol)
— Function File: q = quadl (f, a, b, tol, trace)
— Function File: q = quadl (f, a, b, tol, trace, p1, p2, ...)

Numerically evaluate integral using adaptive Lobatto rule. quadl (f, a, b) approximates the integral of f(x) to machine precision. f is either a function handle, inline function or string containing the name of the function to evaluate. The function f must return a vector of output values if given a vector of input values.

If defined, tol defines the relative tolerance to which to which to integrate f(x). While if trace is defined, displays the left end point of the current interval, the interval length, and the partial integral.

Additional arguments p1, etc, are passed directly to f. To use default values for tol and trace, one may pass empty matrices.

Reference: W. Gander and W. Gautschi, 'Adaptive Quadrature - Revisited', BIT Vol. 40, No. 1, March 2000, pp. 84–101. http://www.inf.ethz.ch/personal/gander/