naginterfaces.library.quad.dim2_​fin

naginterfaces.library.quad.dim2_fin(ya, yb, phi1, phi2, f, absacc, data=None)

dim2_fin attempts to evaluate a double integral to a specified absolute accuracy by repeated applications of the method described by Patterson (1968) and Patterson (1969).

For full information please refer to the NAG Library document for d01da

https://support.nag.com/numeric/nl/nagdoc_30/flhtml/d01/d01daf.html

Parameters

yafloat

$a$, the lower limit of the integral.

ybfloat

$b$, the upper limit of the integral. It is not necessary that $a<b$.

phi1callable retval = phi1(y, data=None)

$phi1$ must return the lower limit of the inner integral for a given value of $y$.

Parameters

yfloat

The value of $y$ for which the lower limit must be evaluated.

dataarbitrary, optional, modifiable in place

User-communication data for callback functions.

Returns

retvalfloat

The lower limit of the inner integral for the given value of $y$.

phi2callable retval = phi2(y, data=None)

$phi2$ must return the upper limit of the inner integral for a given value of $y$.

Parameters

yfloat

The value of $y$ for which the upper limit must be evaluated.

dataarbitrary, optional, modifiable in place

User-communication data for callback functions.

Returns

retvalfloat

The upper limit of the inner integral for the given value of $y$.

fcallable retval = f(x, y, data=None)

$f$ must return the value of the integrand $f$ at a given point.

Parameters

xfloat

The coordinates of the point $(x,y)$ at which the integrand $f$ must be evaluated.

yfloat

The coordinates of the point $(x,y)$ at which the integrand $f$ must be evaluated.

dataarbitrary, optional, modifiable in place

User-communication data for callback functions.

Returns

retvalfloat

The value of $f\left(x\right)$ evaluated at $x$.

absaccfloat

The absolute accuracy requested.

dataarbitrary, optional

User-communication data for callback functions.

Returns

ansfloat

The estimated value of the integral.

nptsint

The total number of function evaluations.

Raises

NagValueError

(errno $imod10\ne 0$)

The outer integral has not converged, and $n$ of the inner integrals have not converged with $255$ points: $n=⟨value⟩$. $ans$ may still contain an approximate estimate of the integral, but its reliability will decrease as $n$ increases.

(errno $imod10=0$)

The outer integral has converged, but $n$ of the inner integrals have not converged with $255$ points: $n=⟨value⟩$. $ans$ may still contain an approximate estimate of the integral, but its reliability will decrease as $n$ increases.

Warns

NagAlgorithmicWarning

(errno $1$)

The outer integral has not converged with $255$ points. However, all the inner integrals have converged, and $ans$ may still contain an approximate estimate of the integral.

Notes

dim2_fin attempts to evaluate a definite integral of the form

$$I={\int }_a^b{\int }_{{\varphi }_1\left(y\right)}^{{\varphi }_2\left(y\right)}f(x,y)dxdy$$

where $a$ and $b$ are constants and ${\varphi }_1\left(y\right)$ and ${\varphi }_2\left(y\right)$ are functions of the variable $y$.

The integral is evaluated by expressing it as

$$I={\int }_a^{b}F\left(y\right)dy\text{, where}F\left(y\right)={\int }_{{\varphi }_1\left(y\right)}^{{\varphi }_2\left(y\right)}f(x,y)dx\text{.}$$

Both the outer integral $I$ and the inner integrals $F\left(y\right)$ are evaluated by the method, described by Patterson (1968) and Patterson (1969), of the optimum addition of points to Gauss quadrature formulae.

This method uses a family of interlacing common point formulae. Beginning with the $3$-point Gauss rule, formulae using $7$, $15$, $31$, $63$, $127$ and finally $255$ points are derived. Each new formula contains all the points of the earlier formulae so that no function evaluations are wasted. Each integral is evaluated by applying these formulae successively until two results are obtained which differ by less than the specified absolute accuracy.

References

Patterson, T N L, 1968, On some Gauss and Lobatto based integration formulae, Math. Comput. (22), 877–881

Patterson, T N L, 1969, The optimum addition of points to quadrature formulae, errata, Math. Comput. (23), 892