## Noname manuscript No.

(will be inserted by the editor)

## The order of convergence of an optimal quadrature formula with derivative in the space W (2 1) , 2

Abdullo R. Hayotov , Rashidjon G. Rasulov ∗

Received: date / Accepted: date

Abstract The present work is devoted to extension of the trapezoidal rule in the space W (2 1) , 2 . The optimal quadrature formula is obtained by minimizing the error of the formula by coefficients at values of the first derivative of a integrand. Using the discrete analog of the operator d 2 dx 2 -1 the explicit formulas for the coefficients of the optimal quadrature formula are obtained. Furthermore, it is proved that the obtained quadrature formula is exact for any function of the set F = span 1 { , x, e x , e -x } . Finally, in the space W (2 1) , 2 the square of the norm of the error functional of the constructed quadrature formula is calculated. It is shown that the error of the obtained optimal quadrature formula is less than the error of the Euler-Maclaurin quadrature formula on the space L (2) 2 .

MSC: 65D30, 65D32.

Keywords: optimal quadrature formula, Hilbert space, the error functional, S.L. Sobolev's method, discrete argument function, the order of convergence.

## 1 Introduction

It is known, that quadrature and cubature formulas, are methods for the approximate evaluation of definite integrals. In addition and even more important, quadrature formulas provide a basic and important tool for the numerical solution of differential and integral equations. The theory of cubature formulas consists mainly of three branches dealing with exact formulas, formulas based on

A.R. Hayotov

∗ Corresponding author: Department of Mathematical Sciences, KAIST, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea

E-mail: hayotov@kaist.ac.kr.

V.I.Romanovskiy Institute of Mathematics, Uzbekistan Academy of Sciences, 81 M.Ulugbek str., Tashkent 100170, Uzbekistan,

E-mail: hayotov@mail.ru.

R.G. Rasulov

V.I.Romanovskiy Institute of Mathematics, Uzbekistan Academy of Sciences, 81 M.Ulugbek str., Tashkent 100170, Uzbekistan E-mail: r.rasulov1990@mail.ru.

functional-analytic methods, and formulas based on probabilistic methods [21,22]. In the functional-analytic methods the error between an integral and corresponding cubature sum is considered as a linear functional on a Banach space and it is estimated by the norm of the error functional in the conjugate Banach space. The norm of the error functional depends on coefficients and nodes of the formula. The problem of finding the minimum of the norm of the error functional by coefficients and by nodes is called S.M.Nikolskii problem , and the obtained formula is called the optimal formula in the sense of Nikolskii (see, for instance, [10]). Minimization of the norm of the error functional by coefficients when the nodes are fixed is called Sards problem . And the obtained formula is called the optimal formula in the sense of Sard . First this problem was studied by A. Sard [11]. Solving these problems in different spaces of differentiable functions various type of optimal formulas of numerical integration are obtained.

There are several methods for constructing the optimal quadrature formulas in the sense of Sard such as the spline method, the φ -function method (see e.g. [3], [13]) and the Sobolev method. It should be noted that the Sobolev method is based on using a discrete analog of a linear differential operator (see e.g. [20, 21,22]). In different spaces based on these methods, the Sard problem was studied by many authors, see, for example, [2,3,4,5,6,7,8,11,13,15,20,21,22,23] and references therein.

Among these formulas the Euler-Maclaurin type quadrature formulas are very important for numerical integration of differentiable functions and are widely used in applications. In different spaces the optimality of the Euler-Maclaurin type quadrature and cubature formulas were studied, for instance, in works [4,8,9,12, 17,18,19,24].

The Euler-Maclaurin quadrature formulas can be viewed as well as an extension of the trapezoidal rule by the inclusion of correction terms. It should be noted that in applications and in solution of practical problems numerical integration formulas are interesting for functions with small smoothness.

The present paper is also devoted to extension of the trapezoidal rule.

We consider a quadrature formula of the form

$$\underset { 0 } { \int } ^ { 1 } \varphi ( x ) d x & \cong \sum _ { \beta = 0 } ^ { N } ( C _ { 0 } [ \beta ] \varphi ( h \beta ) + C _ { 1 } [ \beta ] \varphi ^ { \prime } ( h \beta ) ) \\ \intertext { a l i n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a nd a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a n d a }$$

where C 0 [ β ] are coefficients of the trapezoidal rule, i.e.

$$\begin{array} { c c } C _ { 0 } [ 0 ] & = \frac { h } { 2 }, \\ C _ { 0 } [ \beta ] & = h, \, \beta = 1, 2, \dots, N - 1, \\ C _ { 0 } [ N ] & = \frac { h } { 2 }, \end{array}$$

C 1 [ β ] are unknown coefficients of the formula (1) and they should be found, h = 1 N , N is a natural number. We suppose that an integrand ϕ belongs to W (2 1) , 2 (0 1), , where by W (2 1) , 2 (0 1) we denote the class of all functions , ϕ defined on [0 1] which , posses an absolutely continuous first derivative and whose second derivative is in L 2 (0 1). The class , W (2 1) , 2 (0 1) under the pseudo-inner product ,

$$\langle \varphi, \psi \rangle = \bigotimes _ { 0 } ^ { 1 } ( \varphi ^ { \prime \prime } ( x ) + \varphi ^ { \prime } ( x ) ) ( \psi ^ { \prime \prime } ( x ) + \psi ^ { \prime } ( x ) ) d x$$

