Rombint.m

This function performs the numerical evaluation of an integral using the Romberg method.

This function performs the numerical evaluation of an integral using the Romberg method.

Q = rombint('F',A,B) approximates the integral of F(X) from A to B to within a relative error of 1e-10 using Romberg's method of integration. 'F' is a string containing the name of the function. The function must return a vector of output values if a vector of input values is given.

Q = rombint('F',A,B,DECIMALDIGITS) integrates with accuracy 10^{-DECIMALDIGITS}.

Q = rombint('F',A,B,DECIMALDIGITS,P1,P2,...) allows coefficients P1, P2, ... to be passed directly to function F: G = F(X,P1,P2,...).

Created under MATLAB 5.2, works with all subsequent releases.