Оценка остаточного члена одной кубатурной формулы по чебышевской сетке

Обозначим через $C(Q)$ пространство непрерывных в квадрате $Q=[-1,1]\times[-1,1]$ функций $f(x,y)$ с нормой
\begin{equation} \| f\|=\max(|(f(x,y)|), \quad (x,y)\in Q \end{equation}
Рассматривается кубатурная формула по чебышевской сетке вида
\begin{eqnarray} &\int_{-1}^1\int_{-1}^1\frac{1}{\sqrt{(1-x^2)(1-y^2)}}f(x,y)dxdy= \\ &\frac{\pi^2}{mn}\sum_{i=1}^n\sum_{j=1}^mf\big (\cos\frac{2i-1}{2n}\pi,cos\frac{2j-1}{2m}\pi\big )+R_{m,n}(f) \end{eqnarray}
в некотором классе $H(r_1,r_2)$ функцией $f\in C(Q)$, определяемого при помощи оператора обобщенного сдвига. Для остаточного члена $R_{m,n}(f)$ доказана следующая оценка:
$$ \sup_{f\in H(r_1,r_2)}| R_{m,n}(f) |=O(n^{-r_1+1}+m^{-r_2+1}) $$
где $r_1,r_2>1,\lambda^{-1}\leq n/m\leq\lambda,\lambda>0$; константа, входящая в $O(1)$, зависит от $\lambda$. Библ. 4. Ключевые слова: кубатурная формула, чебышевская сетка, оценка остаточного члена.