Quadrature formula

From Encyclopedia of Mathematics
Revision as of 16:59, 7 February 2011 by (talk) (Importing text file)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

An approximate formula for the calculation of a definite integral:


On the left-hand side is the integral to be calculated. The integrand is written as a product of two functions. The first, , is considered to be fixed for the given quadrature formula and is called the weight function; the function belongs to a fairly-wide class of functions, for example, continuous functions for which the integral on the left-hand side of (1) exists. The sum on the right-hand side of (1) is called the quadrature sum, the numbers are called the nodes of the quadrature formula, while the numbers are called its weights. The determination of an approximate value of the integral by means of formula (1) reduces to the calculation of the quadrature sum; the values of the nodes and the weights are usually taken from tables (see, for example, [3]).

The most widespread quadrature formulas are those based on algebraic interpolation. Let be distinct points (usually although this requirement is not essential) and let be the interpolation polynomial for constructed from its values at these points:

Here is the Lagrangian basis polynomial of the -th node (cf. Lagrange interpolation formula): ( is the Kronecker symbol). The integral over of is approximately replaced by the integral of ; one obtains an approximate equation of the form (1), in which


The existence of the integral in (2) is equivalent to the existence of the moments of the weight function:

(here and in what follows it is assumed that the required moments of exist; in particular, in the case the interval is taken to be finite, cf. Moment problem).

The quadrature formula (1) in which the weights are defined by equations (2) is called an interpolatory quadrature formula. An integer is called the algebraic degree of accuracy of (1) if (1) is exact when is any polynomial of degree not exceeding and if it is not exact for . In order that (1) be an interpolatory quadrature formula it is necessary and sufficient that its algebraic degree of accuracy satisfies the inequality .

Let and be finite. The interpolatory quadrature formula with equally-spaced nodes


where is a positive integer, , is called the Newton–Cotes quadrature formula; this quadrature formula has algebraic degree of accuracy when is odd and when is even. The interpolatory quadrature formula with a single node,

is called the rectangle rule or midpoint rule; its algebraic degree of accuracy when and in the remaining cases.



An interpolatory quadrature formula (1) in which the nodes are the roots of an orthogonal polynomial of degree on with weight function is called a quadrature formula of Gauss type; it is also called a quadrature formula of highest algebraic degree of accuracy, since under the conditions (4) no quadrature formula with nodes can be precise for . The most widely used quadrature formulas of Gauss type are those defined by the following special cases of the weight function and interval :

the Jacobi weight () on with parameter values: a) (the Gauss quadrature formula); b) (the Mehler quadrature formula); c) ; and d) ;

the Hermite weight on (); and

the Laguerre weight () on .

There exist quadrature formulas in which some of the nodes are fixed in advance, while the remaining ones are chosen so that the formula will have highest algebraic degree of accuracy. Such, for example, are the Lobatto quadrature formula and the Radau quadrature formula for the calculation of an integral over with weight function 1. In the first of these the nodes are , while in the second, one of these points is fixed only.

Two quadrature formulas with weight function 1,

are said to be similar if , , , where is defined by the equation . In the case of a finite interval ,


where the are defined by (3). If for the calculation of integrals over the intervals one applies quadrature formulas that are similar to the same quadrature formula, equation (5) leads to a composite quadrature formula for the calculation of the integral at the left-hand side. Such, for example, is the composite rectangle rule:

In the case this quadrature formula is exact for for .

It is possible to consider the interpolatory formulas obtained by integrating the Hermite interpolation polynomials of the function . In the quadrature sum of such a quadrature formula not only the values of the function itself at the nodes do enter, but also the values of its successive derivatives up to some order. The values of the derivatives of the integrand at the ends of the range of integration are also used in the Euler–MacLaurin formula.

For the error of the quadrature formula (1),

there are representations in which the derivative enters. These representations are of little use for the actual estimation of , since one needs an estimate of the derivative . The error is an additive homogeneous functional on the vector space of functions for which it is defined.

Another approach to the construction of a quadrature formula is based on minimization of the norm of the error functional .

Let be the error of a quadrature formula that is exact for all polynomials of degree not exceeding , where and . Let ( and is a positive integer) be the vector space of functions that have on an absolutely-continuous derivative of order and a derivative of order whose -th power is summable. Two functions in are considered to be equivalent if their difference is a polynomial of degree not exceeding . The set of equivalence classes (the quotient space of by the vector space of polynomials of degree not exceeding ) is a vector space, denoted by . A norm can be introduced into be setting for a class ,

where is any function belonging to . The error functional of the quadrature formula can be considered on by setting , . The error functional is continuous on the normed linear space . Its norm characterizes the accuracy of the quadrature formula for all functions in : For any the inequality

holds, and this best possible. It is clear that is a function of the parameters , , of the quadrature formula and it is natural to try to choose them so that has a minimal value. This leads to a quadrature formula (from the class considered) whose error allows a minimal estimate for all functions of the space . Thus, the construction of a quadrature formula reduces to the solution of an extremal problem. This problem, even for the special case above, is extremely complicated and its solution has been obtained only for and .


[1] V.I. Krylov, "Approximate calculation of integrals" , Macmillan (1962) (Translated from Russian)
[2] S.M. Nikol'skii, "Quadrature formulas" , Hindushtan Publ. Comp. , London (1974) (Translated from Russian)
[3] V.I. Krylov, L.T. Shul'gina, "Handbook on numerical integration" , Moscow (1966) (In Russian)



[a1] P.J. Davis, P. Rabinowitz, "Methods of numerical integration" , Acad. Press (1975)
How to Cite This Entry:
Quadrature formula. Encyclopedia of Mathematics. URL:
This article was adapted from an original article by I.P. Mysovskikh (originator), which appeared in Encyclopedia of Mathematics - ISBN 1402006098. See original article