Purpose¶

Integrates a specified function using Gauss-Legendre quadrature. A suite of upper and lower bounds may be calculated in one procedure call.

Format¶

y = intquad2(&f, xlims, ylims[, ...])
Parameters: &f (scalar) – pointer to the procedure containing the function to be integrated. xlims (2x1 or 2xN matrix) – the limits of x. ylims (2x1 or 2xN matrix) – the limits of y. .. (any) – Optional. A variable number of extra scalar arguments to pass to the user function. These arguments will be passed to the user function untouched. y (Nx1 vector) – the estimated integral(s) of $$f(x,y)$$ evaluated between the limits given by xl and yl.

Global Input¶

_intord

scalar, the order of the integration. The larger _intord, the more precise the final result will be. _intord may be set to 2, 3, 4, 6, 8, 12, 16, 20, 24, 32, 40.

Default = 12.

Examples¶

Basic example¶

// Define function to be integrated
proc f(x, y);
retp(x .* sin(x + y));
endp;

// Limits of integration
xlims = { 1, 0 };
ylims = { 1, 0 };

// Calculate integral


After the code above, ans should equal:

0.42892501


Multiple integration limits¶

// Define function to be integrated
proc (1) = myProc(x, y);
retp(x .* sin(x + y));
endp;

// Define multiple integration limits
xlims = {  1  0.5,
0.5    0 };

ylims = {  1  0.5,
0.5  0.3 };

// Calculate integrals


This will integrate the function:

myProc(x) = x.*sin(x+y)


between x = 0 and 0.5, and between y = 0.3 and 0.5 as well as between x = 0.5 and 1, and between y = 0.5 and 1.

The returned variable, ans should be equal to:

0.18352849
0.016593029


Extra arguments to function¶

// Define function to be integrated that takes an additional argument
proc f(x, y, a);
retp(x .* sin(a .* x + y));
endp;

// Limits of integration
xlims = { 1, 0 };
ylims = { 1, 0 };

// Assign extra scalar argument
a = pi/2;

// Calculate integral
ans = intquad2(&f, xlims, ylims, a);


After the code above, ans should equal:

0.44737953


Remarks¶

The user-defined function f must return a vector of function values. intquad2() will pass to user-defined functions a vector or matrix for x and y and expect a vector or matrix to be returned. Use .* and ./ instead of * and /.

intquad2() will expand scalars to the appropriate size. This means that functions can be defined to return a scalar constant. If users write their functions incorrectly (using * instead of .*, for example), intquad2() may not compute the expected integral, but the integral of a constant function.

To integrate over a region which is bounded by functions, rather than just scalars, use intgrat2() or intgrat3().

integral.src

Globals¶

_intord, _intq12, _intq16, _intq2, _intq20, _intq24, _intq3, _intq32, _intq4, _intq40, _intq6, _intq8