“…Piecewise interpolation schemes may be classified as local (see, for example, [1][2][3][4][5]), which may be calculated in any subinterval independently of other subintervals, and global, which calculation involves solving of the system of I or I*k equations, I -number of interpolation knots. Recent spline belongs to the second class, which contains, for example, rational splines [6,7], discrete hyperbolic tension splines [8,9]. It is known that high accuracy interpolations are connected with possible producing of undesirable oscillations near points of interpolated function discontinuities or points of discontinuities of the function first derivative.…”