`quad('expression',ll,ul)`

is modelled after the Octave/Matlab integration
function, but much simpler. Simpson's method is applied with a fixed number of
nodes. This function uses the ``eval''-method rather than symbolic variables.
The function has to be supplied as quoted string, and must be compatible
with vector arguments. Finally, the variablename
must be `x`

. `romberg`

uses symbolic function
definitions, and a symbolic variable has to be supplied.
The maximum number of iterations is set by the variable `rombergit`

(default ) and accuracy by `rombergtol`

(default: ).
>> quad('exp(-x.^2)',0,5) s = 0.88623 >> syms x >> romberg(exp(-x^2),x,0,5) ans = 0.88623

2009-03-15