There are several methods of construction of optimal quadrature formulas in the sense of Sard such as the spline method, ϕ-function method (see e.g. [3,9]) and Sobolev’s method which is based on construction of discrete analogs of a linear differential operator (see e.g. [19,20]) .