On some interpolation process by polynomials of degree 4m+1 on the triangle (Q4254099)

From the finite element method, it is known that the error of polynomial interpolation over an triangle \(\Delta\) depends on geometric properties of \(\Delta\). Let \(f\) be a smooth function of the class \(W^{4m+ 2}M\) \((m\geq 1)\), i.e., all directional derivatives of \(f\) of order \(\leq 4m+ 2\) are continuous in \(\Delta\) and bounded by \(M\). Let \(P_{4m+ 1}(x,y)\) denote a special interpolation polynomial of degree \(\leq 4m+1\) with respect to \(f\in W^{4m+ 2}M\). In this paper, \[ \Biggl|{\partial^n\over\partial x^{n- j}\partial y^j} (f(x, y)- P_{4m+ 1}(x,y))\Biggr|\qquad (0\leq n\leq 4m+ 1) \] are estimated in terms of the longest side and the angles of \(\Delta\). Some of these results are unimprovable in \(W^{4m+ 2}M\). Note that the case \(m=1\) was studied by \textit{Yu. N. Subbotin} [Trudy Inst. Math. Mech. (Ekaterinburg) 2, 110-119 (1992; Zbl 0812.41002)].