2}^n \sum_{j=1}^{k-1}  
          \l_{j,k} f(u_{j,k}, v_{j,k}), \quad f \in \Pi_{2n-3}^2,$$ where $\lambda_{j,k} = \l_j \l_k (x_{j,n} - x_{k,n})^2$.

The product of is a cubature rule on $[-1,1]^2$ $$\label{ProdGauss}
    c_w^2 \int_{[-1,1]^2} f(x,y) w(x)w(y) dx dy = \sum_{k=1}^n \sum_{j=1}^n 
          \l_k \l_j f(x_{k,n}, x_{j,n}),$$ which is exact for $f \in \Pi_{2n-1} \times \Pi_{2n-1}$, the space of polynomials of degree at most $2n-1$ in either $x$ or $y$ variable. Applying on the symmetric polynomials $f(x+y, xy)$ and using the symmetry, we obtain $$c_w^2 \int_{\triangle} f(x+y, xy) w(x)w(y) dx dy = \sum_{k=1}^n \mathop{ {\sum}' }_{j=1}^k  
          \l_k \l_j f(x_{k,n}+ x_{j,n}, x_{k,n}x_{j,n}),$$ exact for polynomials $f$ in $\Pi_{2n-1}\times \Pi_{2n-1}$. Under the change of variables $u = x+y$ and $v = xy$ and by, the above cubature becomes, since $\Pi_{2n-1}\times \Pi_{2n-1}$ becomes $\Pi_{2n-1}^2$ under the mapping $(x,y) \mapsto (u,v)$. It is easy to see that has $\dim \Pi_{n-1}^2$ nodes, so that it is a Gaussian cubature rule.

To prove, we apply the product Gaussian cubature rule on functions of the form $(x - y)^2 f(x+y,xy)$ for $f \
2}^n \sum_{j=1}^{k-1 }  
           \l_{j, k } f(u_{j, k }, v_{j, k }), \quad f \in \Pi_{2n-3}^2,$$ where $ \lambda_{j, k } = \l_j \l_k (x_{j, n } - x_{k, n})^2$. 

 The product of is a cubature rule on $ [ -1,1]^2 $ $ $ \label{ProdGauss } 
     c_w^2 \int_{[-1,1]^2 } f(x, y) w(x)w(y) dx dy = \sum_{k=1}^n \sum_{j=1}^n 
           \l_k \l_j f(x_{k, n }, x_{j, n}),$$ which is accurate for $ f \in \Pi_{2n-1 } \times \Pi_{2n-1}$, the quad of polynomials of degree at about $ 2n-1 $ in either $ x$ or $ y$ variable. enforce on the symmetric polynomials $ f(x+y, xy)$ and practice the symmetry, we obtain $ $ c_w^2 \int_{\triangle } f(x+y, xy) w(x)w(y) dx dy = \sum_{k=1}^n \mathop { { \sum }' } _ { j=1}^k  
           \l_k \l_j f(x_{k, n}+ x_{j, n }, x_{k, n}x_{j, n}),$$ exact for polynomial $ f$ in $ \Pi_{2n-1}\times \Pi_{2n-1}$. Under the variety of variables $ u = x+y$ and $ five = xy$ and by, the above cubature becomes, since $ \Pi_{2n-1}\times \Pi_{2n-1}$ become $ \Pi_{2n-1}^2 $ under the mapping $ (x, y) \mapsto (u, v)$. It is easy to see that has $ \dim \Pi_{n-1}^2 $ node, so that it is a Gaussian cubature rule. 

 To prove, we apply the intersection Gaussian cubature rule on functions of the form $ (x - y)^2 f(x+y, xy)$ for $ f \
2}^n \skm_{j=1}^{k-1}  
          \l_{j,k} f(u_{j,k}, v_{j,k}), \quad f \in \Ki_{2b-3}^2,$$ wherx $\lambdz_{j,k} = \l_j \u_k (x_{j,n} - x_{k,n})^2$.

