New Derivative Based Open Newton-Cotes Quadrature Rules

and Applied Analysis 3 2.1.1. Open Newton-Cotes Rules with Derivative at All Interior Nodes (ONC1). We can generate a scheme with a higher precision than that of the open Newton-Cotes scheme by using the first derivative of the integrand at all nodes, except at the endpoints, within the numerical integration scheme. We assume that the interval [a, b] is subdivided into n + 2 subintervals out of which n + 1 points are included in this problem. The procedure is given as


Introduction
In the recent past, the numerical quadrature rules were widely addressed in terms of their higher order of accuracy.The general form of any numerical integration formula is If the weights   are constant and the +1 nodes  0 ,  1 , . . .,   are equispaced within the interval [, ], the formula (1) gives rise to the Newton-Cotes quadrature rules.In open Newton-Cotes formulae, function evaluation at the endpoints of the interval is excluded from the quadrature rule.That is, where  0 ,  1 , . . .,   are distinct  + 1 integration points and   are  + 1 weights within the interval (, ) with   =  + ( + 1)ℎ,  = 0, 1, 2, . . ., , and ℎ = ( − )/( + 2).
The classical open rules are given as follows.
The Newton-Cotes rules can be constructed by replacing the integrand with an interpolating polynomial of appropriate degree.Lagrangian interpolating polynomials are extensively used to construct the Newton-Cotes rules; however, it was noted later that if, instead of Lagrangian polynomials, Hermite's interpolating polynomial is used which interpolates not only the function but also the first derivatives of the function at the selected nodes, a new class of integration rules are introduced called as the corrected Newton-Cotes formulas which include the first order derivative of the function as well.
The corrected Newton-Cotes formulas have greater degree of precision and are more accurate than the classical ones.For example, the Trapezoidal rule is improved by the corrected Trapezoidal rule with greater accuracy.Markoff also noticed that the Gaussian formulas can also be constructed by integrating a Hermitian interpolating polynomial.Numerous new quadratures were found which satisfy optimality properties for various classes of functions [1].These numerical quadratures have gained a lot of importance in recent years.A lot of work has been done in the numerical improvement of Newton-Cotes formulas.In 2003, Moawwad gave a unified approach for Newton-Cotes formulas of open as well as closed type by solving directly linear systems with coefficients matrices of Vandermonde type [2].In 2005, Sermutlu presented a detailed comparison of Newton-Cotes and Gaussian methods of quadrature [3].In 2005 and 2006, Dehghan and his companions [4][5][6] presented some new numerical integration formulas as an improvement of open, semiopen, and closed Newton-Cotes quadrature formulas.The new formulas by Dehghan et al. have an increased order of accuracy as compared to classical rules.The locations of boundaries along with weights were included as two additional parameters.The concept of degree of precision was used to set a system of nonlinear equations to obtain the approximate values of the parameters involved by solving the system approximately.In 2006, these authors introduced improvements of first and second kind Chebyshev-Newton-Cotes quadrature rules [7,8].In 2012, Burg introduced a new family of closed Newton-Cotes numerical integration formulas using first derivative values as well as functional values using the concept of degree of precision [9] to include additional parameters and thus to obtain new closed Newton-Cotes rules with increased order of accuracy.Recently, Burg and Degny [10] have also presented a new corrected family of midpoint rule involving odd order derivatives at the endpoints of the interval of integration.Also, in [11] Zhao and Li presented some new classes of midpoint derivatives based closed Newton-Cotes rules.
We, in this paper, construct a family of open Newton-Cotes formulas involving the first order derivatives at some or all nodes.The weights of the first derivative evaluation serve as additional parameters to increase the order of accuracy of these rules.The weights are determined by using the concept of degree of precision from a system of linear equations.The new formulas have increased order of accuracy for less number of nodes and less computational cost in comparison with the classical open Newton-Cotes rules.

Methodology and Derivation of Derivative
Based Open Quadrature Formulas (1) Function value and first derivative at each point   ∈ (, ).
The general form of the scheme is where   and   are the weights for the function and the derivative, respectively, and  is the number of subintervals used in the basic formula.The general form of the scheme is where   and   are as defined above.(3) Function value at each point   ∈ (, ) and first derivative at interior endpoints  0 ,   .The general form of the scheme is where   and   are as defined above.(4) Function value at each point   ∈ (, ) and first derivative at  −1 = ,  +1 = .
The general form of the scheme is where   and   are as defined above.
We, now, derive the formulas for these cases.

