## Univariate spline quasi-interpolants and applications to numerical analysis

### Paul Sablonnière INSA and IRMAR, Rennes

Abstract
We describe some new univariate spline quasi-interpolants on uniform partitions of bounded intervals. Then we give some applications to numerical analysis: integration, differentiation and approximation of zeros. AMS classification: 41A15, 65D07, 65D25, 65D32.

1 Introduction

Univariate spline quasi-interpolants (abbr. QIs) can be defined as operators of the form $Qf={\sum }_{j\in J}{\mu }_{j}\left(f\right){B}_{j}$  where $\left\{{B}_{j},j\in J\right\}$  is the B-spline basis of some space of splines, say of degree $d$  , on a bounded interval $I=\left[a,b\right]$  endowed with some partition ${X}_{n}$  of $I$  in $n$  subintervals. We denote by ${\Pi }_{d}$  the space of polynomials of total degree at most $d$  . In general we impose that $Q$  is exact on the space ${\Pi }_{d}$  , i.e. $Qp=p$  for all $p\in {\Pi }_{d}$  . Some authors impose further that $Q$  is a projector on the space of splines itself (see e.g.[6, [11,[12). As a consequence of this property, the approximation order is $O\left({h}^{d+1}\right)$  on smooth functions, $h$  being the maximum steplength of the partition ${X}_{n}$  . The coefficients ${\mu }_{j}$  are local linear functionals which are in general of one of the following types:
(i) differential type : ${\mu }_{j}\left(f\right)$  is a linear combination of values of derivatives of $f$  , of order at most $d$  , at some point in $supp\left({B}_{j}\right)$  (see e.g. [1, [2). The associated quasi-interpolant is called a differential quasi-interpolant (abbr. DQI). (ii) integral type : ${\mu }_{j}\left(f\right)$  is a linear combination of weighted mean values of $f$  , i.e. of quantities ${\int }_{a}^{b}f{w}_{j}$  where ${w}_{j}$  can be, for example, a linear combination of B-splines (see e.g. [1, [2, [18). The associated quasi-interpolant is called an integral quasi-interpolant (abbr. iQI).
(iii) discrete type : ${\mu }_{j}\left(f\right)$  is a linear combination of discrete values of $f$  at some points in the neighbourhood of $supp\left({B}_{j}\right)$  (see e.g. [3, [11, [12, [14). The associated quasi-interpolant is called a discrete quasi-interpolant (abbr. dQI).
The main advantage of QIs is that they have a direct construction without solving any system of linear equations. Moreover,they are local, in the sense that the value of $Qf\left(x\right)$  depends only on values of $f$  in a neighbourhood of $x$  . Finally, they have a rather small infinity norm, so they are nearly optimal approximants.
In this paper, we only consider dQs, neither DQIs nor iQIs. We also restrict our study to splines defined on uniform partitions of $I=\left[a,b\right]$  . Our aim is to give explicit formulas for dQIs of degrees $2\le d\le 5$  and some applications to three classical problems in numerical analysis: approximate integration and derivation and location of zeros of functions. The paper is organised as follows. In sections 2, we recall some facts about splines and quasi-interpolants. In section 3, we describe dQIs of degrees $2\le d\le 5$  and we give their infinity norms and their approximation orders. In sections 4 and 5, we give the associated quadrature formulas and some numerical examples. In section 6, we give the differentiation matrices for quadratic and cubic splines with numerical examples. Finally, in section 7, we show, on a simple polynomial example, how quadratic dQIs can be applied to the location of zeros of functions.

2 Spline spaces on uniform partitions and dQIs

For $I=\left[a,b\right]$  , we denote by ${S}_{d}\left(I,{X}_{n}\right)$  the space of splines of degree $d$  and class ${C}^{d-1}$  on the uniform partition ${X}_{n}=\left\{{x}_{i}=a+ih,0\le i\le n\right\}$  with meshlength $h=\frac{b-a}{n}$  .
A basis of this space is $\left\{{B}_{j},j\in J\right\}$  , with $J=\left\{1,2,\dots ,n+d\right\}$  . With these notations, $supp\left({B}_{j}\right)=\left[{x}_{j-d-1},{x}_{j}\right]$  , and ${\mathcal{N}}_{j}=\left\{{x}_{j-d},\dots ,{x}_{j-1}\right\}$  is the set of the $d$  interior knots in the support of ${B}_{j}$  . As usual, we add multiple knots at the endpoints: $a={x}_{0}={x}_{-1}=\dots ={x}_{-d}$  and $b={x}_{n}={x}_{n+1}=\dots ={x}_{n+d}$  . We recall the represention of monomials in terms of symmetric functions of knots in ${\mathcal{N}}_{j}$  [6,[19.
${e}_{r}\left(x\right)={x}^{r}={\sum }_{j\in J}{\theta }_{j}^{\left(r\right)}{B}_{j}\left(x\right),{\theta }_{j}^{\left(r\right)}={\left(\genfrac{}{}{0}{}{d}{r}\right)}^{-1}sym{m}_{r}\left({\mathcal{N}}_{j}\right),0\le r\le d.$  In particular, the Greville points ${\theta }_{j}={\theta }_{j}^{\left(1\right)}=\frac{1}{d}{\sum }_{\ell =1}^{d}{x}_{j-\ell }$  are the coefficients of ${e}_{1}={\sum }_{j\in J}{\theta }_{j}{B}_{j}$  and the vertices of the control polygon of $S={\sum }_{j\in J}{c}_{j}{B}_{j}$  are the control points $\mathcal{P}=\left\{{\stackrel{~}{c}}_{j}=\left({\theta }_{j},{c}_{j}\right),j\in J\right\}$  . The Schoenberg-Marsden operator is the simplest discrete quasi-interpolant which is exact on the space ${\Pi }_{1}$  :
${S}_{1}f={\sum }_{j\in J}f\left({\theta }_{j}\right){B}_{j},{S}_{1}p=p\forall p\in {\Pi }_{1}.$  —————A discrete quasi-interpolant (abbr. dQI) of degree $d$  is a spline operator of the form:
${Q}_{d}f={\sum }_{j\in J}{\mu }_{j}\left(f\right){B}_{j}$  whose coefficients ${\mu }_{j}\left(f\right)$  are linear combinations of values of $f$  on either the set ${T}_{n}$  ( for $d$  even) or on the set ${X}_{n}$  (for $d$  odd), where ${T}_{n}=\left\{{t}_{j},j\in J\right\},{t}_{j}=\frac{1}{2}\left({x}_{j-2}+{x}_{j-1}\right),{X}_{n}=\left\{{x}_{j},0\le j\le n\right\}.$  Therefore, for $d$  even, we set $f\left({T}_{n}\right)=\left\{{f}_{j}=f\left({t}_{j}\right),j\in J\right\}$  , and for $d$  odd, we set $f\left({X}_{n}\right)=\left\{{f}_{j}=f\left({x}_{j}\right),0\le j\le n\right\}$  . Moreover, ${Q}_{d}$  is exact on ${\Pi }_{d}$  :
${Q}_{d}p=p\forall p\in {\Pi }_{d}.$  For the construction of dQIs, i.e. for the determination of functionals $\left\{{\mu }_{j}\left(f\right),j\in J\right\}$  , the exactness of ${Q}_{d}$  on ${\Pi }_{d}$  amounts to solve a system of linear equations for interior B-splines and a finite number of specific linear systems for boundary B-splines. The determinants of these systems being Vandermonde determinants, there is existence and unicity of dQIs with the above assumptions (see also [3for more general cases). For the sake of completeness, we give below complete formulas for degrees $2\le d\le 5$  . Moreover, we give exact values or upper bounds of $\parallel {Q}_{d}{\parallel }_{\infty }$  and their approximation order on smooth functions. Actually, it is well known (see e.g. [6, chapter 5) that for any subinterval ${I}_{k}=\left[{x}_{k-1},{x}_{k}\right],1\le k\le n$  , and for any function $f$  $\parallel f-{Q}_{d}f{\parallel }_{\infty ,{I}_{k}}\le \left(1+\parallel {Q}_{d}{\parallel }_{\infty }\right){d}_{\infty ,{I}_{k}}\left(f,{\Pi }_{d}\right)$  where the distance of $f$  to polynomials is defined by ${d}_{\infty ,{I}_{k}}\left(f,{\Pi }_{d}\right)=inf\left\{\parallel f-p{\parallel }_{\infty ,{I}_{k}},p\in {\Pi }_{d}\right\}$  Here, as usual, $\parallel f-p{\parallel }_{\infty ,{I}_{k}}=ma{x}_{x\in {I}_{k}}|f\left(x\right)-p\left(x\right)|$  . Therefore, for $f$  smooth enough, e.g. $f\in {C}^{d+1}\left(I\right)$  , this implies that $\parallel f-{Q}_{d}f{\parallel }_{\infty }=O\left({h}^{d+1}\right)$  .

3 Discrete Quasi-Interpolants of degrees $2\le d\le 5$

3.1 ${C}^{1}$  Quadratic dQI

For the ${C}^{1}$  -quadratic dQI ${Q}_{2}f={\sum }_{j=1}^{n+2}{\mu }_{j}\left(f\right){B}_{j}$  , the coefficient functionals are easy to compute (details are given in [15,[16):
${\mu }_{1}\left(f\right)={f}_{1}$  , ${\mu }_{2}\left(f\right)=\frac{1}{6}\left(-2{f}_{1}+9{f}_{2}-{f}_{3}\right)$  , ${\mu }_{n+1}\left(f\right)=\frac{1}{6}\left(-{f}_{n}+9{f}_{n+1}-2{f}_{n+2}\right)$  , ${\mu }_{n+2}\left(f\right)={f}_{n+1}$  , and for $3\le j\le n$  ${\mu }_{j}\left(f\right)=\frac{1}{8}\left(-{f}_{j-1}+10{f}_{j}-{f}_{j+1}\right)$  The exact value $\parallel {Q}_{2}{\parallel }_{\infty }=1.4734$  has been computed in [16. Therefore, for $f\in {C}^{3}\left(I\right)$  for example, we have the following error estimates $\parallel f-{Q}_{2}f{\parallel }_{\infty ,{I}_{k}}\le \frac{5}{2}{d}_{\infty ,{I}_{k}}\left(f,{\Pi }_{2}\right)\text{for}1\le k\le n⟹\parallel f-{Q}_{2}f{\parallel }_{\infty }=O\left({h}^{3}\right).$

3.2 ${C}^{2}$  Cubic dQI

For the ${C}^{2}$  cubic dQI ${Q}_{3}f={\sum }_{j=1}^{n+3}{\mu }_{j}\left(f\right){B}_{j}$  , the coefficient functionals are respectively:
${\mu }_{2}\left(f\right)=\frac{1}{18}\left(7{f}_{0}+18{f}_{1}-9{f}_{2}+2{f}_{3}\right)$  , ${\mu }_{n+2}\left(f\right)=\frac{1}{18}\left(2{f}_{n-3}-9{f}_{n-2}+18{f}_{n-1}+7{f}_{n}\right)$  , ${\mu }_{1}\left(f\right)={f}_{0},{\mu }_{n+3}\left(f\right)={f}_{n}$  , and for $3\le j\le n+1$  ${\mu }_{j}\left(f\right)=\frac{1}{6}\left(-{f}_{j-3}+8{f}_{j-2}-{f}_{j-1}\right).$  As $|{\mu }_{2}{|}_{\infty }=|{\mu }_{n+2}{|}_{\infty }=2$  and $|{\mu }_{j}{|}_{\infty }=\frac{5}{3}$  for $3\le j\le n+1$  , we obtain the upper bound $\parallel {Q}_{3}{\parallel }_{\infty }\le 2$  . It is possible to improve that result by writing the operator in the ”quasi-Lagrange” form:
${Q}_{3}f={\sum }_{j=1}^{n+3}{f}_{j}{L}_{j}$  where the fundamental functions are linear combinations of B-splines, e.g. for $4\le j\le n$  , ${L}_{j}=\frac{1}{6}\left(-{B}_{j+3}+8{B}_{j+2}-{B}_{j+1}\right)$  . It is well known that $\parallel {Q}_{3}{\parallel }_{\infty }$  is equal to the Chebyshev norm of the Lebesgue function:
${\Lambda }_{3}={\sum }_{j=1}^{n+3}|{L}_{j}|.$  In each interval of the uniform partition, ${\Lambda }_{3}$  is bounded above by the cubic polynomial whose Bernstein-Bézier (abbr. BB-) coefficients are sums of absolute values of BB-coefficients of fundamental functions. This allows to see that the maximum of ${\Lambda }_{3}$  is attained in the interval $\left[{x}_{1},{x}_{2}\right]$  and we obtain:
$\parallel {Q}_{3}{\parallel }_{\infty }=\parallel {\Lambda }_{3}{\parallel }_{\infty }\approx 1.631.$  From that we deduce for $f\in {C}^{4}\left(I\right)$  for example, we have the following error estimates $\parallel f-{Q}_{3}f{\parallel }_{\infty ,{I}_{k}}\le \frac{8}{3}{d}_{\infty ,{I}_{k}}\left(f,{\Pi }_{3}\right)\text{for}1\le k\le n⟹\parallel f-{Q}_{3}f{\parallel }_{\infty }=O\left({h}^{4}\right).$

3.3 Quartic dQI

For the ${C}^{3}$  quartic dQI ${Q}_{4}f={\sum }_{j=1}^{n+4}{\mu }_{j}\left(f\right){B}_{j}$  , the coefficient functionals are respectively:
${\mu }_{1}\left(f\right)={f}_{1},{\mu }_{n+4}\left(f\right)={f}_{n+2}$  , ${\mu }_{2}\left(f\right)=\frac{17}{105}{f}_{1}+\frac{35}{32}{f}_{2}-\frac{35}{96}{f}_{3}+\frac{21}{160}{f}_{4}-\frac{5}{224}{f}_{5}$  , ${\mu }_{3}\left(f\right)=-\frac{19}{45}{f}_{1}+\frac{377}{288}{f}_{2}+\frac{61}{288}{f}_{3}-\frac{59}{480}{f}_{4}+\frac{7}{288}{f}_{5}$  , ${\mu }_{4}\left(f\right)=\frac{47}{315}{f}_{1}-\frac{77}{144}{f}_{2}+\frac{251}{144}{f}_{3}-\frac{97}{240}{f}_{4}+\frac{47}{1008}{f}_{5}$  , ${\mu }_{n+1}\left(f\right)=\frac{47}{315}{f}_{n+2}-\frac{77}{144}{f}_{n+1}+\frac{251}{144}{f}_{n}-\frac{97}{240}{f}_{n-1}+\frac{47}{1008}{f}_{n-2}$  , ${\mu }_{n+2}\left(f\right)=-\frac{19}{45}{f}_{n+2}+\frac{377}{288}{f}_{n+1}+\frac{61}{288}{f}_{n}-\frac{59}{480}{f}_{n-1}+\frac{7}{288}{f}_{n-2}$  , ${\mu }_{n+3}\left(f\right)=\frac{17}{105}{f}_{n+2}+\frac{35}{32}{f}_{n+1}-\frac{35}{96}{f}_{n}+\frac{21}{160}{f}_{n-1}-\frac{5}{224}{f}_{n-2}$  , and for $5\le j\le n$  ${\mu }_{j}\left(f\right)=\frac{47}{1152}\left({f}_{j-4}+{f}_{j}\right)-\frac{107}{288}\left({f}_{j-3}+{f}_{j-1}\right)+\frac{319}{192}{f}_{j-2}$  Let us give some details on the computation of functionals ${\mu }_{k},k=2,3,4$  . As ${\mu }_{k}\left({e}_{r}\right)={\theta }_{k}^{\left(r\right)}$  for $0\le r\le 4$  , we determine the five coefficients of the discrete functional ${\mu }_{k}\left(f\right)={\alpha }_{k}{f}_{1}+{\beta }_{k}{f}_{2}+{\gamma }_{k}{f}_{3}+{\delta }_{k}{f}_{4}+{\zeta }_{k}{f}_{5}$  as solutions of the three corresponding linear systems ( $2\le k\le 4$  ) of $5×5$  linear equations ${t}_{1}^{r}{\alpha }_{k}+{t}_{2}^{r}{\beta }_{k}+{t}_{3}^{r}{\gamma }_{k}+{t}_{4}^{r}{\delta }_{k}+{t}_{5}^{r}{\zeta }_{k}={\theta }_{k}^{\left(r\right)},0\le r\le 4$  They have the same Vandermonde determinant ${V}_{5}\left({t}_{1},{t}_{2},{t}_{3},{t}_{4},{t}_{5}\right)\ne 0$  since the ${t}_{i}^{\prime }s$  are distinct. Therefore they have unique solutions. The same technique is applied to the computation of other coefficient functionals.
As $|{\mu }_{2}{|}_{\infty }=|{\mu }_{n+3}{|}_{\infty }\approx 1.77$  , $|{\mu }_{3}{|}_{\infty }=|{\mu }_{n+2}{|}_{\infty }\approx 2.09$  , $|{\mu }_{4}{|}_{\infty }=|{\mu }_{n+1}{|}_{\infty }\approx 2.88$  , $|{\mu }_{j}{|}_{\infty }\approx 2.49$  for $1\le j\le 5$  , we can conclude that $\parallel {Q}_{4}{\parallel }_{\infty }\le 2.88$  and that for $f\in {C}^{5}\left(I\right)$  for example, we have the following error estimates $\parallel f-{Q}_{4}f{\parallel }_{\infty ,{I}_{k}}\le 4{d}_{\infty ,{I}_{k}}\left(f,{\Pi }_{4}\right)\text{for}1\le k\le n⟹\parallel f-{Q}_{4}f{\parallel }_{\infty }=O\left({h}^{5}\right).$

3.4 Quintic dQI

For the ${C}^{4}$  quintic dQI ${Q}_{5}f={\sum }_{j=1}^{n+5}{\mu }_{j}\left(f\right){B}_{j}$  , the coefficient functionals are respectively:
${\mu }_{1}={f}_{0},{\mu }_{n+5}={f}_{n}$  , ${\mu }_{2}=\frac{163}{300}{f}_{0}+{f}_{1}-{f}_{2}+\frac{2}{3}{f}_{3}-\frac{1}{4}{f}_{4}+\frac{1}{25}{f}_{5}$  ${\mu }_{3}=\frac{1}{200}{f}_{0}+\frac{103}{60}{f}_{1}-\frac{73}{60}{f}_{2}+\frac{7}{10}{f}_{3}-\frac{29}{120}{f}_{4}+\frac{11}{300}{f}_{5}$  ${\mu }_{4}=-\frac{41}{400}{f}_{0}+\frac{43}{60}{f}_{1}+\frac{103}{120}{f}_{2}-\frac{7}{10}{f}_{3}+\frac{13}{48}{f}_{4}-\frac{13}{300}{f}_{5}$  ( symmetric formulas for $n+2\le j\le n+4$  ), and for $5\le j\le n$  :
${\mu }_{j}=\frac{13}{240}\left({f}_{j-5}+{f}_{j-1}\right)-\frac{7}{15}\left({f}_{j-4}+{f}_{j-2}\right)+\frac{73}{40}{f}_{j-3}$  As $|{\mu }_{2}{|}_{\infty }=|{\mu }_{n+4}{|}_{\infty }=3.5$  , $|{\mu }_{3}{|}_{\infty }=|{\mu }_{n+3}{|}_{\infty }\approx 3.92$  , $|{\mu }_{4}{|}_{\infty }=|{\mu }_{n+2}{|}_{\infty }\approx 2.69$  , and $|{\mu }_{j}{|}_{\infty }\approx 2.87$  for $5\le j\le n$  , we deduce that $\parallel {Q}_{5}{\parallel }_{\infty }\le 3.92.$  Using a similar technique as for cubics, we find that $\parallel {Q}_{5}{\parallel }_{\infty }\approx 3.106.$  Therefore, for $f\in {C}^{6}\left(I\right)$  for example, we have the following error estimates $\parallel f-{Q}_{5}f{\parallel }_{\infty ,{I}_{k}}\le 4.5{d}_{\infty ,{I}_{k}}\left(f,{\Pi }_{5}\right)\text{for}1\le k\le n⟹\parallel f-{Q}_{5}f{\parallel }_{\infty }=O\left({h}^{6}\right).$

4 Application to numerical integration

Newton-Cotes formulas are obtained by integrating interpolation polynomials (see e.g.
[4,[7,[9). In the same way, integrating spline quasi-interpolants give interesting quadrature formulas (abbr. QF) which are easily deduced from the above computations. We use the notations $\mathcal{ℐ}\left(f\right)={\int }_{a}^{b}f,{\mathcal{ℐ}}_{d}\left(f\right)={\int }_{a}^{b}{Q}_{d}f={\sum }_{j\in J}{\mu }_{j}\left(f\right){\int }_{a}^{b}{B}_{j},{E}_{d}\left(f\right)=\mathcal{ℐ}\left(f\right)-{\mathcal{ℐ}}_{d}\left(f\right).$  As ${\int }_{a}^{b}{B}_{j}=\frac{1}{d}\left({x}_{j}-{x}_{j-d-1}\right)$  and ${\mu }_{j}\left(f\right)$  are known explicitly, we can compute the following quadrature formulas. Moreover, as QIs give the best approximation order, we can conclude that ${E}_{d}\left(f\right)=O\left({h}^{d+1}\right)$  for $f\in {C}^{d+1}\left(I\right)$  , where $h$  is the meshlength. Moreover, as for Newton-Cotes formulas, we get a higher approximation order for even degrees.
(i) QF for quadratics ${\mathcal{ℐ}}_{2}\left(f\right)={\int }_{a}^{b}{Q}_{2}f=h{\sum }_{j=4}^{n-1}{f}_{j}+h\left[\frac{1}{9}\left({f}_{1}+{f}_{n+2}\right)+\frac{7}{8}\left({f}_{2}+{f}_{n+1}\right)+\frac{73}{72}\left({f}_{3}+{f}_{n}\right)\right]$  Error : for $f\in {C}^{4}\left(I\right)$  , ${E}_{2}\left(f\right)=\mathcal{ℐ}\left(f\right)-{\mathcal{ℐ}}_{2}\left(f\right)$  $=\frac{23}{5760}{h}^{4}{D}^{4}f\left(c\right)-\frac{1}{192}{h}^{5}{D}^{4}f\left(\overline{c}\right)$  This result is proved in [16.
Error for Simpson: ${E}_{2}^{*}\left(f\right)=\mathcal{ℐ}\left(f\right)-{\mathcal{ℐ}}_{2}^{*}\left(f\right)=-\frac{1}{180}{h}^{4}{D}^{4}f\left(c\right)$  By comparing the two above errors, we see that the linear combination (extrapolation):
${\stackrel{~}{\mathcal{ℐ}}}_{2}\left(f\right)=\frac{1}{55}\left(32{\mathcal{ℐ}}_{2}\left(f\right)+23{\mathcal{ℐ}}_{2}^{*}\left(f\right)\right)$  is such that $\mathcal{ℐ}\left(f\right)-{\stackrel{~}{\mathcal{ℐ}}}_{2}\left(f\right)=O\left({h}^{5}\right)$  .
(ii) QF for cubics ${\mathcal{ℐ}}_{3}\left(f\right)={\int }_{a}^{b}{Q}_{3}f=h{\sum }_{j=4}^{n-4}{f}_{j}+h\left[\frac{23}{72}\left({f}_{0}+{f}_{n}\right)+\frac{4}{3}\left({f}_{1}+{f}_{n-1}\right)+\frac{19}{24}\left({f}_{2}+{f}_{n-2}\right)+\frac{19}{18}\left({f}_{3}+{f}_{n-3}\right)\right]$  Error: ${E}_{3}\left(f\right)=\mathcal{ℐ}\left(f\right)-{\mathcal{ℐ}}_{3}\left(f\right)=O\left({h}^{4}\right)$  for $f\in {C}^{4}\left(I\right)$  . Numerical experiments show that this formula is not as good as the preceding one.
(iii) QF for quartics ${\int }_{a}^{b}{Q}_{4}f={\mathcal{ℐ}}_{4}\left(f\right)=h{\sum }_{j=6}^{n-3}{f}_{j}+h\left[\frac{206}{1575}\left({f}_{1}+{f}_{n+2}\right)+\frac{107}{128}\left({f}_{2}+{f}_{n+1}\right)+\frac{6019}{5760}\left({f}_{3}+{f}_{n}\right)\right]$  $+h\left[\frac{9467}{9600}\left({f}_{4}+{f}_{n-1}\right)+\frac{13469}{13440}\left({f}_{5}+{f}_{n-2}\right)\right]$  Error: $\mathcal{ℐ}\left(f\right)-{\mathcal{ℐ}}_{4}\left(f\right)=O\left({h}^{6}\right)$  for $f\in {C}^{6}\left(I\right)$  . This is a remarkable formula, which can be compared to the Newton-Cotes formula of the same order. Numerical experiments show that the error for the former QF has also the opposite sign of the error for the latter, as in the quadratic case. The proof will be given elsewhere.
(iv) QF for quintics ${\mathcal{ℐ}}_{5}\left(f\right)={\int }_{a}^{b}{Q}_{5}f=h{\sum }_{j=6}^{n-6}{f}_{j}+h\left[\frac{157}{480}\left({f}_{0}+{f}_{n}\right)+\frac{961}{720}\left({f}_{1}+{f}_{n-1}\right)+\frac{133}{180}\left({f}_{2}+{f}_{n-2}\right)\right]+$  $h\left[\frac{271}{240}\left({f}_{3}+{f}_{n-3}\right)+\frac{1393}{1440}\left({f}_{4}+{f}_{n-4}\right)+\frac{361}{360}\left({f}_{5}+{f}_{n-5}\right)\right]$  Error: ${E}_{5}\left(f\right)={\mathcal{ℐ}}_{\left(f\right)}-{\mathcal{ℐ}}_{5}\left(f\right)=O\left({h}^{6}\right)$  for $f\in {C}^{6}\left(I\right)$  . Numerical experiments show that this formula is not as good as the preceding one.

5 Numerical examples

We compare numerical results on QF applied to the computation of $\mathcal{ℐ}\left({f}_{1}\right)={\int }_{-1}^{1}\frac{1}{1+16{x}^{2}}dx\text{and}\mathcal{ℐ}\left({f}_{2}\right)={\int }_{-1}^{1}{e}^{-x}sin\left(5\pi x\right)dx.$  (i) QF/dQI degrees 2 and 3 ${E}_{2}\left(f\right)=\mathcal{ℐ}\left(f\right)-{\mathcal{ℐ}}_{2}\left(f\right)=O\left({h}^{4}\right)$  , ${E}_{3}\left(f\right)=\mathcal{ℐ}\left(f\right)-{\mathcal{ℐ}}_{3}\left(f\right)=O\left({h}^{4}\right)$  for $f\in {C}^{4}\left(I\right)$  Simpson's QF ${E}_{2}^{*}\left(f\right)=\mathcal{ℐ}\left(f\right)-{\mathcal{ℐ}}_{2}^{*}\left(f\right)=O\left({h}^{4}\right)$  for $f\in {C}^{4}\left(I\right)$  Example1: $\mathcal{ℐ}\left({f}_{1}\right)$
 $\begin{array}{cccc}n& {E}_{2}^{*}& {E}_{2}& {E}_{3}\end{array}$
 $\begin{array}{cccc}& & & \end{array}$
 $\begin{array}{cccc}128& 0.73\left(-9\right)& -0.55\left(-9\right)& -0.44\left(-8\right)\end{array}$
 $\begin{array}{cccc}256& 0;45\left(-10\right)& -0.33\left(-10\right)& -0.26\left(-9\right)\end{array}$
 $\begin{array}{cccc}512& 0.28\left(-11\right)& -0.21\left(-11\right)& -0.15\left(-10\right)\end{array}$
 $\begin{array}{cccc}1024& 0.18\left(-12\right)& -0.13\left(-12\right)& -0.95\left(-12\right)\end{array}$
Example 2: $\mathcal{ℐ}\left({f}_{2}\right)$
 $\begin{array}{cccc}n& {E}_{2}^{*}& {E}_{2}& {E}_{3}\end{array}$
 $\begin{array}{cccc}& & & \end{array}$
 $\begin{array}{cccc}128& 0.14\left(-6\right)& -0.11\left(-6\right)& -0.92\left(-6\right)\end{array}$
 $\begin{array}{cccc}256& 0.90\left(-8\right)& -0.67\left(-8\right)& -0.52\left(-7\right)\end{array}$
 $\begin{array}{cccc}512& 0.56\left(-9\right)& -0.41\left(-9\right)& -0.31\left(-8\right)\end{array}$
 $\begin{array}{cccc}1024& 0.73\left(-9\right)& -0.52\left(-9\right)& -0.37\left(-8\right)\end{array}$
(ii) QF/dQI degree 4 ${E}_{4}\left(f\right)=\mathcal{ℐ}\left(f\right)-{\mathcal{ℐ}}_{4}\left(f\right)=O\left({h}^{6}\right)$  for $f\in {C}^{6}\left(I\right)$  .
Newton-Cotes QF of degree 4: ${E}_{4}^{*}\left(f\right)=\mathcal{ℐ}\left(f\right)-{\mathcal{ℐ}}_{4}^{*}\left(f\right)=O\left({h}^{6}\right)$  for $f\in {C}^{6}\left(I\right)$  .
Example 1: $\mathcal{ℐ}\left({f}_{1}\right)$
 $\begin{array}{ccc}n& {E}_{4}& {E}_{4}^{*}\end{array}$
 $\begin{array}{ccc}& & \end{array}$
 $\begin{array}{ccc}128& -0.83\left(-12\right)& 1.10\left(-12\right)\end{array}$
 $\begin{array}{ccc}256& -0.12\left(-13\right)& 0.24\left(-13\right)\end{array}$
 $\begin{array}{ccc}512& -0.18\left(-15\right)& 0.37\left(-15\right)\end{array}$
 $\begin{array}{ccc}1024& -0.29\left(-17\right)& 0.59\left(-17\right)\end{array}$
Example 2: $\mathcal{ℐ}\left({f}_{2}\right)$
 $\begin{array}{ccc}n& {E}_{4}& {E}_{4}^{*}\end{array}$
 $\begin{array}{ccc}& & \end{array}$
 $\begin{array}{ccc}128& 0.23\left(-7\right)& -0.68\left(-7\right)\end{array}$
 $\begin{array}{ccc}256& 0.44\left(-9\right)& -1.04\left(-9\right)\end{array}$
 $\begin{array}{ccc}512& 0.73\left(-11\right)& -1.62\left(-11\right)\end{array}$
 $\begin{array}{ccc}1024& 0.12\left(-12\right)& -0.25\left(-12\right)\end{array}$
(iii) QF/dQI degree 4: ${E}_{4}\left(f\right)=\mathcal{ℐ}\left(f\right)-{\mathcal{ℐ}}_{4}\left(f\right)=O\left({h}^{6}\right)$  for $f\in {C}^{6}\left(I\right)$  .
QF/dQI degree 5: ${E}_{5}\left(f\right)=\mathcal{ℐ}\left(f\right)-{\mathcal{ℐ}}_{5}\left(f\right)=O\left({h}^{6}\right)$  for $f\in {C}^{6}\left(I\right)$  , Example 1: $\mathcal{ℐ}\left({f}_{1}\right)$
 $\begin{array}{ccc}n& {E}_{4}& {E}_{5}\end{array}$
 $\begin{array}{ccc}& & \end{array}$
 $\begin{array}{ccc}128& -0.83\left(-12\right)& 0.95\left(-11\right)\end{array}$
 $\begin{array}{ccc}256& -0.12\left(-13\right)& 0.14\left(-12\right)\end{array}$
 $\begin{array}{ccc}512& -0.18\left(-15\right)& 0.21\left(-14\right)\end{array}$
 $\begin{array}{ccc}1024& -0.29\left(-17\right)& 0.32\left(-16\right)\end{array}$
Example 2: $\mathcal{ℐ}\left({f}_{2}\right)$
 $\begin{array}{ccc}n& {E}_{4}& {E}_{5}\end{array}$
 $\begin{array}{ccc}& & \end{array}$
 $\begin{array}{ccc}128& 0.23\left(-7\right)& -0.27\left(-6\right)\end{array}$
 $\begin{array}{ccc}256& 0.44\left(-9\right)& -0.50\left(-8\right)\end{array}$
 $\begin{array}{ccc}512& 0.73\left(-11\right)& -0.83\left(-10\right)\end{array}$
 $\begin{array}{ccc}1024& 0.12\left(-12\right)& -0.13\left(-11\right)\end{array}$

6 Application to numerical differentiation

Differentiating interpolation polynomials leads to classical finite differences for the approximate computation of derivatives. Therefore, it seems natural to approximate derivatives of $f$  by derivatives of ${Q}_{d}f$  as long as it is possible, i.e. up to the order $d-1$  . The general theory will be developed elsewhere. Here we only give results for the first derivative and $d=2,3$  . We evaluate $\left({Q}_{d}f{\right)}^{\prime }={\sum }_{j\in J}{\mu }_{j}\left(f\right){B}_{j}^{\prime }$  at points ${T}_{n}$  for $d$  even and at points ${X}_{n}$  for $d$  odd.
(i) Differentiation matrix for quadratics The derivation matrix ${\mathcal{D}}_{2}\in {\mathbb{R}}^{\left(n+2\right)×\left(n+2\right)}$  is defined as follows: setting $y\in {\mathbb{R}}^{n+2}$  for the vector with components ${y}_{j}=f\left({t}_{j}\right),j\in J$  and ${y}^{\prime }\in {\mathbb{R}}^{n+2}$  for the vector with components ${y}_{j}^{\prime }=\left({Q}_{2}f{\right)}^{\prime }\left({t}_{j}\right),j\in J$  , we simply write:
${y}^{\prime }={\mathcal{D}}_{2}y$  ${\mathcal{D}}_{2}=\left(\begin{array}{ccccccccc}-8/3& 3& -1/3& 0& 0& 0& ...& 0& 0\\ -7/6& 11/16& 13/24& -1/16& 0& 0& ...& 0& 0\\ 1/6& -3/4& 1/48& 5/8& -1/16& 0& ...& 0& 0\\ 0& 1/16& -5/8& 0& 5/8& -1/16& ...& 0& 0\\ ...& ...& ...& ...& ...& ...& ...& ...& ...\\ 0& 0& ...& 1/16& -5/8& 0& 5/8& -1/16& 0\\ 0& 0& ...& 0& 1/16& -5/8& -1/48& 3/4& -1/6\\ 0& 0& ...& 0& 0& 1/16& -13/24& -11/16& 7/6\\ 0& 0& ...& 0& 0& 0& 1/3& -3& 8/3\end{array}\right)$  (ii) Differentiation formula for cubics The derivation matrix ${\mathcal{D}}_{3}\in {\mathbb{R}}^{\left(n+1\right)×\left(n+1\right)}$  is defined as follows: setting $y\in {\mathbb{R}}^{n+1}$  for the vector with components ${y}_{j}=f\left({x}_{j}\right),0\le j\le n$  and ${y}^{\prime }\in {\mathbb{R}}^{n+1}$  for the vector with components ${y}_{j}^{\prime }=\left({Q}_{3}f{\right)}^{\prime }\left({x}_{j}\right),0\le j\le n$  , we obtain:
${y}^{\prime }={\mathcal{D}}_{3}y$  ${\mathcal{D}}_{3}=\left(\begin{array}{ccccccccc}-11/6& 3& -3/2& 1/3& 0& 0& ...& 0& 0\\ -1/3& -1/2& 1& -1/6& 0& 0& ...& 0& 0\\ 1/12& -2/3& 0& 2/3& -1/12& 0& ...& 0& 0\\ 0& 1/12& -2/3& 0& 2/3& -1/12& ...& 0& 0\\ ...& ...& ...& ...& ...& ...& ...& ...& ...\\ 0& 0& ...& 1/12& -2/3& 0& 2/3& -1/12& 0\\ 0& 0& ...& 0& 1/12& -2/3& 0& 2/3& -1/12\\ 0& 0& ...& 0& 0& 1/6& -1& 1/2& 1/3\\ 0& 0& ...& 0& 0& -1/3& 3/2& -3& 11/6\end{array}\right)$  (iii) Some numerical results Again we use the two functions ${f}_{1}\left(x\right)=\frac{1}{1+16{x}^{2}}$  and ${f}_{2}\left(x\right)={e}^{-x}sin\left(5\pi x\right)$  on the interval $I=\left[-1,1\right]$  . For $p=1,2$  , we set ${\varepsilon }_{p}={max}_{v\in V}|{f}_{p}^{\prime }\left(v\right)-\left({Q}_{d}{f}_{p}{\right)}^{\prime }\left(v\right)|$  where $V={T}_{n}$  (resp. $V={X}_{n}$  ) for $d=2$  (resp. $d=3$  ) and ${\varepsilon }_{p}^{*}={max}_{v\in V}|{f}_{p}^{\prime }\left(v\right)-\delta {f}_{p}\left(v\right)|$  where $\delta {f}_{p}\left(v\right)$  is the classical centered approximation of ${f}_{p}^{\prime }\left(v\right)$  of order $2$  (with standard modifications at the endpoints).
For quadratics, we obtain the following results
 $\begin{array}{ccccc}n& {\varepsilon }_{1}& {\varepsilon }_{1}^{*}& {\varepsilon }_{2}& {\varepsilon }_{2}^{*}\end{array}$
 $\begin{array}{ccc}& & \end{array}$
 $\begin{array}{ccccc}64& 0.014009& 0.047853& 0.016143& 0.046317\end{array}$
 $\begin{array}{ccccc}128& 0.003138& 0.012079& 0.003674& 0.011606\end{array}$
 $\begin{array}{ccccc}256& 0.000767& 0.003036& 0.000872& 0.002902\end{array}$
 $\begin{array}{ccccc}512& 0.000190& 0.000759& 0.000212& 0.00725\end{array}$
 $\begin{array}{ccccc}1024& 0.0000475& 0.0001899& 0.000052& 0.000181\end{array}$
We see that the orders are all $O\left({h}^{2}\right)$  . However, the errors for the derivatives of the quadratic QI ( ${\varepsilon }_{1}$  and ${\varepsilon }_{2}$  ) are between $3$  and $4$  times less than the errors for the centered finite differences ( ${\varepsilon }_{1}^{*}$  and ${\varepsilon }_{2}^{*}$  ).
For cubics, we obtain the following results
 $\begin{array}{ccccc}n& {\varepsilon }_{1}& {\varepsilon }_{1}^{*}& {\varepsilon }_{2}& {\varepsilon }_{2}^{*}\end{array}$
 $\begin{array}{ccc}& & \end{array}$
 $\begin{array}{ccccc}64& 3.0\left(-3\right)& 4.7\left(-2\right)& 1.0\left(-2\right)& 4.7\left(-2\right)\end{array}$
 $\begin{array}{ccccc}128& 2.0\left(-4\right)& 1.2\left(-2\right)& 1.4\left(-3\right)& 1.2\left(-2\right)\end{array}$
 $\begin{array}{ccccc}256& 1.3\left(-5\right)& 3.0\left(-3\right)& 1.8\left(-4\right)& 2.9\left(-3\right)\end{array}$
 $\begin{array}{ccccc}512& 8.0\left(-7\right)& 7.6\left(-4\right)& 2.4\left(-5\right)& 7.2\left(-4\right)\end{array}$
 $\begin{array}{ccccc}1024& 5.0\left(-8\right)& 1.9\left(-4\right)& 3.0\left(-6\right)& 1.8E\left(-4\right)\end{array}$
Of course, ${\varepsilon }_{1}^{*}$  and ${\varepsilon }_{2}^{*}$  are both $O\left({h}^{2}\right)$  and ${\varepsilon }_{1}$  and ${\varepsilon }_{2}$  are both at least $O\left({h}^{3}\right)$  . However, for the function ${f}_{1}$  , a superconvergence phenomenon occurs because we have ${\varepsilon }_{1}=O\left({h}^{4}\right)$  instead of $O\left({h}^{3}\right)$  . We shall study this kind of results in a further paper.
————–

7 Approximating zeros of a function by those of a quadratic dQI

Let $f$  be a continuous function defined on $I=\left[a,b\right]$  . In order to locate the zeros of $f$  in this interval, we approximate $f$  by its ${C}^{1}$  quadratic dQI $g={Q}_{2}f$  and we compute the exact zeros of $g$  : this is quite possible because $g$  is piecewise quadratic. The complete study will be done elsewhere. Here we take a simple example: we want to approximate the zeros of the Legendre polynomial ${P}_{8}\left(x\right)=$  in the interval $I=\left[-1,1\right]$  . The five zeros of ${P}_{8}$  are respectively $\left\{±{x}_{1},±{x}_{2},±{x}_{3},±{x}_{4}\right\}$  , with ${x}_{1}=.1834346425,{x}_{2}=.5255324099,{x}_{3}=.7966664774,{x}_{4}=.9602898565.$  The following array gives the errors ${\varepsilon }_{k}={x}_{k}-{\overline{x}}_{k},1\le k\le 4$  where ${\overline{x}}_{k}$  is the zero of $g$  nearest to ${x}_{k}$  .
 $\begin{array}{ccccc}n& {\varepsilon }_{1}& {\varepsilon }_{2}& {\varepsilon }_{3}& {\varepsilon }_{4}\end{array}$
 $\begin{array}{ccccc}& & & & \end{array}$
 $\begin{array}{ccccc}16& .000543& .003784& .013753& -.007841\end{array}$
 $\begin{array}{ccccc}32& -.000043& .000210& .000556& -.001017\end{array}$
 $\begin{array}{ccccc}64& -.000013& -.000012& .000043& .000026\end{array}$
Acknowledgements: the author thanks Professor Catterina Dagnino and the Department of Mathematics of the University of Turin for their kind invitation to deliver this seminar during his stay from January 12 to 20, 2005.
References

1. C.de Boor, A practical guide to splines, Revised edition. Springer-Verlag, New-York (2001).
2. C.de Boor, Splines as linear combinations of B-splines, Approximation Theory II, G.G. Lorentz et al. (eds), Academic Press, New-York (1976), 1-47.
3. G.Chen, C.K. Chui, M.J. Lai, Construction of real-time spline quasi-interpolation schemes, Approx. Theory Appl. 4 (1988), 61-75.
4. P.J. Davis, P. Rabonowitz, Numerical integration. Blaisdell, Waltham (1967).
5. S.A. De Swardt, J.M. De Villiers, Gregory type quadrature based on quadratic nodal spline interpolation. Numer. Math. 85 (2000), 129-153.
6. R.A. DeVore, G.G. Lorentz, Constructive approximation, Springer-Verlag, Berlin (1993)
8. W. Gautschi: Orthogonal polynomials: applications and computation. Dans Acta Numerica 1996, A. Iserles (ed.), CUP 1996.
9. A. Krommer, Ch. W. Ueberhuber, Computational integration. SIAM, Philadelphia (1998).
10. V. Lampret, An invitation to Hermite's integration and summation: a comparison between Hermite's and Simpson's rules. SIAM Review 46, No 2 (2004), 329-345.
11. B.G. Lee, T. Lyche, L.L. Schumaker, Some examples of quasi-interpolants constructed from local spline projectors. In Math methods for CAGD Oslo II, 243-252.
12. T. Lyche, L.L. Schumaker, Local spline approximation, J. Approx. Theory 15 (1975), 294-325.
13. M.J.D. Powell, Approximation theory and methods , Cambridge University Press (1981).
14. P. Sablonnière: Quasi-interpolantes splines sobre particiones uniformes. First Meeting in Approximation Theory of the University of Jaén (Ubeda, June 29-July 2, 2000). Prépublication IRMAR 00-38, Rennes (June 2000).
15. P. Sablonnière: On some multivariate quadratic spline quasi-interpolants on bounded domains. In Modern Developments in Multivariate Approximation, W. Haussmann et al. (eds), ISNM Vol. 145, Birkhäuser Verlag (2003), 263-278.
16. P. Sablonnière: Quadratic spline quasi-interpolants on bounded domains of ${\mathbb{R}}^{d},d=1,2,3$  . In Spline and radial functions, Rend. Sem. Univ. Pol. Torino, Vol. 61 (2003), 61-78.
17. P. Sablonnière: A quadrature formula associated with a univariate quadratic spline quasi-interpolant. Prépublication IRMAR, Rennes, April 2005 (submitted).
18. P. Sablonnière, D. Sbibih: Integral spline operators exact on polynomials. Approx. Theory Appl. 10, No 3 (1994), 56-73.
19. L.L. Schumaker, Spline functions: basic theory, John Wiley and Sons, New-York (1981).