Quadrature

The word quadrature has (at least) three incompatible meanings. Integration by quadrature either means solving an integral analytically (i.e., symbolically in terms of known functions), or solving of an integral numerically (e.g., Gaussian quadrature, Newton-Cotes formulas). Ueberhuber (1997, p. 71) uses the word "quadrature" to mean numerical computation of a univariate integral, and "cubature" to mean numerical computation of a multiple integral.

The word quadrature is also used to mean squaring: the construction of a square using only compass and straightedge which has the same area as a given geometric figure. If quadrature is possible for a plane figure, it is said to be quadrable.

For a function f(x) tabulated at given values x_i (so the abscissas cannot be chosen at will), write the function phi as a sum of orthonormal functions p_j satisfying

 int_a^bp_i(x)p_j(x)W(x)dx=delta_(ij)
(1)

as

 phi(x)=sum_(j=0)^inftya_jp_j(x),
(2)

and plug into the Lagrange interpolating polynomial of f(x) through the m points (as is done in Gaussian quadrature)

int_a^bphi(x)W(x)dx=int_a^bsum_(j=1)^(m)(pi(x)W(x))/((x-x_j)pi^'(x_j))dxf(x_j)
(3)
=sum_(j=1)^(m)w_jf(x_j),
(4)

where

 pi(x)=product_(j=1)^m(x-x_j),
(5)

giving

 int_a^bsum_(j=0)^inftya_jp_j(x)W(x)dx=sum_(i=1)^nw_i[sum_(j=0)^inftya_jp_j(x_i)].
(6)

But we wish this to hold for all degrees of approximation, so

 a_jint_a^bp_j(x)W(x)dx=a_jsum_(i=1)^nw_ip_j(x_i)
(7)
 int_a^bp_j(x)W(x)dx=sum_(i=1)^nw_ip_j(x_i).
(8)

Setting i=0 in (◇) gives

 int_a^bp_0(x)p_j(x)W(x)dx=delta_(0j).
(9)

The zeroth order orthonormal function can always be taken as p_0(x)=1, so (9) becomes

int_a^bp_j(x)W(x)dx=delta_(0j)
(10)
=sum_(i=1)^(n)w_ip_j(x_i),
(11)

where (◇) has been used in the last step. We therefore have the matrix equation

 [p_0(x_1) ... p_0(x_n); p_1(x_1) ... p_1(x_n); | ... |; p_(n-1)(x_1) ... p_(n-1)(x_n)][w_1; w_2; |; w_n]=[1; 0; |; 0],
(12)

which can be inverted to solve for the w_is (Press et al. 1992).

Wolfram Web Resources

Mathematica »

The #1 tool for creating Demonstrations and anything technical.

Wolfram|Alpha »

Explore anything with the first computational knowledge engine.

Wolfram Demonstrations Project »

Explore thousands of free applications across science, mathematics, engineering, technology, business, art, finance, social sciences, and more.

Computerbasedmath.org »

Join the initiative for modernizing math education.

Online Integral Calculator »

Solve integrals with Wolfram|Alpha.

Step-by-step Solutions »

Walk through homework problems step-by-step from beginning to end. Hints help you try the next step on your own.

Wolfram Problem Generator »

Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.

Wolfram Education Portal »

Collection of teaching and learning tools built by Wolfram education experts: dynamic textbook, lesson plans, widgets, interactive Demonstrations, and more.

Wolfram Language »

Knowledge-based programming for everyone.