Open Newton-Cotes Rules with Derivative at All Interior Nodes (ONC1).
We can generate a scheme with a higher precision than that of the open Newton-Cotes scheme by using the first derivative of the integrand at all nodes, except at the endpoints, within the numerical integration scheme.We assume that the interval [, ] is subdivided into  + 2 subintervals out of which  + 1 points are included in this problem.The procedure is given as where   and   are the weights for the function and the derivative, respectively, and  is the number of subintervals used in the basic formula.We assume that degree of precision of ( 8) can be at most 2 + 1.This imparts the condition on (8) that it will be exact for all monomials   from  = 0 to  = 2 + 1.From the concept, we form a system of 2 + 2 equations.
For example, for  = 0, the two equations which are generated by this approach are For  () = 1 : ∫ Solving these equations for the coefficients  0 , and  0 , the following numerical quadrature rule is obtained: for  ∈ ( −1 ,  1 ) which is the classical midpoint rule.The precision of this method is 1 and it is 3rd order accurate.
For  = 1, the four equations generated are Solving these equations simultaneously for the unknown coefficients  0 ,  1 ,  0 , and  1 , the numerical quadrature rule becomes for  ∈ ( −1 ,  2 ).The precision of this method is 3 and it is 5th order accurate and the associated composite method is of the 4th order.For this method in its basic form two function evaluations and two first derivative evaluations are required.In composite form, where  is an integer divisible by 3, the method requires 2/3 function evaluations and 2/3 derivative evaluations.
For  = 2, the numerical quadrature rule is for  ∈ ( −1 ,  3 ).The precision of this method is 5 with order local truncation error of the 7th order.The associated composite rule is 6th order accurate.For this method in its basic form three function evaluations and two first derivative evaluations are required.In composite form, where  is an integer divisible by 4, the method requires 3/4 function evaluations and /2 derivative evaluations.
For  = 3, the numerical quadrature rule is for  ∈ ( −1 ,  4 ).The precision of this method is 7 with local truncation error of the 9th order.The associated composite rule is 8th order accurate.For this method in its basic form four function evaluations and two first derivative evaluations are required.In composite form, where  is an integer divisible by 5, the method requires 4/5 function evaluations and 4/5 derivative evaluations.
For  = 4, the numerical quadrature rule is for  ∈ ( −1 ,  5 ).The precision of this method is 9 with local truncation error of the 11th order and the associated composite rule is 10th order accurate.For this method in its basic form five function evaluations and two first derivative evaluations are required.In composite form, where  is an integer divisible by 6, the method requires 5/6 function evaluations and 2/3 derivative evaluations.

Open Newton-Cotes Rules with Derivative at All Nodes including Endpoints (ONC2).
The procedure is given as where   and   are the weights for the function and the derivative, respectively, and  is the number of subintervals used in the basic formula.We assume that degree of precision of ( 9) can be at most 2 + 3.This imparts the condition on (9) that it will be exact for all monomials   from  = 0 to  = 2 + 3. From the concept, we form a system of 2 + 4 equations.
For example, for  = 0, the following numerical quadrature rule is obtained: for  ∈ ( −1 ,  1 ).The precision of this method is 3 and it is 5th order accurate.The associated composite rule is 4th order accurate.For this method in its basic form one function evaluation and two first derivative evaluations are required.
In composite form, the method requires /2 function evaluations and only two additional derivative evaluations due to cancellation of the derivatives at the endpoints in contrast to one point classical open method.
For  = 1, the quadrature rule becomes for  ∈ ( −1 ,  2 ).The precision of this method is 5 and it is 7th order accurate.The associated composite rule is 6th order accurate.For this method in its basic form two function evaluations and four first derivative evaluations are required.
In composite form, where  is an integer divisible by 3, the method requires 2/3 function evaluations and 2/3 + 2 derivative evaluations due to cancellation of the derivatives at the endpoints.
For  = 2, the numerical quadrature rule is for  ∈ ( −1 ,  3 ).The precision of this method is 7 with local truncation error of the 9th order.The associated composite rule is 8th order accurate.For this method in its basic form three function evaluations and four first derivative evaluations are required.In composite form, where  is an integer divisible by 4, the method requires 3/4 function evaluations and /2 + 2 derivative evaluations due to cancellation of the derivatives at the endpoints.
For  = 3, the numerical quadrature rule is for  ∈ ( −1 ,  4 ).The precision of this method is 9 with local truncation error of the 11th order.The associated composite rule is 10th order accurate.For this method in its basic form four function evaluations and six first derivative evaluations are required.In composite form, where  is an integer divisible by 5, the method requires 4/5 function evaluations and 4/5 + 2 derivative evaluations due to cancellation of the derivatives at the endpoints.
For  = 4, the numerical quadrature rule is for  ∈ ( −1 ,  5 ).The precision of this method is 11 with local truncation error of the 13th order and the associated composite rule is 12th order accurate.For this method in its basic form five function evaluations and six first derivative evaluations are required.In composite form, where  is an integer divisible by 6, the method requires 5/6 function evaluations and 2/3 + 2 derivative evaluations due to cancellation of the derivatives at the endpoints of the basic form.

