naginterfaces.library.quad.dim1_​data

naginterfaces.library.quad.dim1_data(x, y)

dim1_data integrates a function which is specified numerically at four or more points, over the whole of its specified range, using third-order finite difference formulae with error estimates, according to a method due to Gill and Miller (1972).

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

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

Parameters

xfloat, array-like, shape $\left(n\right)$

The values of the independent variable, i.e., the $x_1,x_2,\dots ,x_n$.

yfloat, array-like, shape $\left(n\right)$

The values of the dependent variable $y_{\text{i}}$ at the points $x_{\text{i}}$, for $\text{i}=1,2,\dots ,n$.

Returns

ansfloat

The estimated value of the integral.

erfloat

An estimate of the uncertainty in $ans$.

Raises

NagValueError

(errno $1$)

On entry, $n=⟨value⟩$.

Constraint: $n\ge 4$.

(errno $2$)

On entry, $x\left[0\right]=⟨value⟩$ and $x\left[1\right]=⟨value⟩$, $x[I-2]=⟨value⟩$, $x[I-1]=⟨value⟩$ and $I=⟨value⟩$.

Constraint: either $x\left[0\right]<x\left[1\right]<\cdots <x[n-1]$ or $x\left[0\right]>x\left[1\right]>\cdots >x[n-1]$.

(errno $3$)

On entry, $I=⟨value⟩$, $x[I-1]=⟨value⟩$ and $x[I-2]=⟨value⟩$.

Constraint: $x[I-1]\ne x[I-2]$.

Notes

dim1_data evaluates the definite integral

$$I={\int }_{x_1}^{x_n}y\left(x\right)dx\text{,}$$

where the function $y$ is specified at the $n$-points $x_1,x_2,\dots ,x_n$, which should be all distinct, and in either ascending or descending order. The integral between successive points is calculated by a four-point finite difference formula centred on the interval concerned, except in the case of the first and last intervals, where four-point forward and backward difference formulae respectively are employed. If $n$ is less than $4$, the function fails. An approximation to the truncation error is integrated and added to the result. It is also returned separately to give an estimate of the uncertainty in the result. The method is due to Gill and Miller (1972).

References

Gill, P E and Miller, G F, 1972, An algorithm for the integration of unequally spaced data, Comput. J. (15), 80–83