This section introduces some of the mathematical concepts used with PV‑WAVE.

The first function in this chapter, INTFCN, is designed to compute approximations to integrals of the following form:

or:

The weight function w is used to incorporate known singularities (either algebraic or logarithmic) or to incorporate oscillations. The default action of this function assumes univariate quadrature, a weight function w(x) = 1, and the existence of endpoint singularities. Even if no endpoint singularities exist, the default method is still effective for general-purpose integration. If more efficiency is desired, then a more specialized method can be specified through the use of specific parameter and keyword combinations. The available methods can be summarized as follows, where the description refers to subsections of the documentation for the function INTFCN:

Function INTFCN returns an estimated answer R and provides keywords to specify a requested absolute error ε, the requested relative error ρ, and a named variable in which to return an estimate of the error E. These numbers are related in the equation:

One situation that arises in univariate quadrature concerns the approximation of integrals when only tabular data is given. The functions above do not directly address this question. However, the standard method for handling this problem is to interpolate the data then integrate the interpolant. This can be accomplished by using the PV-WAVE IMSL Mathematics spline interpolation functions with the spline integration function, SPINTEG Function on page 134, which can be found in Chapter 4: Interpolation and Approximation.

Two functions, INTFCN and INTFCNHYPER, have been included in this chapter that can be used to approximate certain multivariate integrals.

INTFCN can be called with additional parameters and keywords to return an approximation to a two-dimensional iterated integral of the form:

Function INTFCNHYPER returns an approximation to the integral of a function of n variables over a hyper-rectangle as shown in the equation:

When working with two-dimensional, tensor-product tabular data, use the PV-WAVE IMSL Mathematics spline interpolation function BSINTERP, followed by the spline integration function SPINTEG.

For a fixed number of nodes, N, the Gauss quadrature rule is the unique rule that integrates polynomials of degree less than 2N. These quadrature rules can be easily computed using procedure GQUAD, which produces the points {xi} and weights {wi} for i = 1, ..., N that satisfy:

for all functions f that are polynomials of degree less than 2N. The weight functions w can be selected from Table 5-1: Weight Functions:

Where permissible, GQUAD also computes Gauss-Radau and Gauss-Lobatto quadrature rules.