NAG Library Function Document
nag_1d_quad_vals (d01gac) integrates a function which is specified numerically at four or more points, over the whole of its specified interval using third-order finite difference formulae with error estimates, according to a method due to Gill and Miller.
||nag_1d_quad_vals (Integer n,
const double x,
const double y,
nag_1d_quad_vals (d01gac) evaluates the definite integral
where the function
is specified at the
, 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
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)
Gill P E and Miller G F (1972) An algorithm for the integration of unequally spaced data Comput. J. 15 80–83
n – IntegerInput
On entry: the number of points, .
x[n] – const doubleInput
On entry: the values of the independent variable, i.e., .
either or .
y[n] – const doubleInput
On entry: the values of the dependent variable at the points
, for .
ans – double *Output
On exit: the estimate of the integral.
er – double *Output
: an estimate of the uncertainty in ans
fail – NagError *Input/Output
The NAG error argument (see Section 3.6
in the Essential Introduction).
6 Error Indicators and Warnings
- No error is reported arising from the relative magnitudes of ans and er on return, due to the difficulty when the true answer is zero.
On entry, .
The sequence x
is not strictly decreasing:
The sequence x
is not strictly increasing:
The sequence x
has first two points equal:
No accuracy level is specified before calling nag_1d_quad_vals (d01gac) but on return
is an approximation to, but not necessarily a bound for,
. If on exit fail
is not equal to NE_NOERROR, both ans
are returned as zero.
The time taken by nag_1d_quad_vals (d01gac) depends on the number of points supplied, .
In their paper, Gill and Miller (1972)
do not add the quantity er
before return. However, extensive tests have shown that a dramatic reduction in the error often results from such addition. In other cases, it does not make an improvement, but these tend to be cases of low accuracy in which the modified answer is not significantly inferior to the unmodified one. You have the option of recovering the Gill–Miller answer by subtracting er
on return from the function.
This example evaluates the integral
reading in the function values at 21 unequally-spaced points.
9.1 Program Text
Program Text (d01gace.c)
9.2 Program Data
Program Data (d01gace.d)
9.3 Program Results
Program Results (d01gace.r)