 Solution

Use the standard knot-domain convention in which the spline space contains every polynomial of degree at most $k-1$ on $[a,b]$, with completed or clamped endpoint knots. The De Boor–Fix functional extracts a local spline coefficient, and its Hahn-Banach theorem extension is linear on continuous local data with the stated bound. Only that bound and spline reproduction are needed here.

For $x$ in a knot cell $[t_j,t_{j+1}]$, support and Gram bandwidth of B-splines implies that only indicescan contribute, restricted also to $1\le i\le n$. Every corresponding local support $[t_i,t_{i+k}]$ lies insideNonnegativity and subpartition of unity for B-splines therefore giveTaking the supremum over the cell proves the local stability estimate for spline quasi-interpolation:Endpoint limits follow by continuity for $k\ge 2$; the order-one version uses the usual one-sided interval convention.

Let $h={\max }_i|t_{i+1}-t_i|$ and $f\in C^k[a,b]$. Choose any point $y_j$ in the stencil and form its degree-$k-1$ Taylor polynomial $p_j$. The Taylor theorem remainder givesThe polynomial is itself a spline, so $Q{p}_j=p_j$. By linearity and the local stability estimate,Taking the maximum over cells yieldsThe constant depends only on $k$ and the supplied functional bound, not on mesh ratios or the number of knots.

The endpoint convention is essential to the global assertion. The domain can be the basic knot interval, with any needed exterior data extended smoothly, or the full approximation interval can use clamped endpoint knots. If one instead takes a finite, unclamped collection on the entire outer support interval $[t_1,t_{n+k}]$, every order-$k\ge 2$ basis function vanishes at $t_1$, so $Q1(t_1)=0$ and uniform approximation to the constant one there is impossible. This does not affect the local estimate; it explains the usual boundary completion implicit in the global spline approximation statement.