Abstract

This study examined the characteristics of a variable three-point Gauss quadrature using a variable set of weighting factors and corresponding optimal sampling points. The major findings were as follows. The one-point, two-point, and three-point Gauss quadratures that adopt the Legendre sampling points and the well-known Simpson’s 1/3 rule were found to be special cases of the variable three-point Gauss quadrature. In addition, the three-point Gauss quadrature may have out-of-domain sampling points beyond the domain end points. By applying the quadratically extrapolated integrals and nonlinearity index, the accuracy of the integration could be increased significantly for evenly acquired data, which is popular with modern sophisticated digital data acquisition systems, without using higher-order extrapolation polynomials.

1. Introduction

Numerical integration methods may be grouped in two categories. One is the rule for discrete data and the other is for function of continuous data. The Gauss-Legendre quadrature [16] is a well-known rule belonging to the latter category. In the former category for discrete data, the Newton-Cotes [79] method is applicable with many orders of integration. The 1st-order rule is the same as the trapezoidal rule [7] which can be extended to the Romberg integration method [7]. The 2nd-order rule is the same as Simpson’s 1/3rd rule, and 3rd-order rule is the same as Simpson’s 3/8th rule [7]. The 4th-order rule of Newton-Cotes is known as Boole’s rule using 5 data points [7]. As the order of integration in the Newton-Cotes methods increases, Runge’s phenomenon arises and the accuracy of integration becomes worse owing to the fluctuation of the interpolated higher-order polynomials [10].

The Gauss-Legendre quadrature uses function values at interior sampling points with corresponding best weights to result in a very accurate result in spite of the relatively small number of sampling points. However, this quadrature is inapplicable to discrete data points because it does not use boundary point data. In this study, we provide a lemma with a formula for the new 3-point Gauss quadrature of variable sampling points which include the Legendre point as well.

An examination of the effect of these varied sampling points found that the one-point, two-point, and three-point Gauss quadratures adopting the Legendre sampling points and Simpson’s 1/3 rule were actually special cases of the variable three-point Gauss quadrature. The order of the polynomial that can be integrated precisely by conventional three-point Gauss quadrature is 5. On the other hand, the use of varied sampling points with a variable three-point Gauss quadrature does not allow exact integration of a polynomial order 5.

Despite their reduced accuracy, variable Gauss quadratures can be applied effectively to special situations, for example, the shear-locking problems that arise when using the finite element method for a plate/shell when the ratio of the thickness to the width is quite small, as reported previously [11, 12]. Another effective application is integration of discrete data, which is the same as Simpson’s rule, rather than the conventional Gauss quadrature which is inapplicable owing to its use of inner-range data. For a weighting factor of at the central sampling point, the optimal sampling positions for the outer sampling points were obtained as end points. Therefore, when using this specific group of weighting factor and sampling points, the end point Gauss quadrature becomes equal to Simpson’s 1/3 rule. If the central weighting factor of is increased further to or , this produces a 1st extended end point integration or 2nd extended end point integration, where the outer sampling points are located at the extended outer end points or .

Accordingly, this study examined the characteristics of various groups of weighting factors and sampling points and tested the performance of the extended end points quadrature using the outer out-of-domain sampling points. A new method, adopting quadratically extrapolated integrals and a nonlinearity index, using integrals of the variable three-point Gauss integrations of the 1st and 2nd extended end points and conventional end point integration, was applied to the integration of evenly acquired discrete data to obtain new four kinds of numerical integration formulae.

2. Variable Three-Point Gauss Integration Formula with Variable Weighting Factors and Sampling Points

The modification of the Gauss integration formula with a near-zero center-weight factor was included in a previous study [11]. This paper presents a comprehensive interpretation of the variable three-point Gauss integration formula with effective applications. Essentially, to vary the Gauss integration weights, the weighting factors are represented by a central weighting factor , and any changes in the central weighting factor are reflected in the weights and locations of outer sampling points for the optimal integration. The detailed procedure to derive the variable three-point Gauss integration formula with the weighting factors is as follows [11]: where , , and are the sampling points within the range of . In addition, the integration weights (, , and ) are given by the following expressions: where is the Legendre polynomial formula,

The general expressions for the sampling points and weights of the Gauss integration can be driven as follows. In (2), if it is assumed that the weight at and the other weights are , at and , respectively, then the optimum results can be obtained as follows:

From (4), the relationship between the central weighting factor and the corresponding optimal location of outer sampling point is as follows:

We present a lemma regarding three-point Gauss quadrature of variable sampling point.

Lemma 1. Integral from −1 to 1 of function is obtained numerically by the equation where which is exact for the polynomial of order up to 3.

Proof. Without loss of generality, we assumed the range of integration , which is a mapped variable, from −1 to 1. For the odd monomials with odd index such as and , the integrals are 0 obviously. For the constant integrand , the integral is exactly , because sum of the weights is 2. We only need to verify that the integral of the integrand is .
Equation (6) for the integrand of is

Equation (7) shows the general expression for the sampling points and weights of the variable three-point Gauss integration rule. The weighting factor (=) can be selected arbitrarily within the range of 0~2.

3. Classification of a Set of Weighting Factors and Sampling Points and Their Integration Characteristics

In the variable three-point Gauss integration formula, the term of “variable” was adopted because the weighting factor can be selected arbitrarily within the range of 0~2. In addition, any change in the central weighting factor is also reflected in the weighting factors and locations of outer sampling points, thereby changing the numerical integration characteristics ((5) and (7)). These different groups of integration weights and sampling points can be classified according to their weight ratio .

Table 1 lists nine fundamental weight ratios. The type (a) ratio of corresponds to the conventional two-point Gauss-Legendre quadrature, which is a special case of the variable three-point Gauss quadrature. The type (b) ratio of corresponds to the quasi two-point rule, which is applied in FEM to prevent numerical instability [11]. The type (c) ratio of corresponds to the even weight three-point rule that might be applied effectively to some problems in the future. The type (d) ratio of corresponds to the conventional three-point Gauss-Legendre quadrature. The type (e) ratio of corresponds to the double-center-weight three-point rule that might also be applied effectively to some problems in the future. The type (f) ratio of corresponds to the end point rule, which is the same as Simpson’s 1/3 rule, where the outer sampling points are located at the end points of a dimensionless domain of . The type (g) ratio of corresponds to the 1st extended end point rule, where the outer sampling points are located at the end points of an extended dimensionless domain of . The type (h) ratio of corresponds to the 2nd extended end point rule, where the outer sampling points are located at the end points of a doubly extended dimensionless domain of . Finally, the type (i) ratio of corresponds to the extended end point rule, where the outer sampling points are located at the end points of an infinitely extended dimensionless domain of , which is actually the conventional one-point Gauss-Legendre quadrature.

The variable three-point rule of a Gauss quadrature includes the conventional Gauss-Legendre quadrature of the one-point rule and two-point rule , as well as the three-point rule . In addition, it also includes the well-known Simpson’s 1/3 rule when . This can also be expanded further to the cases of ,  ,  ,  , and , respectively. Each case has particular meaning related to the weighting factors or sampling points. Figure 1 shows the proportional change in the integration weights and sampling points according to the weight ratio.

We present three corollaries regarding three-point Gauss quadrature of variable sampling point.

Corollary 2. The three-point Gauss quadrature of Lemma 1 (6)~(7) reduces to one-point Gauss-Legendre quadrature with .

Corollary 3. The three-point Gauss quadrature of Lemma 1 (6)~(7) reduces to two-point Gauss-Legendre quadrature with .

Corollary 4. The three-point Gauss quadrature of Lemma 1 (6)~(7) reduces to Simpson’s 1/3 rule with .

4. Extrapolated Quadratic Polynomial of the Integrals Using Data from End Point Integral and Extended End Point Integrals for Evenly Acquired Discrete Data

To test the characteristics of variable three-point Gauss integration formulae, the error ratios of the integral from 0 to 2 for the 5 types of monomial integrand were compared (Table 2). Most types exhibited accurate results up to the cubic monomial except in the case of a weight ratio (conventional one-point rule). The type of weight ratio (conventional three-point rule of Legendre sampling points) results in very accurate result up to the monomial, as expected.

According to the formula of Gauss quadrature of variable three sampling points with corresponding best weights, one can see that the sampling points may be located at the boundary points (integral ) which means that the method is applicable to discrete data points. Moreover, the outer sampling points may be located out of region such as 2nd next discrete points (integral ) and 3rd next discrete points (integral ). Using the integrals , , and , we can extrapolate the integral function , like Richardson’s extrapolation to interval [7], to the case of Legendre sampling point of , where there is no data.

In other words, when one integrates a function from −1~1 in local coordinate , the best sampling points are the Legendre points of , ., and . However, data are only on the discrete points of Thus, we compute the following three integrals first.