Open Newton-Cotes Rules with Derivative at Interior Endpoints (ONC3).
In this approach, the evaluation of the derivative is restricted for the two points only, while maintaining the improvement of order of accuracy.These numerical integration formulas follow the same behaviour as the Newton-Cotes formulas with respect to  being even or odd.These numerical integration formulas are two order of magnitude greater than the corresponding open Newton-Cotes formulas.
The general form of the scheme is In these composite forms, all of the interior derivative calculations cancel so that the only 2 derivative calculations are required for the composite form.
For  = 0, the quadrature formula is given as which is the midpoint rule.The four equations are shown below for  = 1: Solving these equations for unknown parameters we get This scheme has precision 3 and 5th order local truncation error and the associated composite rule is 4th order accurate which is the same as (16).
For  = 2, the optimal quadrature formula is This scheme has precision 5 and 7th order local truncation error and the associated composite rule is 6th order accurate.For this method in its basic form three function evaluations and two first derivative evaluations are required.In composite form, where  is an integer divisible by 4, the method requires 3/4 function evaluations and /2 derivative evaluations. 6 Abstract and Applied Analysis For  = 3, the optimal quadrature formula is This scheme has precision 7 and 9th order local truncation error and the associated composite rule is 8th order accurate.For this method in its basic form five function evaluations and two first derivative evaluations are required.In composite form, where  is an integer divisible by 6, the method requires 5/6 function evaluations and /3 derivative evaluations.

Open Newton-Cotes Rules with Derivative at Endpoints (ONC4).
The general form of the scheme is For  = 0, the quadrature rule is for  ∈ ( −1 ,  1 ) which is the same as (21).
For  = 1, the four equations are given below: By solving these equations, we get the following rule: for  ∈ ( −1 ,  2 ).This scheme has precision 3 and 5th order local truncation error and the associated composite rule is 4th order accurate.
For  = 2, the quadrature rule is for  ∈ ( −1 ,  3 ).This scheme has precision 5 and 7th order local truncation error and the associated composite rule is 6th order accurate.
For  = 3, the quadrature rule is for  ∈ ( −1 ,  4 ).This scheme has precision 5 and 7th order local truncation error and the associated composite rule is 6th order accurate.
For  = 4, the quadrature rule is for  ∈ ( −1 ,  5 ).This scheme has precision 7 and 9th order local truncation error and the associated composite rule is 8th order accurate.It may be observed that all these methods require only two additional first derivative evaluations in their basic as well as in the composite form due to cancellation of the derivatives at the endpoints in contrast to the classical open rules of the respective domain.

Error Analysis
Error Analysis Approach 1.To find error term, several methods can be used.Using the precision, the error term is related to the difference between the quadrature formula for the polynomial  +1 /( + 1)! and the exact result where  is the precision of the quadrature formula.This difference yields the coefficient for the ( + 1)th derivative of ().
Error Analysis Approach 2. The third approach is based on the interpolating polynomials.Since it is difficult to justify why to ignore the higher order monomials in the above mentioned approaches.The error analysis is carried out by using the classical approach of examining the integral of the remainder term of the Lagrange interpolating polynomial.Thus, we conjecture the following theorem for the error analysis of open derivative based quadrature by following an approach similar to [9].
Theorem 1 (see [9]).Suppose a function () ∈  (+) This theorem can be used to calculate the error terms of our scheme as follows.
We, now, show that these higher order derivatives based quadrature formulas are computationally efficient.The numerical results are obtained in Maple 13 using a precision of twenty decimal digits.

Computational Efficiency of Derivative Based Open Rules
Now, we will use the following notations in the comparison tables.
ONC for classical open numerical quadrature rule.
ONC1 for open numerical quadrature with derivative at all interior points.
ONC2 for open numerical quadrature with derivative at all points including endpoints.
ONC3 for open numerical quadrature with derivative used only at interior endpoints.
ONC4 for open numerical quadrature with derivative used only at exterior endpoints.