is a Hilbert space if we identify functions that differ by a linear combination of a constant and e -x (see, for example, [1]). Here, in the Hilbert space W (2 1) , 2 (0 1), , we consider the corresponding norm

$$\| \varphi | W _ { 2 } ^ { ( 2, 1 ) } ( 0, 1 ) \| & = \left [ \int _ { 0 } ^ { 1 } ( \varphi ^ { \prime \prime } ( x ) + \varphi ^ { \prime } ( x ) ) ^ { 2 } d x \right ] ^ { 1 / 2 }. \\ \intertext { n o d f o r o n o c }$$

The difference

$$( \ell, \varphi ) = \bigotimes _ { 0 } ^ { 1 } \varphi ( x ) d x - \sum _ { \beta = 0 } ^ { N } ( C _ { 0 } [ \beta ] \varphi ( h \beta ) + C _ { 1 } [ \beta ] \varphi ^ { \prime } ( h \beta ) ) \\ \intertext { m a n d a n d a n d a n d }$$

is called the error and

$$\ell ( x ) = \varepsilon _ { [ 0, 1 ] } ( x ) - \sum _ { \beta = 0 } ^ { N } ( C _ { 0 } [ \beta ] \delta ( x - h \beta ) - C _ { 1 } [ \beta ] \delta ^ { \prime } ( x - h \beta ) ), \quad \quad ( 5 ) \\ \cdots \, \text{and set the same} \, \text{and set the } \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \text{$\delta$} \, \end{cases}$$

is said to be the error functional of the quadrature formula (1), where ε [0 1] , ( x ) is the indicator of the interval [0 1] and , δ is Dirac's delta function. The value ( /lscript, ϕ ) of the error functional /lscript at a function ϕ is defined as ( /lscript, ϕ ) = ∞ -∞ /lscript ( x ϕ x dx. ) ( )

∫ In order that the error functional (5) is defined on the space W (2 1) , 2 (0 1) it is , necessary to impose the following conditions for the functional /lscript

$$( \ell, 1 ) \coloneqq \frac { 1 } { \sum _ { \beta = 0 } ^ { N } C _ { 0 } [ \beta ] } = 0, \\ \vdots \frac { N } { \vdots }$$

$$( \ell, e ^ { - x } ) \coloneqq \, \underset { 0 } { \overset { 1 } { \int } } e ^ { - x } d x - \sum _ { \beta = 0 } ^ { N } ( C _ { 0 } [ \beta ] e ^ { - h \beta } - C _ { 1 } [ \beta ] e ^ { - h \beta } ) = 0. \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \cdots \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \text{$\quad \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \, $ \} }$$

The last two equations mean that the quadrature formula (1) is exact for any constant and e -x . We have chosen the coefficients C 0 [ β ], β = 0 1 , , ..., N such that the equality (6) is fulfilled. Therefore we have only condition (7) for coefficients C 1 [ β ], β = 0 1 , , ..., N .

The error functional /lscript of the formula (1) is a linear functional in W (2 1) , ∗ 2 (0 1), , where W (2 1) , ∗ 2 (0 1) is the conjugate space to the space , W (2 1) , 2 (0 1). ,

By the Cauchy-Schwarz inequality we have the following

$$| ( \ell, \varphi ) | \leq \| \varphi | W _ { 2 } ^ { ( 2, 1 ) } ( 0, 1 ) \| \cdot \| \ell | W _ { 2 } ^ { ( 2, 1 ) * } ( 0, 1 ) \|.$$

Hence we conclude that the error (4) of the formula (1) is estimated by the norm

$$\left \| \ell | W _ { 2 } ^ { ( 2, 1 ) * } ( 0, 1 ) \right \| = \sup _ { \left \| \varphi | W _ { 2 } ^ { ( 2, 1 ) } ( 0, 1 ) \right \| = 1 } | ( \ell, \varphi ) | & & ( 8 ) \\ \text{or functional } ( 5 ).$$

of the error functional (5).

The main aim of this work is to find the minimum of the absolute value of the error (4) by coefficients C 1 [ β ] for given C 0 [ β ] in the space W (2 1) , 2 . That is the problem is to find the coefficients C 1 [ β ] that satisfy the following equality

$$\left \| \mathring { \ell } | W _ { 2 } ^ { ( 2, 1 ) * } \right \| = \inf _ { C _ { 1 } [ \beta ] } \left \| \ell | W _ { 2 } ^ { ( 2, 1 ) * } \right \|.$$

The coefficients C 1 [ β ] which satisfy the last equation are called optimal and are denoted as ˚ [ C 1 β ].

Thus, to obtain the optimal quadrature formula of the form (1) in the sense of Sard in the space W (2 1) , 2 (0 1), we need to solve the following problems. ,

Problem 1. Find the norm of the error functional (5) of the quadrature formula (1) in the space W (2 1) , ∗ 2 .

Problem 2. Find the coefficients ˚ [ C 1 β ] that satisfy equality (9).

Here we solve Problems 1 and 2 by Sobolev's method using the discrete analog of the differential operator d 2 dx 2 -1.

