Given a B-spline M on R, s≥ 1 we consider a classical discrete quasi-interpolant Qd written in the form Qdf = ∑ i∈Zs f (i)L(· − i), where L(x) :=∑ j∈JcjM(x− j) for some finite subset J ⊂ Z and cj ∈R. This fundamental function is determined to produce a quasi-interpolation operator exact on the space of polynomials of maximal total degree included in the space spanned by the integer translates o...