Numerical integration rules based on B-Spline bases (Q6580299)

It assumed that \(f \in \mathcal{C}^{p+1}(\mathbb{R})\), \(h\) a constant and \(B_p(x)\) the \(p\)-degree B-spline supported on \(I_p= \Big[-\frac{p+1}{2}, \frac{p+1}{2}\Big]\). Let the vector \(f_{n,p}=\big(f_{n-\lfloor \frac{p}{2}\rfloor},\dots, f_{n+\lfloor \frac{p}{2}\rfloor}\big)\) be defined on the equidistant knots \(S_p=\Big\{ -\frac{p+1}{2},\dots, \frac{p+1}{2}\Big\} \), with \(f_i=f(ih)\) and where \(\lfloor \cdot \rfloor \) is the floor function. The function \(f\) approximates by the quasi-interpolation operator of the form\N\[\NQ_p(f)(x)=\sum_{n \in \mathbb{Z}} L_p(f_{n,p}) B_p(\frac{x}{h}-n), \quad \text{where} \quad L_p(f_{n.p})=\sum_{j=-\lfloor \frac{p}{2}\rfloor}^{\lfloor \frac{p}{2}\rfloor} c_{p,j} f_{n+j}.\N\]\NThe following formula\N\[\N\int\limits_a^b f(x) d x =h \sum_{i=-2\lfloor \frac{p}{2}\rfloor}^{-1} \xi_{p,i} (f^a_i+f^a_{N-i}-f^a_{-i}-f^a_{N-i})+\frac{h}{2}(f^a_0+f^a_N)+h \sum_{i=1}^{N-1} f^a_i,\N\]\Nwhere \(f^a_i=f(a+ih)\) is obtained in the article. The table of coefficients \(\xi_{p,i}\) is given in the paper. The resulting rules are defined as a perturbation of the trapezoidal integration method. Examples are considered in the article.