The paper is organized as follows: in Section 2 using the extremal function of the error functional /lscript the norm of this functional is calculated, i.e. Problem 1 is solved; Section 3 is devoted to solution of Problem 2. Here the system of linear equations for the coefficients C 1 [ β ] of the optimal quadrature formulas (1) is obtained in the space W (2 1) , 2 (0 1). In Subsection 3.1 using the discrete analog , of the operator d 2 dx 2 -1 the explicit formulas for the coefficients C 1 [ β ] of optimal quadrature formula of the form (1) are obtained. Furthermore, it is proved that the obtained quadrature formula of the form (1) is exact for any function of the set F = span 1 { , x, e x , e -x } . Finally, in Subsection 3.2 in the space W (2 1) , 2 the square of the norm of the error functional of the constructed quadrature formula is calculated. It is shown that the error of the obtained optimal quadrature formula is less than the error of the Euler-Maclaurin quadrature formula on the space L (2) 2 .

## 2 The norm of the error functional (5)

In this section we study Problem 1. To calculate the norm of the error functional (5) in the space W (2 1) , ∗ 2 (0 1) we use , the extremal function for this functional (see, [21,22]) which satisfies the equality

$$( \ell, \psi _ { \ell } ) & = \left \| \ell | W _ { 2 } ^ { ( 2, 1 ) * } \right \| \cdot \left \| \psi _ { \ell } | W _ { 2 } ^ { ( 2, 1 ) } \right \|. \\ \colon \text{in } [16 ] & \text{for a linear functional } \ell \text{ defined on } t$$

∥ ∥ ∥ ∥ We note that in [16] for a linear functional /lscript defined on the Hilbert space W ( m,m -1) 2 the extremal function was found and it was shown that the extremal function ψ /lscript is the solution of the boundary value problem

$$\psi _ { \ell } ^ { ( 2 m ) } ( x ) - \psi _ { \ell } ^ { ( 2 m - 2 ) } ( x ) = ( - 1 ) ^ { m } \ell ( x ),$$

$$& \tilde { \left ( \psi _ { \ell } ^ { ( m + s ) } ( x ) - \psi _ { \ell } ^ { ( m + s - 2 ) } ( x ) \right ) | _ { x = 0 } ^ { x = 1 } } = 0, \ \ s = 0, 1, \dots, m - 1, \quad \quad ( 1 1 ) \\ & \left ( \psi _ { \ell } ^ { ( m ) } ( x ) + \psi _ { \ell } ^ { ( m - 1 ) } ( x ) \right ) | _ { x = 0 } ^ { x = 1 } } = 0.$$

That is for the extremal function ψ /lscript the following was proved.

$$\overset { \cdot } { \left ( \psi _ { \ell } ^ { ( m ) } ( x ) + \psi _ { \ell } ^ { ( m - 1 ) } ( x ) \right ) } | _ { x = 0 } ^ { x = 1 } = 0. & & ( 1 2 ) \\ \text{That is for the extremal function $\psi_{e}$ the following was proved.}$$

Theorem 1 (Theorem 2.1 of [16]) The solution of the boundary value problem (10)-(12) is the extremal function ψ /lscript of the error functional /lscript and has the following form

$$\psi _ { \ell } ( x ) = ( - 1 ) ^ { m } \ell ( x ) * G _ { m } ( x ) + P _ { m - 2 } ( x ) + d e ^ { - x },$$

where

$$G _ { m } ( x ) = \frac { \text{sgn} x } { 2 } \left ( \frac { e ^ { x } - e ^ { - x } } { 2 } - \sum _ { k = 1 } ^ { m - 1 } \frac { x ^ { 2 k - 1 } } { ( 2 k - 1 )! } \right ) \quad \quad ( 1 3 )$$

is the solution of the eqution G (2 m ) m ( x ) -G (2 m -2) m ( x ) = δ ( x , d ) is any real number and P m -2 ( x ) is a polynomial of degree m -2 .

Furthermore, there were shown that ‖ | /lscript W ( m,m -1) 2 ‖ = ‖ ψ /lscript | W ( m,m -1) ∗ 2 ‖ and

$$( \ell, \psi _ { \ell } ) = \left \| \ell | W _ { 2 } ^ { ( m, m - 1 ) * } \right \| ^ { 2 }. \quad & ( 1 4 ) \\ 1, \text{ in the case } m = 2, \, \text{ we get the extremal function } \psi _ { \ell } \text{ for the } \,.$$

From Theorem 1, in the case m = 2, we get the extremal function ψ /lscript for the error functional (5) and it has the form

$$\psi _ { \ell } ( x ) = \ell ( x ) * G _ { 2 } ( x ) + p _ { 0 } + d e ^ { - x },$$

where

$$G _ { 2 } ( x ) = \frac { s g n x } { 2 } \left ( \frac { e ^ { x } - e ^ { - x } } { 2 } - x \right ),$$

d and p 0 are any real numbers.

Then, from (14), in the case m = 2, using (5) and (15), taking into account equations (6) and (7), we get