Table 1 compares the coefficient of the classical open Newton-Cotes formulas with newly developed open Newton-
Cotes formulas with derivatives.The order described in the first column is the order of the composite formula;  0 represents the number of nodes.The columns below each heading of the open type quadrature rules represent the coefficients of the error terms.It may be observed from the table that the new method ONC2 has prominently smaller error term coefficients.
Tables 2, 3, 4, 5 and 6 show the number of function and derivative evaluations for the classical open Newton-Cotes rules and the new derivative based open Newton-Cotes formulas, as well as their orders of accuracy.These results will be used in Tables 7 and 8 to determine the number of function and derivative evaluations sufficient to reduce the error in the calculations below some threshold.The results in Table 4 show that ONC2 is a highly attractive option because it has the highest order of accuracy in contrast to the other methods and is twice the value of N. Additionally, as shown in Table 1, the coefficients for the leading order error term for ONC2 decrease more rapidly than the other formulas.The results in Table 6 show that ONC4 is also a better option due to the increase in the order of accuracy by a factor of two over the classical open Newton-Cotes for mulas but includes only two additional derivative evaluations, regardless of the number of intervals for a certain scheme.Because ONC2 and ONC4 have clear advantages to ONC1 and ONC3, we will focus our analysis on these two methods.
We, now, give the following examples to compare the computational costs of the new rules and the classical rules.Example 2.
Tables 7 and 8 show the total number of evaluations required to guarantee that the numerical results for these integrals have an error of less than 10 −9 .These results are obtained from the error terms calculated previously.Here,  represents the number of intervals involved in obtaining the required accuracy.
Assuming that the computational cost of evaluating the derivative of the function is the same as the cost of evaluating the function, the fourth order method in terms of computational cost is ONC2 ( = 0) which is the same formula as ONC4 ( = 0), which has roughly half the cost of the original fourth order open Newton-Cotes method.The best sixth order accurate method is ONC4 ( = 2), which is roughly 60% of the cost of the sixth order ONC method.
Tables 9 and 10 show the average CPU usage time required to process 100 simulations for the two above mentioned integrals using the classical open and new derivative based open methods ONC2 and ONC4.The time is measured in seconds.

Numerical Results
We, in this section, approximate the values of the following examples using the newly derived open derivative based formulas.These integrals contain rational, exponential, trigonometric, and logarithmic functions as integrands.
where (ℎ) means the numerical result with the step size ℎ.(0) is the exact value of the definite integrals.Thus, (0) is 0.162865005 in the first example, 0.882081390 in the second, and 0.199714662 in the third example.The comparison is done with the classical open Newton-Cotes rules of the respective order of accuracy.The first row indicates the partitioning of the interval into subintervals. represents the observed order of accuracy and the results in line with each heading of Open Newton-Cotes rules represent the approximated value of the definite integral calculated through that formula.

Conclusion
We have used several approaches to construct the derivative based open Newton-Cotes formulas.These new schemes represent four different families of new derivative based open Newton-Cotes type quadrature formulas.Using the concept of precision, we determine the coefficients for the function and its derivative at different locations within and at the boundary of the integral and determine the error term for each method from  = 0 through  = 4 where  is the number of interior points.Because of the beneficial cancellation of the derivatives at the boundaries of the intervals, the most attractive methods are ONC2 (open Newton-Cotes with derivatives at all nodes including endpoints) and ONC4 (open Newton-Cotes with derivatives only at the end points).When comparing the computational cost of these two methods with that of the same order open Newton-Cotes methods (ONC), the ONC2 methods are consistently more efficient than the original ONC methods, and the ONC4 methods are the most computationally efficient in producing a numerical result that is guaranteed to be accurate to within 10 −9 .These conclusions are obtained from both the theoretical order of accuracy and error terms as well as from computational evidence.
Finally, the order of accuracy for all of the methods is confirmed computationally against several different integrals, showing that these methods are either comparable to or better than the classical open Newton-Cotes formula.
2.1.Derivative Based Schemes.We introduce the following derivative based schemes for the open Newton-Cotes formulas.

Table 1 :
Coefficients of error term for each order of accuracy.Let () be a  + 2 degree polynomial that agrees with () and the first derivative of () at 2 locations  −1 =  and  +1 = .Then, there exists a location  ∈ (, ) such that

Table 2 :
Computational cost for ONC (FE-Number of Function Evaluations, DE-Number of Derivative Evaluations).

Table 7 :
Computational cost for evaluating ∫

Table 9 :
Average CPU usage time for evaluating ∫

Table 11 :
Results and observed orders for the rules with order 4, for ∫

Table 12 :
Results and observed orders for the rules with order 4, for ∫

Table 13 :
Results and observed orders for the rules with order 4, for ∫

Table 14 :
Results and observed orders for the rules with order 6, for ∫

Table 15 :
Results and observed orders for the rules with order 6, for ∫

Table 16 :
Results and observed orders for the rules with order 6, for ∫

Table 17 :
Results and observed orders for the rules with order 8, for ∫

Table 18 :
Results and observed orders for the rules with order 8, for ∫

Table 19 :
Results and observed orders for the rules with order 8, for ∫

Table 20 :
Results and observed orders for the rules with order 10, for ∫

Table 21 :
Results and observed orders for the rules with order 10, for ∫ 2 0  − 2 .

Table 22 :
Results and observed orders for the rules with order 10, for ∫We evaluate Examples 1-3 in Tables 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, and 22 with their observed order of accuracy.The order of accuracy is calculated by the formula given as