Abstract

The derivative-based trapezoid rule for the Riemann-Stieltjes integral is presented which uses 2 derivative values at the endpoints. This kind of quadrature rule obtains an increase of two orders of precision over the trapezoid rule for the Riemann-Stieltjes integral and the error term is investigated. At last, the rationality of the generalization of derivative-based trapezoid rule for Riemann-Stieltjes integral is demonstrated.

1. Introduction

In mathematics, the Riemann-Stieltjes integral is a kind of generalization of the Riemann integral, named after Bernhard Riemann and Thomas Stieltjes. It is Stieltjes that first gave the definition of this integral [1] in 1894. It serves as an instructive and useful precursor of the Lebesgue integral. It is known that the Riemann-Stieltjes integral has wide applications in the field of probability theory [2, 3], stochastic process [4], and functional analysis [5], especially in original formulation of F. Riesz’s theorem and the spectral theorem for self-adjoint operators in a Hilbert space.

Definite integration is one of the most important and basic concepts in mathematics. And it has numerous applications in fields such as physics and engineering. In several practical problems, we need to calculate integrals. As is known to all, as for , once the primitive function of integrand is known, the definite integral of over the interval is given by Newton-Leibniz formula, that is,

However, the explicit primitive function is not available or its primitive function is not easy to obtain, such as , , and . Moreover, the integrand is only available at certain points , . How to get high-precision numerical integration formulas becomes one of the challenges in fields of mathematics [6].

The methods of quadrature are usually based on the interpolation polynomials and can be written in the following form: where there are distinct integration points at within the interval and weights . If the integration points are uniformly distributed over the interval, so in which .

These can be derived in several different ways [79]. One is interpolate at the points , using the Lagrange polynomials and then integrating the foresaid polynomials to obtain (2).

The other is based on the precision of a quadrature formula. Select the , , so that the error, is exactly zero for , . Using the method of undetermined coefficients, this approach generates a system of linear equations for weights . Since the monomials are linearly independent, the linear system of equations has a unique solution.

The trapezoidal rule is the most well-known numerical integration rule of this type. Trapezoidal rule for classical Riemann integral is where .

In spite of the many accurate and efficient methods for numerical integration being available in [79], recently Mercer has obtained trapezoid rule for Riemann-Stieltjes integral which engender a generalization of Hadamard’s integral inequality [10]. Then he develops Midpoint and Simpson’s rules for Riemann-Stieltjes integral [11] by using the concept of relative convexity. Burg has proposed derivative-based closed Newton-Cotes numerical quadrature [12] which uses both the function value and the derivative value on uniformly spaced intervals. Zhao and Li have proposed midpoint derivative-based closed Newton-Cotes quadrature [13] and numerical superiority has been shown. Recently, Simos and his partners have made a contribution to the Newton-Cotes formula for the Riemann integral and its applications [1421], especially the connection between closed Newton-Cotes, trigonometrically fitted differential methods, symplectic integrators, and efficient solution of the Schrodinger equation [1721].

Motivation for the research presented here lies in construction of derivative-based trapezoid rule for the Riemann-Stieltjes integral, which is generalization of the results in [1012]. In this paper, the derivative-based trapezoid rule for the Riemann-Stieltjes integral is presented. This new scheme is investigated in Section 2. In Section 3, the error term is presented. Finally, conclusions are drawn in Section 4.

2. Derivative-Based Trapezoid Rule for the Riemann-Stieltjes Integral

In this section, by adding the derivatives at the endpoints, derivative-based trapezoid rule for the Riemann-Stieltjes integral is presented.

Theorem 1. Suppose that and are continuous on and is increasing there. The derivative-based trapezoid rule for the Riemann-Stieltjes integral is

Proof. Looking for the derivative-based trapezoid rule for the Riemann-Stieltjes integral, we seek numbers , , , such that is equality for . That is,
Therefore,
Solving simultaneous (8) for , we obtain
So we have the derivative-based trapezoid rule for the Riemann-Stieltjes integral as desired.

Corollary 2. The precision of the derivative-based trapezoid rule for the Riemann-Stieltjes integral is 3.

Proof. From the construction of , , , , we obtain that the derivative-based trapezoidal rule for the Riemann-Stieltjes integral has degree of precision not less than 3.
In Section 3 Theorem 3, we can easily see that the quadrature is not equality for . So the precision of this method is 3.

3. The Error Term for Riemann-Stieltjes Derivative-Based Trapezoid Rule

In this section, the error term of the derivative-based trapezoid rule for the Riemann-Stieltjes is investigated. The error term can be found in mainly 3 different ways [8, 9].

Here, we use the concept of precision to calculate the error term, where the error term is related to the difference between the quadrature formula for the monomial and the exact value , where is the precision of the quadrature formula.

Theorem 3. Suppose that and are continuous on and is increasing there. The derivative-based trapezoid rule for the Riemann-Stieltjes integral with the error term is where , . And the error term of this method is

Proof. Let .
So
By Theorem 1, we have By Equations (12)-(13), we obtain This implies that

Corollary 4. Conditions are the same as for Theorem 3. When , (10) reduces to the derivative-based trapezoid rule (see [12]) for the classical Riemann integral.

Proof. It is easy to obtain
By Theorem 3,

Remark 5. From Corollary 4, we know that the results in Theorem 3 possess the reducibility. When , formula (10) reduces to the derivative-based trapezoid rule for the classical Riemann integral. So Theorem 3 is a reasonable generalization of the results in [12].

4. Conclusions

We briefly summarize our main conclusions in this paper as follows. (1)The derivative-based trapezoid rule for the Riemann-Stieltjes integral is presented which uses 2 derivative values at the endpoints. (2)This kind of quadrature rule obtains an increase of two orders of precision over the trapezoid rule for the Riemann-Stieltjes integral. (3)The error term for Riemann-Stieltjes derivative-based trapezoid rule is investigated. And the rationality of the generalization of derivative-based trapezoid rule for Riemann-Stieltjes integral is demonstrated.

The derivative-based midpoint and Simpson’s rules for the Riemann-Stieltjes integral will be achieved by further research.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the Scientific Research Foundation of Civil Aviation University of China (no. 2013QD01X), the Fundamental Research Funds for the Central Universities (no. 3122014C023), and the National Natural Science Foundation of China (no. 71171190).