$$| | \ell | | ^ { 2 } = ( \ell, \psi _ { \ell } ) = \sum _ { \beta = 0 } ^ { N } \sum _ { \gamma = 0 } ^ { N } \left ( C _ { 0 } [ \beta ] C _ { 0 } [ \gamma ] G _ { 2 } ( h \beta - h \gamma ) - C _ { 1 } [ \beta ] C _ { 1 } [ \gamma ] G _ { 2 } ^ { \prime \prime } ( h \beta - h \gamma ) \right ) + \\ + 2 \sum _ { \beta = 0 } ^ { N } C _ { 1 } [ \beta ] \left ( \int _ { 0 } ^ { 1 } G _ { 2 } ^ { \prime } ( x - h \beta ) d x + \sum _ { \gamma = 0 } ^ { N } C _ { 0 } [ \gamma ] G _ { 2 } ^ { \prime } ( h \beta - h \gamma ) \right ) - \\ - 2 \sum _ { \beta = 0 } ^ { N } C _ { 0 } [ \beta ] \int _ { 0 } ^ { 1 } G _ { 2 } ( x - h \beta ) d x + \int _ { 0 } ^ { 1 } \int _ { 0 } ^ { 1 } G _ { 2 } ( x - y ) d x d y, \\ \text{where } G _ { 2 } ( x ) \text{ is defined by } ( 1 6 ), G _ { 2 } ^ { \prime } ( x ) \text{ and } G _ { 2 } ^ { \prime \prime } ( x ) \text{ are derivatives of } G _ { 2 } ( x ), i. e.$$

where G 2 ( x ) is defined by (16), G ′ 2 ( x ) and G ′′ 2 ( x ) are derivatives of G 2 ( x ), i.e.

$$G _ { 2 } ^ { \prime } ( x ) = \frac { \text{sgnx} } { 2 } \left ( \frac { e ^ { x } + e ^ { - x } } { 2 } - 1 \right ) \text{ and } G _ { 2 } ^ { \prime \prime } ( x ) = \frac { \text{sgnx} } { 2 } \left ( \frac { e ^ { x } - e ^ { - x } } { 2 } \right ). \quad ( 1 8 )$$

It is easy to see from (13) and (18) that

$$G _ { 1 } ( x ) = G _ { 2 } ^ { \prime \prime } ( x ).$$

Thus Problem 1 is solved.

In the next section we study Problem 2.

## 3 Minimization of the norm (17)

Now we consider the minimization problem of the expression (17) by the coefficients C 1 [ β ] under the condition (7). For this we use the Lagrange method of conditional extremum.

Consider the Lagrange function

$$\Psi ( C _ { 1 } [ 0 ], C _ { 1 } [ 1 ], \dots, C _ { 1 } [ N ], d ) = \| \ell \| ^ { 2 } + 2 d ( \ell, e ^ { - x } ).$$

Taking partial derivatives from the function Ψ by C 1 [ β ], β = 0 1 , , ..., N then equating them to 0 and using the condition (7), we get the following system of N +2 linear equations with N +2 unknowns

$$\sum _ { \gamma = 0 } ^ { N } C _ { 1 } [ \gamma ] G _ { 2 } ^ { \prime \prime } ( h \beta - h \gamma ) + d e ^ { - h \beta } = F ( h \beta ), \ \beta = 0, 1, 2, \dots, N, \quad \ ( 2 0 )$$

$$\sum _ { \gamma = 0 } ^ { \cdot } C _ { 1 } [ \gamma ] e ^ { - h \gamma } = g,$$

where

$$F ( h \beta ) = \int _ { 0 } ^ { 1 } G _ { 2 } ^ { \prime } ( x - h \beta ) d x + \sum _ { \gamma = 0 } ^ { N } C _ { 0 } [ \gamma ] G _ { 2 } ^ { \prime } ( h \beta - h \gamma ),$$

$$g = e ^ { - 1 } - 1 + \sum _ { \gamma = 0 } ^ { N } C _ { 0 } [ \gamma ] e ^ { - h \gamma }. \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \quad.$$

Here C 0 [ γ ], γ = 0 1 , , ..., N are defined by (2), C 1 [ β ], β = 0 1 , , ..., N and d are unknowns.

The system (20)-(21) has a unique solution for any fixed natural number N and this solution gives the minimum to the expression (17). Here we omit the proof of the existence and uniqueness of the solution of this system. These statements can be proved similarly as the proof of the existence and uniqueness of the solution of the discrete Wiener-Hopf type system for the optimal coefficients of quadrature formulas with the form ∫ 1 0 f ( x dx ) ∼ = ∑ N β =0 C β f β [ ] [ ] in the space L ( m ) 2 (0 1) (see , [20,21,22]).

## 3.1. The coefficients of the optimal quadrature formula (1)

In this subsection we solve the system (20)-(21) and find the explicit forms for optimal coefficients C 1 [ β ], β = 0 1 , , ..., N .

Here we use the concept of discrete argument functions and operations on them. The theory of discrete argument functions is given in [21,22]. We give some definitions about functions of discrete argument.

Suppose that ϕ and ψ are real-valued functions of real variable x and are defined in real line R .

A function ϕ hβ ( ) is called a function of discrete argument if it is defined on some set of integer values of β .

The inner product of two discrete functions ϕ hβ ( ) and ψ hβ ( ) is defined as the following number

$$& \text{er} \\ & [ \varphi ( h \beta ), \psi ( h \beta ) ] = \sum _ { \beta = - \infty } ^ { \infty } \varphi ( h \beta ) \cdot \psi ( h \beta ), \\ & \text{the ri\shorthat hand side of the last equality converse $s$}$$

if the series on the right hand side of the last equality converges absolutely.

The convolution of two discrete argument functions ϕ hβ ( ) and ψ hβ ( ) is the inner product

$$\varphi ( h \beta ) * \psi ( h \beta ) = [ \varphi ( h \gamma ), \psi ( h \beta - h \gamma ) ] = \sum _ { \gamma = - \infty } ^ { \infty } \varphi ( h \gamma ) \cdot \psi ( h \beta - h \gamma ).$$

