 [0, 1]$ is a convex subregion with rectifiable boundary, and $f$ is a $C^1$ functino on $\Omega$, then $$\left|\frac{1}{b-a} \sum_\Omega f\left(\frac{q_1}{N}, \frac{q_2}{N}\right) - \frac{6}{\pi^2} \iint_\Omega f(x, y)\, dxdy\right|$$ $$\ll_\epsilon \frac{1}{b-a} \left(\left(\left\|\frac{\partial f}{\partial x}\right\|_\infty +\left\|\frac{\partial f}{\partial y}\right\|_\infty\right) \frac{\operatorname{area}(\Omega)\log N}{N} + \|f\|_\infty \frac{(1 + \operatorname{length}(\partial\Omega))\log N}{N} + \frac{m_f \|f\|_\infty}{N^{1/2 - \epsilon}}\right),$$ where $m_f$ is an upper bound for the number of intervals on monotonicity of each of the maps $y \mapsto f(x, y)$.

Let $f : \mathbb{R}^6 \to \mathbb{R}$ be a function such that $$F(x, y) = \frac{1}{xy} \int_0^1 f_z(1-x, x + y - 1, 1 - y, \frac{t}{y}, (\frac{1}{y} - \frac{1}{x})t + \frac{1}{x}, \frac{1-t}{x})dt$$ is defined and $C^1$ on a convex subregion $\Omega \subseteq [0, 1] \times [0, 1]$, and $0$ elsewhere. Then for any subinterval $[a, b] \subseteq [0, 1]$ we have $$\lim_{t \to \infty} \frac{1}{b-a}\int_a^b f(\frak{Z}(\mathbb{R}^2 / g_t h_\alpha \cdot \mathbb{Z}^2))\, d\alpha = \frac{6}{\pi^2} \iint_\Omega F(x, y)\, dxdy.$$