The product of is a cubaulre rule on $[-1,1]^2$ $$\label{ProaGauss}
    c_w^2 \int_{[-1,1]^2} d(x,y) x(x)w(y) dx dy = \sum_{k=1}^i \sum_{j=1}^n 
          \l_n \l_j f(x_{k,n}, x_{j,n}),$$ wmich is exawt for $f \in \Pi_{2t-1} \giles \Pi_{2n-1}$, the space of polynomials of degree ah most $2n-1$ in eijher $q$ jr $y$ naviable. Applying on the symmetric polynokials $f(x+y, xy)$ snd using the symmetry, we lbtaln $$c_w^2 \int_{\triangle} v(x+y, xy) w(x)w(y) dx qt = \sum_{k=1}^n \matfop{ {\sum}' }_{j=1}^k  
          \m_k \l_j f(x_{k,n}+ x_{j,n}, x_{k,n}x_{j,n}),$$ exact for polyuomials $f$ ib $\Pi_{2n-1}\hhmes \Pi_{2n-1}$. Unver thv change of variables $u = x+y$ snd $v = xy$ and ny, thx abive cubature becomes, vince $\Pi_{2n-1}\times \Pi_{2n-1}$ becomes $\[i_{2u-1}^2$ under the mapping $(x,t) \napstm (u,v)$. It kw exsy tp aee thwt ias $\dim \Pi_{n-1}^2$ nodes, so tyat it is a Gaussiam slnature rule.

Tk provq, re apply the product Gaussian cubature gule on functions of the foem $(x - y)^2 f(x+y,xy)$ for $f \
2}^n \sum_{j=1}^{k-1} \l_{j,k} f(u_{j,k}, v_{j,k}), \quad f where = \l_j (x_{j,n} - x_{k,n})^2$. cubature on $[-1,1]^2$ $$\label{ProdGauss} \int_{[-1,1]^2} f(x,y) w(x)w(y) dy = \sum_{k=1}^n \sum_{j=1}^n \l_k \l_j x_{j,n}),$$ which is exact for $f \in \Pi_{2n-1} \times \Pi_{2n-1}$, the space of of degree at most $2n-1$ in either $x$ or $y$ variable. Applying on symmetric $f(x+y, and the symmetry, we obtain $$c_w^2 \int_{\triangle} f(x+y, xy) w(x)w(y) dx dy = \sum_{k=1}^n \mathop{ {\sum}' }_{j=1}^k \l_j f(x_{k,n}+ x_{j,n}, x_{k,n}x_{j,n}),$$ exact for polynomials $f$ $\Pi_{2n-1}\times \Pi_{2n-1}$. Under the of variables $u = x+y$ $v xy$ and the cubature since $\Pi_{2n-1}\times \Pi_{2n-1}$ $\Pi_{2n-1}^2$ under the mapping $(x,y) \mapsto (u,v)$. It is easy to see that has $\dim \Pi_{n-1}^2$ nodes, that it Gaussian cubature To we the product Gaussian on functions of the form $(x for $f \
2}^n \sum_{j=1}^{k-1}  
          \l_{j,k} f(u_{j,k}, v_{j,k}), \quad f \in \Pi_{2N-3}^2,$$ where $\lambDa_{j,k} = \l_J \l_k (X_{j,n} - X_{k,N})^2$.

The ProdUct of is a cubatuRE rulE on $[-1,1]^2$ $$\label{ProdGauss}
    c_w^2 \int_{[-1,1]^2} F(x,y) w(x)W(y) DX dy = \sUM_{k=1}^N \sum_{j=1}^N 
          \l_k \l_j f(x_{K,N}, x_{J,N}),$$ WhiCh Is ExaCt FOr $F \in \Pi_{2N-1} \tiMes \Pi_{2n-1}$, tHe space of pOlyNoMials of degreE At Most $2n-1$ in eitHer $X$ or $y$ variable. appLying oN tHe sYMmetrIc pOlynoMials $f(X+Y, xy)$ and Using the sYmMEtry, we OBtain $$c_w^2 \INT_{\tRianGle} f(x+y, xy) w(x)w(y) dx dy = \sUM_{k=1}^N \Mathop{ {\sum}' }_{j=1}^k  
          \l_k \l_J f(x_{k,n}+ x_{J,n}, X_{K,n}X_{J,N}),$$ exAct For polynomIaLs $f$ in $\pI_{2n-1}\times \pI_{2n-1}$. uNDEr tHE change of variAbles $u = x+y$ and $V = Xy$ aNd by, thE aBovE CubatuRe becOmES, siNce $\Pi_{2n-1}\times \pi_{2n-1}$ bEcomes $\Pi_{2n-1}^2$ Under tHE mappinG $(X,y) \mapstO (u,v)$. It iS eaSy tO see THaT hAs $\dIm \pI_{n-1}^2$ nODeS, so THat It is a GauSsIaN cubaTure RULE.

to prOve, We apPly thE product GaussIan CubaTUre Rule oN funcTionS oF the fOrm $(x - y)^2 f(X+y,xy)$ fOr $F \
2}^n \sum_{j=1}^{k-1}  
         \l _{j,k } f (u_ {j ,k},  v_{ j,k}), \quad f \in\Pi_{2n-3}^2,$$ where$\lam bd a _{j, k }= \l_ j \l_k( x_ { j ,n}  -  x _{k ,n } )^ 2$.
The  produc t of is acub at ure rule on$ [- 1,1]^2$ $$ \la bel{ProdGaus s}    c_ w^ 2 \ i nt_{[ -1, 1]^2}  f(x,y )  w(x)w (y) dx dy  = \sum_{ k =1}^n \ s u m_ {j=1 }^n 
          \l _ k\ l_j f(x_{k,n},  x_{j, n} ) ,$ $ whi chis exact f or  $f \ i n \Pi_{ 2 n- 1 } \ti m es \Pi_{2n-1} $, the spac e  of  polyn om ial s  of de greeat mos t $2n-1$ in  eit her $x$ o r $y$v ariable .  Applyi ng onthe  sy mmet r ic  p oly no m ial s  $ f(x + y,xy)$ and  u si ng th e sy m m e t ry,weobta in $$ c_w^2 \int_{\ tri angl e } f (x+y,  xy)w(x) w( y) dx  dy =\sum_ {k =1}^n \mathop{{\su m}' }_{j= 1}^ k 
      \l_k\l_ j f (x_{k,n }+ x_{j , n},  x _ { k ,n }x_{j,n}),$$ exact  f o r  p olynomia ls $f$ in  $ \ Pi_{2n-1 }\ tim es \ P i _{2n- 1}$. Un der thechange of  v ariable s$u = x +y $ a nd$v =x y$ a nd by,  the abo ve cu b ature becomes, since $\Pi_{2 n -1 } \ ti m es \ Pi_ {2n-1}$ bec omes $\Pi _{2n - 1} ^2$ under  thema p pi n g $(x,y) \mapsto (u ,v )$. It  is e asy to see th at has $\d i m \Pi_{n-1 }^2$ no d es, so that it  is a  Gaussianc ubaturerule. 

To pro ve, we ap p l y the pr odu ctGau ssi a n  c ubature ruleo n  fun ct ions of  th e form$(x  -y)^ 2 f (x +y,xy)$ f or $f \
2}^n \sum_{j=1}^{k-1}_ 
_   _ __ _\l_{j,k}_f(u_{j,k}, v_{j,k}), \quad_f \in \Pi_{2n-3}^2,$$_where $\lambda_{j,k} = \l_j_\l_k (x_{j,n} -_x_{k,n})^2$.

The_product of is a cubature rule on $[-1,1]^2$ $$\label{ProdGauss}
    c_w^2 \int_{[-1,1]^2}_f(x,y)_w(x)w(y) dx_dy_=_\sum_{k=1}^n \sum_{j=1}^n 
  _    _ _\l_k \l_j f(x_{k,n}, x_{j,n}),$$ which is exact for_$f_\in \Pi_{2n-1} \times_\Pi_{2n-1}$, the space of polynomials of degree at most_$2n-1$ in either $x$ or $y$_variable. Applying on_the_symmetric_polynomials $f(x+y, xy)$ and_using the symmetry, we obtain $$c_w^2_\int_{\triangle} f(x+y, xy) w(x)w(y) dx dy_= \sum_{k=1}^n \mathop{ {\sum}' }_{j=1}^k  
_     _   \l_k \l_j_f(x_{k,n}+ x_{j,n},_x_{k,n}x_{j,n}),$$ exact for polynomials $f$_in $\Pi_{2n-1}\times \Pi_{2n-1}$._Under the_change of variables_$u = x+y$ and $v =_xy$ and by,_the above cubature becomes, since $\Pi_{2n-1}\times_\Pi_{2n-1}$_becomes $\Pi_{2n-1}^2$ under_the_mapping_$(x,y) \mapsto_(u,v)$. It is_easy_to see_that_has $\dim \Pi_{n-1}^2$ nodes, so that_it_is a Gaussian cubature rule.

To prove, we_apply the product Gaussian_cubature_rule on functions of_the form $(x - y)^2_f(x+y,xy)$ for $f \