Furthermore, for finding the coefficients C 1 [ β ] of the optimal quadrature formula (1) we need the discrete analog of the differential operator d 2 dx 2 -1. It should be noted that in the work [14] the discrete analog of the differential operator d 2 m dx 2 m -d 2 m -2 dx 2 m -2 was constructed. In particular, when m = 1 from the result of the work [14] we get the following

Theorem 2 The discrete analog D 1 ( hβ ) of the differential operator d 2 dx 2 -1 satisfying the equation

$$D _ { 1 } ( h \beta ) * G _ { 1 } ( h \beta ) = \delta _ { d } ( h \beta )$$

has the form

$$m \\ D _ { 1 } ( h \beta ) = \frac { 1 } { 1 - e ^ { 2 h } } \begin{cases} 0, & | \beta | \geq 2, \\ - 2 e ^ { h }, & | \beta | = 1, \\ 2 ( 1 + e ^ { 2 h } ), \, \beta = 0, \end{cases} \\ h \beta ) = \frac { \text{sgn} ( h \beta ) } { 2 } \left ( \frac { e ^ { h \beta } - e ^ { - h \beta } } { 2 } \right ) \ a n d \ \delta _ { d } ( h \beta ) = \begin{cases} 1, \beta = 0, \\ 0, \beta \neq 0. \end{cases}$$

/negationslash