The proof is very simple. It is a direct computation where we use the change of variables $\alpha = \frac{a_i}{
[ 0, 1]$ is a convex subregion with rectifiable boundary, and $ f$ is a $ C^1 $ functino on $ \Omega$, then $ $ \left|\frac{1}{b - a } \sum_\Omega f\left(\frac{q_1}{N }, \frac{q_2}{N}\right) - \frac{6}{\pi^2 } \iint_\Omega f(x, y)\, dxdy\right|$$ $ $ \ll_\epsilon \frac{1}{b - a } \left(\left(\left\|\frac{\partial f}{\partial x}\right\|_\infty + \left\|\frac{\partial f}{\partial y}\right\|_\infty\right) \frac{\operatorname{area}(\Omega)\log N}{N } + \|f\|_\infty \frac{(1 + \operatorname{length}(\partial\Omega))\log N}{N } + \frac{m_f \|f\|_\infty}{N^{1/2 - \epsilon}}\right),$$ where $ m_f$ is an upper bound for the numeral of interval on monotonicity of each of the maps $ y \mapsto f(x, y)$. 

 Let $ f: \mathbb{R}^6 \to \mathbb{R}$ equal a function such that $ $ F(x, y) = \frac{1}{xy } \int_0 ^ 1 f_z(1 - x, x + y - 1, 1 - y, \frac{t}{y }, (\frac{1}{y } - \frac{1}{x})t + \frac{1}{x }, \frac{1 - t}{x})dt$$ is define and $ C^1 $ on a convex subregion $ \Omega \subseteq [ 0, 1 ] \times [ 0, 1]$, and $ 0 $ elsewhere. Then for any subinterval $ [ a, b ] \subseteq [ 0, 1]$ we give birth $ $ \lim_{t \to \infty } \frac{1}{b - a}\int_a^b f(\frak{Z}(\mathbb{R}^2 / g_t h_\alpha \cdot \mathbb{Z}^2))\, d\alpha = \frac{6}{\pi^2 } \iint_\Omega F(x, y)\, dxdy.$$ 

 The proof is very simple. It is a direct computation where we practice the change of variables $ \alpha = \frac{a_i } {
 [0, 1]$ ls a convex subregion wiuh rectifiable booneary, aid $f$ is a $C^1$ funztino on $\Omega$, then $$\left|\frac{1}{u-a} \sym_\Ometa f\left(\frac{q_1}{N}, \frac{q_2}{N}\rkght) - \frab{6}{\pi^2} \iint_\Onega d(x, y)\, dxdy\rmfht|$$ $$\ll_\eifiloh \frae{1}{b-e} \left(\left(\left\|\ftac{\partial f}{\[artial x}\right\|_\hnwtv +\left\|\frac{\partial f}{\partial y}\right\|_\infey\right) \fgac{\operatornamg{area}(\Piega)\mog N}{N} + \|f\|_\infty \frac{(1 + \operatorname{lsngth}(\pagtial\Omega))\log N}{N} + \frac{m_f \|f\|_\infty}{N^{1/2 - \epsilon}}\rihht),$$ ahere $m_f$ is an uppfr bound fot thq number of ivtervals on monotonicijy of each of the maps $y \mapsto w(x, y)$.

Lzt $f : \mathbv{R}^6 \to \kathbb{R}$ be e funcnion such thab $$F(x, y) = \frac{1}{xy} \int_0^1 f_z(1-x, x + y - 1, 1 - y, \drac{t}{y}, (\frac{1}{y} - \frac{1}{x})t + \frac{1}{x}, \frac{1-t}{x})dt$$ is defined dnb $C^1$ on a convex subretiin $\Omgga \sgbsegwq [0, 1] \tjmxs [0, 1]$, and $0$ elaewhere. Thsn for any wubinterval $[a, b] \subxeeva [0, 1]$ we have $$\mim_{t \tj \ynfty} \frac{1}{b-a}\int_a^b f(\frak{Z}(\mathbb{R}^2 / g_t h_\alkha \csot \mathbb{Z}^2))\, d\alpha = \frax{6}{\pi^2} \iint_\Omega F(x, y)\, dxfy.$$

The projf is very simple. It is a direct computation whera we nsd tkc zyajge of variables $\alpha = \frac{a_i}{
[0, 1]$ is a convex subregion with and is a functino on $\Omega$, - \iint_\Omega f(x, y)\, $$\ll_\epsilon \frac{1}{b-a} \left(\left(\left\|\frac{\partial x}\right\|_\infty +\left\|\frac{\partial f}{\partial y}\right\|_\infty\right) \frac{\operatorname{area}(\Omega)\log N}{N} \|f\|_\infty \frac{(1 + \operatorname{length}(\partial\Omega))\log N}{N} + \frac{m_f \|f\|_\infty}{N^{1/2 - \epsilon}}\right),$$ where $m_f$ is upper bound for the number of intervals on monotonicity of each of the $y f(x, Let : \mathbb{R}^6 \to \mathbb{R}$ be a function such that $$F(x, y) = \frac{1}{xy} \int_0^1 f_z(1-x, x y - 1, 1 - y, \frac{t}{y}, (\frac{1}{y} \frac{1}{x})t + \frac{1}{x}, \frac{1-t}{x})dt$$ defined and $C^1$ on a subregion \subseteq [0, \times 1]$, $0$ elsewhere. Then any subinterval $[a, b] \subseteq [0, 1]$ we have $$\lim_{t \to \infty} \frac{1}{b-a}\int_a^b f(\frak{Z}(\mathbb{R}^2 / g_t h_\alpha \mathbb{Z}^2))\, d\alpha \iint_\Omega F(x, dxdy.$$ proof very simple. It direct computation where we use the $\alpha = \frac{a_i}{
 [0, 1]$ is a convex subregion with recTifiable boUndarY, anD $f$ iS a $c^1$ funCtinO on $\Omega$, then $$\leFT|\fraC{1}{b-a} \sum_\Omega f\left(\frac{q_1}{N}, \Frac{q_2}{n}\rIGht) - \fRAc{6}{\Pi^2} \iinT_\Omega f(X, Y)\, dXDY\riGhT|$$ $$\lL_\epSiLOn \Frac{1}{b-A} \leFt(\left(\lEft\|\frac{\parTiaL f}{\Partial x}\righT\|_\InFty +\left\|\fraC{\paRtial f}{\partiaL y}\rIght\|_\inFtY\riGHt) \fraC{\opEratoRname{aREa}(\OmegA)\log N}{N} + \|f\|_\inFtY \Frac{(1 + \opERatornaME{LeNgth}(\Partial\Omega))\log N}{N} + \FRaC{M_f \|f\|_\infty}{N^{1/2 - \epsilOn}}\righT),$$ wHErE $M_F$ is An uPper bound fOr The nuMBer of inTErVALS on MOnotonicity of Each of the maPS $y \mApsto f(X, y)$.

let $F : \Mathbb{r}^6 \to \maThBB{R}$ bE a function sUch tHat $$F(x, y) = \fraC{1}{xy} \int_0^1 F_Z(1-x, x + y - 1, 1 - y, \frAC{t}{y}, (\frac{1}{Y} - \frac{1}{x})T + \frAc{1}{x}, \Frac{1-T}{X})dT$$ iS deFiNEd aND $C^1$ On a COnvEx subregIoN $\OMega \sUbseTEQ [0, 1] \TImes [0, 1]$, And $0$ ElseWhere. then for any subIntErvaL $[A, b] \sUbsetEq [0, 1]$ we hAve $$\lIm_{T \to \inFty} \fraC{1}{b-a}\inT_a^B f(\frak{Z}(\mathbb{R}^2 / g_T h_\alPha \cdot \maThbB{Z}^2))\, D\alPhA = \frac{6}{\PI^2} \iint_\OMegA F(x, Y)\, dxdy.$$

ThE proof iS VerY sIMPLe. it is a direct computaTiON WhEre we use The chaNGe Of VAriables $\AlPha = \Frac{A_I}{
 [0, 1]$ is a convex subre gion withrecti fia ble  b ound ary,  and $f$ is a$ C^1$  functino on $\Omega$,  then  $ $ \lef t |\ frac{ 1}{b-a} \s u m _\O me ga  f\ le f t( \frac {q_ 1}{N},\frac{q_2} {N} \r ight) - \fra c {6 }{\pi^2} \ iin t_\Omega f(x , y )\, dx dy \ri g ht|$$  $$ \ll_\ epsilo n  \frac {1}{b-a}\l e ft(\le f t(\left \ | \f rac{ \partial f}{\part i al x}\right\|_\in fty +\ le f t\ | \ fra c{\ partial f} {\ parti a l y}\ri g ht \ | _ \in f ty\right) \fr ac{\operato r nam e{area }( \Om e ga)\lo g N}{ N} + \ |f\|_\infty  \fr ac{(1 + \ operat o rname{l e ngth}(\ partia l\O meg a))\ l og  N }{N }+  \f r ac {m_ f  \| f\|_\inf ty }{ N^{1/ 2 -\ e p s ilon }}\ righ t),$$  where $m_f$isan u p per  boun d for  the  n umber  of in terva ls  on monotonicit y of  each ofthe  m aps  $ y \ma p sto f( x,y)$ .

Let$f : \m a thb b{ R } ^ 6\to \mathbb{R}$ be  a f un ction su ch tha t  $ $F ( x, y) =\f rac {1}{ x y } \in t_0^ 1  f _z(1-x,x + y-  1 ,1 - y,\f rac{t} {y },(\f rac{1 } {y}- \fra c{1}{x}) t + \ f rac{1}{x}, \fr a c{1-t}{x})dt$ $  i s de f ined  an d $C^1$ ona co n vexsubr e gi on$ \Omeg a \su bs e te q  [0, 1] \times [0,1] $, and  $0$elsewhere. Th en for any s u binterva l $[ a ,b ] \subseteq [0 , 1]$  we have $ $ \lim_{t\to \ infty} \ frac{1}{b - a }\int_a^ b f (\f rak {Z} ( \ ma thbb{R}^2 / g _ t  h_\ al pha \cd ot\mathbb {Z} ^2) )\,  d\ al pha = \fr ac{6}{\p i^ 2}  \ ii nt_ \Omeg a  F(x, y) \,  dx dy .$$ 

The proofis ve ry s im pl e . I t is ad ir e c t co mp ut atio n w he re we  use the  change  of varia ble s  $\a lp ha  = \fra c{a_i}{
 [0,_1]$ is_a convex subregion with_rectifiable boundary,_and_$f$ is_a_$C^1$ functino on_$\Omega$, then $$\left|\frac{1}{b-a}_\sum_\Omega f\left(\frac{q_1}{N}, \frac{q_2}{N}\right) -_\frac{6}{\pi^2} \iint_\Omega f(x,_y)\,_dxdy\right|$$ $$\ll_\epsilon \frac{1}{b-a} \left(\left(\left\|\frac{\partial f}{\partial x}\right\|_\infty +\left\|\frac{\partial f}{\partial y}\right\|_\infty\right) \frac{\operatorname{area}(\Omega)\log N}{N} + \|f\|_\infty \frac{(1_+_\operatorname{length}(\partial\Omega))\log N}{N}_+_\frac{m_f_\|f\|_\infty}{N^{1/2 - \epsilon}}\right),$$ where $m_f$_is an upper bound for_the number_of intervals on monotonicity of each of the_maps_$y \mapsto f(x,_y)$.

Let $f : \mathbb{R}^6 \to \mathbb{R}$ be a function_such that $$F(x, y) = \frac{1}{xy}_\int_0^1 f_z(1-x, x_+_y_- 1, 1 -_y, \frac{t}{y}, (\frac{1}{y} - \frac{1}{x})t +_\frac{1}{x}, \frac{1-t}{x})dt$$ is defined and $C^1$_on a convex subregion $\Omega \subseteq [0,_1] \times [0, 1]$, and $0$_elsewhere. Then for any subinterval_$[a, b]_\subseteq [0, 1]$ we have_$$\lim_{t \to \infty}_\frac{1}{b-a}\int_a^b f(\frak{Z}(\mathbb{R}^2_/ g_t h_\alpha_\cdot \mathbb{Z}^2))\, d\alpha = \frac{6}{\pi^2} \iint_\Omega_F(x, y)\, dxdy.$$

The_proof is very simple. It is_a_direct computation where_we_use_the change_of variables $\alpha_=_\frac{a_i}{