Then, we extrapolate these integrals to compute the integral , which cannot be obtained directly as follows because there is no data on the sampling point :

The process of extrapolation is presented here. The integral, which is a function of outer sampling point , may be treated as a quadratic polynomial and its coefficients are determined through the interpolating points , , and . , , and are given in (9)~(11).

On the other hand, the error of the integral choosing the outer sampling point as , 2, or 3 increases explosively. Therefore, a nonlinearity index, , is used in a quadratically interpolated integral, which is a function of the absolute position of the outer sampling point . Hence, the quadratic polynomial of the integral, depending on the position of the outer sampling point of the three-point Gauss quadrature, can be represented as follows: where nonlinearly mapped variable of the absolute position of the outer sampling point , : end point integration ,  : 1st extended end point integration , : 2nd extended end point integration , and nonlinearity index.

After obtaining the coefficients, , , and , of the polynomial, the best expected integral can be predicted with the nonlinearly mapped extrapolated outer sampling point , where the point is the Legendre sampling point of the conventional three points of Gauss quadrature, and that point originally could not be sampled for this type of discrete data. This whole procedure can be repeated for the linear interpolation of the integral , where we use only two data of points and . and are given in (9)~(10).

We present four theorems regarding extrapolated Gauss quadrature of extended end points.

Theorem 5. Integral from −1 to 1 of discrete data is obtained numerically by formula I, which is exact for the polynomial of order up to 3 with the nonlinearity index (see Table 3).

Proof. As in the proof of Lemma 1, we only need to verify that the integral of the integrand is 2/3.
Equation (14) for the integrand is

Theorem 6. Integral from −1 to 1 of discrete data is obtained numerically by formula II, which is exact for the polynomial of order up to 5 with the nonlinearity index (see Table 4).

Proof. As in the proof of Theorem 5, we need to verify that the integral of the integrand is 2/3 and integral of is 2/5.
Equation (16) for the integrand is
Equation (16) for the integrand is

Theorem 7. Integral from −1 to 1 of discrete data is obtained numerically by formula III, which is exact for the polynomial of order up to 5 with the nonlinearity index (see Table 5).

Proof. As in the proof of Theorem 6, we need to verify that the integral of the integrand is 2/3 and integral of is 2/5.
Equation (19) for the integrand is
Equation (19) for the integrand is

Theorem 8. Integral from −1 to 1 of discrete data is obtained numerically by formula IV, which is exact for the polynomial of order up to 5 with the nonlinearity index (see Table 6).

As in the proof of Theorem 7, we need to verify that the integral of the integrand is 2/3 and integral of is 2/5.

Equation (22) for the integrand is

Equation (22) for the integrand is

Through this extrapolation to obtain the integral , we can integrate polynomial up to the order of 5 exactly and up to the order of 7 approximately, without any fluctuating phenomenon of Runge [10].

We compared the four formulae with Simpson’s 1/3 rule for integration of the discrete data from 1 to 3. The discrete data are given at integer point of from the monomials , , and exp(), respectively, as shown in Table 7.

Table 8 compares the four formulae with Boole’s rule for integration of the discrete data from 1 to 5. The discrete data are given at integer points of from the monomials , , and exp(), respectively. Formula I results in some error for the integrand of , because it is for the polynomial of order up to 3. Formula II and formula IV still give better result compared with Boole’s.

5. Conclusion

(1) A lemma for the variable three-point Gauss quadrature has been presented. Based on that, comprehensive sets of weighting factors and corresponding optimal sampling points were presented.

Three corollaries are presented to show that one-point and two-point Gauss-Legendre quadrature and Simpson’s 1/3 rule are special cases of the variable three-point Gauss quadrature.

The variable three-point Gauss quadrature can be applied to end point integration of discrete data, such as Simpson’s rule, and even to the 1st and 2nd extended end point integrations that use out of region discrete data.

Four theorems of extrapolation method with nonlinearity index ( or 2) are presented. The integrals of variable three-point Gauss integrations of the end point integration and extended end point integration are extrapolated to predict the value at dimensionless point.

Formula II and formula IV given in the theorems exhibit excellent accuracies compared with Simpson’s rule or Boole’s rule. The formulae can be applied effectively to evenly acquired data, which is the case in modern sophisticated digital data acquisition systems, without using higher-order extrapolation polynomials which might result in numerical instability.

Acknowledgment

This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF-2012R1A1A2008903).