Furthermore, it is easy to check that

 where G 1 ( hβ ) = sgn( hβ ) 2 ( e hβ -e -hβ 2 ) and δ d ( hβ ) = { 1 , β = 0 , 0 , β = 0 .

$$D _ { 1 } ( h \beta ) * e ^ { h \beta } = 0 \text{ and } D _ { 1 } ( h \beta ) * e ^ { - h \beta } = 0.$$

Now we turn to get the solution of the system (20)-(21) using (24). Suppose C 1 [ β ] = 0 when β &lt; 0 and β &gt; N . Then we rewrite the system (20)-(21) in the following convolution form

$$C _ { 1 } [ \beta ] * G _ { 2 } ^ { \prime \prime } ( h \beta ) + d e ^ { - h \beta } = F ( h \beta ), \ \beta = 0, 1, \dots, N, \quad \quad ( 2 6 )$$

$$\sum _ { \gamma = 0 } ^ { \, ^ { \ast \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \rho \, \rho } } ^ { \, ^ { \ast \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \rho \, \rho } } ^ { \, ^ { \ast \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rangle } } ^ { \, ^ { \ast \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho } } ^ { \, ^ { \ast \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \,\nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \, \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \rho \langle \, \rangle \, \rangle \, \rangle \, \rangle \, \rangle \, \rangle \, \rangle \, \rangle \, \rangle \, \rangle \, \rangle \, \rangle \, \rangle \, \rangle \, \rangle \, \rangle \rangle \, \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \langle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle \rangle$$

Here, calculating the right hand sides of equalities (22) and (23), for F hβ ( ) and g we get

$$F ( h \beta ) = \sum _ { 1 } ^ { \ } \left ( \frac { e ^ { h \beta } } { 8 } ( e ^ { - 1 } + 1 ) - \frac { e ^ { - h \beta } } { 8 } ( e + 1 ) \right ) \left ( \frac { h ( e ^ { h } + 1 ) } { e ^ { h } - 1 } - 2 \right ), \quad \ \left ( 2 8 \right )$$

$$g = \frac { \lambda } { 2 } ( 1 - e ^ { - 1 } ) \left ( \frac { h ( e ^ { h } + 1 ) } { e ^ { h } - 1 } - 2 \right ). \begin{array} { c c c } \lambda & \lambda & \lambda \\ & \lambda & \lambda & \lambda \\ & & \lambda & \lambda \\ \end{array}$$

Taking into account (18) and (19), using Theorem 2 we get

$$D _ { 1 } ( h \beta ) * G _ { 2 } ^ { \prime \prime } ( h \beta ) = \delta _ { d } ( h \beta ),$$

Denoting by

$$\dot { \nu } _ { 0 } u ( h \beta ) = C _ { 1 } [ \beta ] * G _ { 2 } ^ { \prime \prime } ( h \beta ) + d e ^ { - h \beta }$$

the left hand side of the equation (26) we get

$$C _ { 1 } [ \beta ] = D _ { 1 } ( h \beta ) * u ( h \beta ).$$

Indeed, if the function u hβ ( ) is defined at all integer values of β , then taking into account Theorem 2 and using properties (25) of the function D 1 ( hβ ), we have

$$D _ { 1 } ( h \beta ) * u ( h \beta ) & = D _ {$$

Thus, if we find the function u hβ ( ) for all integer values of β then the optimal coefficients C 1 [ β ] will be found from the equality (32).

The following is true.

Theorem 3 The coefficients of the optimal quadrature formula of the form (1) in the sense of Sard in the space W (2 1) , 2 (0 1) , have the following form

$$\begin{array} {$$

Proof. From equality (26) taking into account (31) we get that

$$u ( h \beta ) = F ( h \beta ) \text{ for } \beta = 0$$

Now we find the function u hβ ( ) for β &lt; 0 and β &gt; N .

Let β ≤ 0 then from (31), using the form (18) of the function G ′′ 2 ( x ) and equality (27), we have

$$u ( h \beta ) = - \frac { 1 } { 4 } e ^ { h \beta$$

Similarly, for β ≥ N we obtain

$$u ( h \beta ) = \frac { 1 } { 4 } e ^ { h \beta }$$

Then, keeping in mind the last two equalities and denoting by

$$D = \frac { 1 } { 4 } \sum _ { \gamma = 0 } ^ { N } C _ { 1 } [ \gamma ] e ^ { h \gamma },$$

for u hβ ( ) we get the following

$$\cdot \quad \cdot \quad \cdot \quad \cdot \\ u ( h \beta ) = \begin{cases} - \frac { 1 } { 4 } e ^ { h \beta } g + ( d + D ) e ^ { - h \beta }, \, \beta \leq 0, \\ F ( h \beta ), & 0 \leq \beta \leq N, \\ \frac { 1 } { 4 } e ^ { h \beta } g + ( d - D ) e ^ { - h \beta }, & \beta \geq N. \end{cases} \\ \text{in the equality (34), $d$ and $D$ are unknowns. These unknowns can be found} \\ \text{the conditions of consistency of values of the function $u(h\beta)$ at the points}$$

Here, in the equality (34), d and D are unknowns. These unknowns can be found from the conditions of consistency of values of the function u hβ ( ) at the points β = 0 and β = N . Therefore from (34) when β = 0 and β = N we obtain the system of linear equations for d and D . Then, using (28) and (29), after some calculations, we have

$$\begin{cases} d = 0, & ( 3 5 ) \\ D = \frac { 1 } { 8 } \left ( \frac { h ( e ^ { h } + 1 ) } { e ^ { h } - 1 } - 2 \right ) ( 1 - e ). \end{cases} \\ \mathbb { Z } \text{2} \right ) \text{for} \, \beta = 0, 1, \dots, N, \, \text{using} \, ( 2 4 ) \, \text{and} \, ( 3 4 ) \, \text{and taking into account} \\ \text{v calculations we set} \, ( 3 3 ) \text{ Theorem $3$ is nrowed} \quad \text{$\sqcap$}$$

Finally, from (32) for β = 0 1 , , ..., N , using (24) and (34) and taking into account (35), by directly calculations we get (33). Theorem 3 is proved. ✷

Remark 1. Using (2) and (33), one can get that ( /lscript, x ) = 0 and ( /lscript, e x ) = 0. These equalities mean that the optimal quadrature formula of the form (1) with the coefficients (2) and (33) is also exact for functions x and e x . Therefore, keeping in mind equalities (6) and (7), we conclude that the optimal quadrature formula of the form (1.1) with coefficients (2) and (33) is exact for any linear combinations of functions 1 , x, e x and e -x , i.e. it is exact for elements of the set F = span 1 { , x, e x , e -x } .

## 3.2. The norm of the error functional ˚ /lscript of the optimal quadrature formula (1)

In this subsection we study the order of convergence of the optimal quadrature formula of the form (1) with coefficients (2) and (33), i.e. we calculate the square of the norm (17) of the error functional for the optimal quadrature formula (1). The following holds

Theorem 4 Square of the norm of the error functional (5) for the optimal quadrature formula (1) with coefficients (2) and (33) on the space W (2 1) , 2 (0 1) , has the form

$$\left \| \hat { \ell } | W _ { 2 } ^ { ( 2, 1 ) * } ( 0, 1 ) \right \| ^ { 2 } & = - \sum _ { n = 4 } ^ { \infty } \frac { B _ { n } h ^ { n } } { n! } \\ & = \frac { 1 } { 7 2 0 } h ^ { 4 } - \frac { 1 } { 3 0 2 4 0 } h ^ { 6 } + O ( h ^ { 8 } ), \\ _ { - } \quad _ { - } \quad.$$

where B n are Bernoulli numbers.

Proof. We rewrite the expression (17) as follows

$$\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ \| \hat { \ell } \| ^ { 2 } & = - \sum _ { \beta = 0 } ^ { N } \hat { C } _ { 1 } [ \beta ] \left ( \sum _ { \gamma = 0 } ^ { N } \hat { C } _ { 1 } [ \gamma ] G _ { 2 } ^ { \prime \prime } ( h \beta - h \gamma ) - F ( h \beta ) \right ) + \sum _ { \beta = 0 } ^ { N } \hat { C } _ { 1 } [ \beta ] F ( h \beta ) \\ & \quad + \sum _ { \beta = 0 } ^ { N } \sum _ { \gamma = 0 } ^ { N } C _ { 0 } [ \beta ] C _ { 0 } [ \gamma ] G _ { 2 } ( h \beta - h \gamma ) - 2 \sum _ { \beta = 0 } ^ { N } C _ { 0 } [ \beta ] \int _ { 0 } ^ { 1 } G _ { 2 } ( x - h \beta ) d x + \int _ { 0 } ^ { 1 } \int _ { 0 } ^ { 1 } G _ { 2 } ( x - y ) d x d y, \\ \text{where } F ( h \beta ) \text{ is defined by } ( 2 2 ).$$

where F hβ ( ) is defined by (22).

Hence, taking into account (35), we have

$$\| \mathring { \ell } \| ^ { 2 } = A _ { 1 } + A _ { 2 } - 2 A _ { 3 } + A _ { 4 },$$

here

$$\stackrel { \dots } { \dots } A _ { 1 } & = \sum _ { \beta = 0 } ^ { N } \stackrel { \ast } { C } _ { 1 } [ \beta ] \ F ( h \beta ), \quad \ A _ { 2 } = \sum _ { \beta = 0 } ^ { N } \sum _ { \gamma = 0 } ^ { N } C _ { 0 } [ \beta ] C _ { 0 } [ \gamma ] G _ { 2 } ( h \beta - h \gamma ), \\ A _ { 3 } & = \sum _ { \beta = 0 } ^ { N } C _ { 0 } [ \beta ] \int _ { 0 } ^ { 1 } G _ { 2 } ( x - h \beta ) d x, \, A _ { 4 } = \int _ { 0 } ^ { 1 } \int _ { 0 } ^ { 1 } G _ { 2 } ( x - y ) d x d y. \\ \text{Now we need the following sums which are obtained by using (2) and (33)}$$

Now we need the following sums which are obtained by using (2) and (33)

$$\lambda \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \\ \sum _ { \beta = 0 } ^ { N } C _ { 0 } [ \beta ] = 1, \ \sum _ { \beta = 0 } ^ { N } C _ { 0 } [ \beta ] ( h \beta ) = \frac { 1 } { 2 }, \ \sum _ { \beta = 0 } ^ { N } C _ { 0 } [ \beta ] ( h \beta ) ^ { 2 } = \frac { h ^ { 2 } } { 6 } + \frac { 1 } { 3 }, \\ \sum _ { \beta = 0 } ^ { N } \hat { C } _ { 1 } [ \beta ] e ^ { - h \beta } = \frac { 1 } { 2 } ( 1 - e ^ { - 1 } ) \left ( \frac { h ( e ^ { h } + 1 ) } { e ^ { h - 1 } } - 2 \right ), \ \sum _ { \beta = 0 } ^ { N } \hat { C } _ { 1 } [ \beta ] e ^ { h \beta } = \frac { 1 } { 2 } ( 1 - e ) \left ( \frac { h ( e ^ { h } + 1 ) } { e ^ { h - 1 } } - 2 \right ). \\ \text{Taking into account (28) and (16), using (38) for } A _ { 1 }, A _ { 2 }, A _ { 3 } \text{ and } A _ { 4 } \text{ we get}$$

Taking into account (28) and (16), using (38) for A 1 , A 2 , A 3 and A 4 we get

$$A _ { 1 } & = \frac { 1 - e ^ { 2 } } { 8 e } \left ( \frac { h ( e ^ { h } + 1 ) } { e ^ { h } - 1 } - 2 \right ) ^ { 2 }, \, A _ { 2 } & = \frac { h ^ { 2 } ( e ^ { h } + 1 ) ^ { 2 } ( e ^ { 2 } - 1 ) } { 8 e ( e ^ { h } - 1 ) ^ { 2 } } - \frac { h ^ { 2 } + 2 } { 1 2 } - \frac { h ( e ^ { h } + 1 ) } { 2 ( e ^ { h } - 1 ) }, \\ A _ { 3 } & = \frac { h ( e ^ { h } + 1 ) ( e ^ { 2 } - 1 ) } { 4 e ( e ^ { h } - 1 ) }, \, \quad \, A _ { 4 } & = \frac { e ^ { 2 } - 1 } { 2 e } - \frac { 7 } { 6 }.$$

Further, putting the last equalities to (37) and after some simplifications we have

$$\| \mathring { \ell } \| ^ { 2 } = 1 - \frac { 1 } { 2 } h + \frac { 1 } { 1 2 } h ^ { 2 } - \frac { h } { e ^ { h } - 1 }.$$

Hence, using well known formula x e x -1 = ∞ ∑ n =0 B n n ! x n , | x &lt; | 2 π , we get (36). Theorem 4 is proved

✷

Remark 2. It should be noted that optimality of the classical Euler-Maclaurin was proved and the error of this quadrature formula was calculated in L ( m ) 2 , where L ( m ) 2 is the space of functions which are square integrable with m -th generalized derivative (see, for instance, [4,12,17]). In particular, when m = 2 from Corollary 5.1 of the work [17] we get optimality of the Euler-Maclaurin formula

$$\int _ { 0 } ^ { 1 } \varphi ( x ) d x \cong h \left ( \frac { 1 } { 2 } \varphi ( 0 ) + \varphi ( h ) + \varphi ( 2 h ) + \dots + \varphi ( h ( N - 1 ) ) + \frac { 1 } { 2 } \varphi ( 1 ) \right ) + \frac { h ^ { 2 } } { 1 2 } ( \varphi ^ { \prime } ( 0 ) - \varphi ^ { \prime } ( 1 ) )$$

in the space L (2) 2 (0 1). Furthermore for the square of the norm of the error func-, tional the following is valid

$$\| \mathring { \ell } | L _ { 2 } ^ { ( 2 ) * } ( 0, 1 ) \| ^ { 2 } = \frac { h ^ { 4 } } { 7 2 0 }.$$

Comparison of equalities (36) and (40) shows that the error of the optimal quadrature formula of the form (1.1) on the space W (2 1) , 2 (0 1) is less than the error of , the Euler-Maclaurin quadrature formula (39) on the space L (2) 2 (0 1). ,

## Acknowledgments

This work has been done while A.R.Hayotov was visiting Department of Mathematical Sciences at KAIST, Daejeon, Republic of Korea. A.R.Hayotov is very grateful to professor Chang-Ock Lee and his research group for hospitality. A.R. Hayotov's work was supported by the 'Korea Foundation for Advanced Studies'/'Chey Institute for Advanced Studies' International Scholar Exchange Fellowship for academic year of 2018-2019

## References

- 1. J.H. Ahlberg, E.N. Nilson, J.L. Walsh, The Theory of Splines and Their Applications, Academic Press, New York - London, 1967.
- 2. I. Babuˇka, Optimal quadrature formulas (Russian), Dokladi Akad. Nauk SSSR. 149 (1963) s 227-229.
- 3. P. Blaga, Gh. Coman, Some problems on optimal quadrature, Stud. Univ. Babe¸ s-Bolyai Math. 52, no. 4 (2007) 21-44.
- 4. T. Catina¸ s, Gh. Coman, Optimal quadrature formulas based on the φ -function method, Stud. Univ. Babe¸ s-Bolyai Math. 51, no. 1 (2006) 49-64.
- 5. A.R. Hayotov, G.V. Milovanovi´ c, Kh.M. Shadimetov, On an optimal quadrature formula in the sense of Sard. Numerical Algorithms, v.57, no. 4, (2011) 487-510.
- 6. A.R. Hayotov, G.V. Milovanovi´, Kh.M. Shadimetov, Optimal quadratures in the sense of c Sard in a Hilbert space. Applied Mathematics and Computation, 259 (2015) 637-653.
- 7. P. K¨hler, On the weights of Sard's quadrature formulas, Calcolo, 25 (1988) 169-186. o
- 8. F. Lanzara, On optimal quadrature formulae, J. Ineq. Appl. 5 (2000) 201-225.
- 9. C.A. Micchelli, Best quadrature formulas at equally spaced nodes, J. Math. Anal. Appl. 47 (1974) 232-249.
- 10. S.M. Nikolskii, Quadrature Formulas, Nauka, Moscow, 1988.(in Russian).
- 11. A. Sard, Best approximate integration formulas; best approximation formulas, Amer. J. Math. 71 (1949) 80-91.
- 12. I.J. Schoenberg, On monosplines of least deviation and best quadrature formulae, J. Soc. Indust. Appl. Math. Ser. B Numer. Anal. 2 (1965) 144-170.
- 13. I.J. Schoenberg, S.D. Silliman, On semicardinal quadrature formulae. Math. Comp. 28 (1974) 483-497.
- 14. Kh.M. Shadimetov, A.R. Hayotov, Construction of the discrete analogue of the differential operator d 2 m dx 2 m -d 2 m -2 dx 2 m -2 , Uzbek mathematical journal, 2004, no.2, pp. 85-95.
- 16. Kh.M. Shadimetov, A.R. Hayotov, Optimal quadrature formulas in the sense of Sard in W ( m,m -1) 2 space, Calcolo 51 (2014) 211-243.
- 15. Kh.M. Shadimetov, A.R. Hayotov, Optimal quadrature formulas with positive coefficients in L ( m ) 2 (0 , 1) space, J. Comput. Appl. Math. 235 (2011) 1114-1128.
- 17. Kh.M. Shadimetov, A.R. Hayotov, F.A. Nuraliev, On an optimal quadrature formula in Sobolev space L ( m ) 2 (0 , 1), J. Comput. Appl. Math. 243 (2013) 91-112.
- 18. Kh.M. Shadimetov, A.R. Hayotov, F.A. Nuraliev, Optimal quadrature formulas of EulerMaclaurin type, Applied Mathematics and Computation 276 (2016) 340-355.

- 19. Kh.M. Shadimetov, F.A. Nuraliev, Optimal formulas of numerical integration with derivatives in Sobolev space, Journal of Siberian Federal University. Math. and Phys. 2018, 11 (6) 764-775.
- 20. S.L. Sobolev, The coefficients of optimal quadrature formulas, Selected Works of S.L. Sobolev, Springer, (2006) 561-566.
- 21. S.L. Sobolev, Introduction to the Theory of Cubature Formulas (Russian), Nauka, Moscow, 1974.
- 22. S.L. Sobolev, V.L. Vaskevich, The Theory of Cubature Formulas, Kluwer Academic Publishers Group, Dordrecht, 1997.
- 23. F.Ya. Zagirova, On construction of optimal quadrature formulas with equal spaced nodes (Russian). Novosibirsk (1982), 28 p. (Preprint No. 25, Institute of Mathematics SD of AS of USSR)
- 24. A.A. Zhensikbaev, Monosplines of minimal norm and the best quadrature formulas, Uspekhi Mat. Nauk. 1981, 36, 107-159.
Abdullo R. Hayotov Rashidjon G. Rasulov
The order of convergence of an optimal quadrature formula with derivative in the space $W_2^{(2,1)}$
The present work is devoted to extension of the trapezoidal rule in the space
$W_2^{(2,1)}$. The optimal quadrature formula is obtained by minimizing the
error of the formula by coefficients at values of the first derivative of a
integrand. Using the discrete analog of the operator $\frac{d^2}{dx^{2}}-1$ the
explicit formulas for the coefficients of the optimal quadrature formula are
obtained. Furthermore, it is proved that the obtained quadrature formula is
exact for any function of the set
$\mathbf{F}=\mathrm{span}\{1,x,e^{x},e^{-x}\}$. Finally, in the space
$W_2^{(2,1)}$ the square of the norm of the error functional of the constructed
quadrature formula is calculated. It is shown that the error of the obtained
optimal quadrature formula is less than the error of the Euler-Maclaurin
quadrature formula on the space $L_2^{(2)}$.