Iterative Shepard operator of least squares thin-plate spline type

We propose a modified bivariate Shepard operator of least squares thin-plate spline type based on an iterative method, introduced by A.V. Masjukov and V.V. Masjukov. This iterative method does not require an artificial setup of the parameters and it is based on successive scaling. For the numerical examples, using an idea of J. R. McMahon, we construct some representative sets of knot points for the initial sets of the interpolation nodes.