Formulas for Bivariate Hyperosculatory Interpolation

For a given function f(x, y), bivariate hyperosculatory interpolation formulas are obtained by employing a suitably constructed binary nic pn(x, y) that is fitted to the values of f(x, y) and its first and second partial derivatives at the m points (xi, y;) of