“…In a nutshell, a Shepard approximation scheme employs rational formations of shifts of an appropriately-selected base function to approximate a target function, and its efficiency epitomizes in the reproduction of constants. Because of the singularity of the base function at zero, the original Shepard approximant interpolates values of the target function at x j , 1 ≤ j ≤ n. That is, S Φ,n (x j ) = f (x j ), 1 ≤ j ≤ n. However, for most other choices of base functions, the interpolation feature is lost, and the resulted approximants are called "quasi-interpolants" in the literature; see [4], [17], [18]. Except for the cases in which the base function is compactly supported (see [16], [17], [18]), deriving optimal order error estimates for a Shepard approximation scheme has been uncommon.…”