 has $\left(\begin{smallmatrix}n+1\\ j\\\end{smallmatrix}\right)+1$ terms, where $\left(\begin{smallmatrix}n+1\\j\\\end{smallmatrix}\right)$ is the size of the orbit of $\w_j$.

An efficient way to find the decompositions is to work with products of Weyl group orbits, rather than with orbit functions. Their decomposition has been studied, and many examples have been described in [@HLP]. It is useful to be aware of the congruence class of each product, because all of the orbits in its decomposition necessarily belong to that class. The *congruence number* $\#$ of an orbit $\lambda$ of $A_n$, which is also the congruence number of the orbit functions $C_\lambda$ and $S_\lambda$, specifies the class. It is calculated as follows, $$\begin{gathered}
\#(C_{(a_1,a_2,\dots,a_n)}(x))=\#(S_{(a_1,a_2,\dots,a_n)}(x))=\sum_{k=1}^n ka_k\mod (n+1).\end{gathered}$$ In particular, each $X_j$, where $j=1,2,\dots,n$, is in its own congruence class. During the multiplication, congruence numbers add up $\mod n+1$.

Polynomials in two and three variables originating from orbit functions of the simple Lie algebras $A_2$, $C_2$, $G_2$, $A_3$, $B_3$, and $C_3$ are obtained in the forthcoming paper [@NesterenkoPatera2009].

Exponential substitution
------------------------

There is another approach to multivariate orthogonal polynomials, which is also based on orbit functions. Such polynomials can be constructed by the continuous and invertible change of variables $$\begin{gathered}
\label{subst_exp}
y_j=e^{2\pi i x_j}, \quad x_j\in\R,\quad j=1,2,\dots,n.\end{gathered}$$ Consider an $A_n$ orbit function $C_\lambda(x)$, $S_\lambda(x)$ or $E_\lambda(x)$, when $\lambda$ is given in the $\w$-basis and $x$ is given in the $\alpha$-basis. Each
has $ \left(\begin{smallmatrix}n+1\\ j\\\end{smallmatrix}\right)+1 $ terms, where $ \left(\begin{smallmatrix}n+1\\j\\\end{smallmatrix}\right)$ is the size of the orbit of $ \w_j$. 

 An efficient way to witness the decomposition is to work with products of Weyl group orbits, quite than with orbit functions. Their decay has been studied, and many examples have been described in [ @HLP ]. It is utilitarian to be mindful of the congruence class of each intersection, because all of the orbits in its decomposition necessarily belong to to that class. The * congruence numeral * $ \#$ of an orbit $ \lambda$ of $ A_n$, which is also the congruity number of the orbit functions $ C_\lambda$ and $ S_\lambda$, specify the class. It is calculated as follow, $ $ \begin{gathered } 
 \#(C_{(a_1,a_2,\dots, a_n)}(x))=\#(S_{(a_1,a_2,\dots, a_n)}(x))=\sum_{k=1}^n ka_k\mod (n+1).\end{gathered}$$ In finical, each $ X_j$, where $ j=1,2,\dots, n$, is in its own congruence course. During the multiplication, congruence numbers add up $ \mod n+1$. 

 polynomial in two and three variables originating from orbit functions of the simple Lie algebras $ A_2 $, $ C_2 $, $ G_2 $, $ A_3 $, $ B_3 $, and $ C_3 $ are obtained in the forthcoming paper   [ @NesterenkoPatera2009 ]. 

 Exponential substitution 
 ------------------------ 

 There is another approach to multivariate orthogonal polynomials, which is also based on orbit functions. Such polynomial can be construct by the continuous and invertible change of variable $ $ \begin{gathered } 
 \label{subst_exp } 
 y_j = e^{2\pi i x_j }, \quad x_j\in\R,\quad j=1,2,\dots, n.\end{gathered}$$ think an $ A_n$ orbit function $ C_\lambda(x)$, $ S_\lambda(x)$ or $ E_\lambda(x)$, when $ \lambda$ is given in the $ \w$-basis and $ x$ is give in the $ \alpha$-basis. Each
 had $\left(\begin{smallmatrix}n+1\\ m\\\end{smallmatrix}\rntht)+1$ tecms, whede $\left(\bdgin{smallmatrix}n+1\\j\\\end{smallmatcix}\rught)$ us the size of the orbkt of $\w_j$.

Aj efficiwnt xay to find the vscomposlcions ls to xork with produgts of Weyl group orbits, saghzr than with orbit functions. Their dqcomposotlon has been sjudiec, and many examples have been described in [@HLP]. It is useful to be aware of the congrufnce class of each profuct, becausg alj of the orbigs in its decompositioh necessarily belong to that clxss. Tke *congruenxe nulter* $\#$ of an irbit $\lambda$ of $A_k$, which is alsp the congruenge nukbee of the orbit functimns $C_\lambda$ and $S_\lwmbda$, spewiyies the class. It is xaoculajed av fouoowr, $$\btgii{gafhered}
\#(F_{(a_1,a_2,\vots,a_n)}(x))=\#(S_{(a_1,a_2,\dkts,a_n)}(x))=\sum_{k=1}^n ka_k\mod (n+1).\end{gatherec}$$ Yb particular, sach $X_t$, rhere $j=1,2,\dots,n$, is in its own congruence cnasa. During the multiplicarion, congruence numbets add up $\iod n+1$.

Polynomials in two and three variables origitatinj wron jfvih functions of the simple Lie algebras $A_2$, $C_2$, $G_2$, $W_3$, $B_3$, akd $C_3$ are obtaineb in the forthcokijg kaper [@NesterenkuPaterc2009].

Sxlonential substituhion
------------------------

Thete is qnother akproavh to multivariate orthogonql polynomiajw, which is also baded on orbic funcjions. Xuch polynomials can be conatructed by the contjvuous and invertkblv chdnge of variables $$\begin{gatrered}
\labeo{subxt_exp}
y_j=d^{2\pi o x_j}, \qtad x_j\in\R,\qkad j=1,2,\dots,n.\end{gathered}$$ Cojsidet an $A_t$ orbit fujction $C_\lambda(x)$, $S_\lambda(x)$ or $E_\lambda(x)$, when $\lamndd$ iv given nn the $\w$-basis and $x$ is given in tke $\alpha$-yasis. Dach
has $\left(\begin{smallmatrix}n+1\\ j\\\end{smallmatrix}\right)+1$ terms, where $\left(\begin{smallmatrix}n+1\\j\\\end{smallmatrix}\right)$ is of orbit of An efficient way to with products of group orbits, rather with orbit functions. Their decomposition has studied, and many examples have been described in [@HLP]. It is useful to aware of the congruence class of each product, because all of the orbits its necessarily to class. The *congruence number* $\#$ of an orbit $\lambda$ of $A_n$, which is also the congruence of the orbit functions $C_\lambda$ and $S_\lambda$, specifies class. It is calculated follows, $$\begin{gathered} \#(C_{(a_1,a_2,\dots,a_n)}(x))=\#(S_{(a_1,a_2,\dots,a_n)}(x))=\sum_{k=1}^n ka_k\mod (n+1).\end{gathered}$$ particular, $X_j$, where is its congruence class. During multiplication, congruence numbers add up $\mod n+1$. Polynomials in two and three variables originating from orbit functions the simple $A_2$, $C_2$, $A_3$, and are obtained in paper [@NesterenkoPatera2009]. Exponential substitution ------------------------ There to multivariate orthogonal polynomials, which is also based orbit functions. polynomials can be constructed by the and invertible change of variables $$\begin{gathered} \label{subst_exp} y_j=e^{2\pi x_j}, \quad x_j\in\R,\quad j=1,2,\dots,n.\end{gathered}$$ Consider an $A_n$ orbit function $C_\lambda(x)$, $S_\lambda(x)$ or $E_\lambda(x)$, when $\lambda$ in the $\w$-basis and is given in $\alpha$-basis.
 has $\left(\begin{smallmatrix}n+1\\ j\\\End{smallmaTrix}\rIghT)+1$ teRmS, wheRe $\leFt(\begin{smallmaTRix}n+1\\J\\\end{smallmatrix}\right)$ is The siZe OF the ORbIt of $\w_J$.

An effiCIeNT Way To FiNd tHe DEcOmposItiOns is to Work with prOduCtS of Weyl group ORbIts, rather tHan With orbit funCtiOns. TheIr DecOMposiTioN has bEen stuDIed, and Many exampLeS Have beEN descriBED iN [@HLP]. it is useful to be awaRE oF The congruence cLass of EaCH pRODucT, beCause all of ThE orbiTS in its dECoMPOSitIOn necessarily Belong to thaT ClaSs. The *cOnGruENce numBer* $\#$ of An ORbiT $\lambda$ of $A_n$, WhicH is also thE congrUEnce numBEr of the Orbit fUncTioNs $C_\lAMbDa$ And $s_\lAMbdA$, SpEciFIes The class. it Is CalcuLateD AS FOlloWs, $$\bEgin{GatheRed}
\#(C_{(a_1,a_2,\dots,a_n)}(x))=\#(s_{(a_1,a_2,\Dots,A_N)}(x))=\sUm_{k=1}^n kA_k\mod (N+1).\end{GaThereD}$$ In parTiculAr, Each $X_j$, where $j=1,2,\dotS,n$, is In its own cOngRuEncE cLass. DURing thE muLtiPlicatiOn, congrUEncE nUMBErS add up $\mod n+1$.

PolynomiAlS IN tWo and thrEe variABlEs ORiginatiNg FroM orbIT FunctIons OF tHe simple lie algEBrAs $a_2$, $C_2$, $G_2$, $A_3$, $B_3$, anD $C_3$ Are obtAiNed In tHe forTHcomIng papEr [@NesterEnkoPATera2009].

ExponentiaL Substitution
------------------------

THErE IS aNOtheR apProach to mulTivaRIate OrthOGoNal POlynoMials, WhICh IS also based on orbit fuNcTions. SUch poLynomials can bE constructED BY the contInuoUS aND invertible chaNge of Variables $$\bEGin{gatheRed}
\laBel{subst_Exp}
y_j=e^{2\pi i X_J}, \Quad x_j\in\r,\quAd j=1,2,\DotS,n.\eND{GaThered}$$ ConsideR AN $A_n$ oRbIt functIon $c_\lambda(X)$, $S_\lAmbDa(x)$ Or $E_\LaMbda(x)$, when $\Lambda$ is GiVeN iN tHe $\w$-Basis ANd $x$ is givEn In tHe $\AlpHa$-basIS. Each
 has $\left(\begin{smallma trix}n+1\\  j\\\ end {sm al lmat rix} \right)+1$ ter m s, w here $\left(\begin{sma llmat ri x }n+1 \ \j \\\en d{small m at r i x}\ ri gh t)$  i s  t he si zeof theorbit of $ \w_ j$ .

An effici e nt  way to fi ndthe decompos iti ons is  t o w o rk wi thprodu cts of Weyl g roup orbi ts ,  rathe r  than w i t horbi t functions. Thei r  d e composition ha s been  s t ud i e d,and  many exam pl es ha v e beend es c r i bed in [@HLP]. It  is usefult o b e awar eoft he con gruen ce cla ss of eachprod uct, beca use al l  of the orbitsin its  de com posi t io nnec es s ari l ybel o ngto thatcl as s. Th e *c o n g r uenc e n umbe r* $\ #$ of an orbi t $ \lam b da$  of $ A_n$,  whi ch  is a lso th e con gr uence number of  the  orbit fu nct io ns$C _\lam b da$ an d $ S_\ lambda$ , speci f ies  t h e cl ass. It is calcula te d as  follows , $$\b e gi n{ g athered} 
\ #(C _{(a _ 1 ,a_2, \dot s ,a _n)}(x)) =\#(S_ { (a _1 ,a_2,\d ot s,a_n) }( x)) =\s um_{k = 1}^n  ka_k\ mod (n+1 ).\en d {gathered}$$ I n  particular,e ac h $X _ j$,whe re $j=1,2,\ dots , n$,is i n  i tso wn co ngrue nc e  c l ass. During the mul ti plicat ion,congruence nu mbers addu p $\mod n+ 1$.Po l ynomials in tw o and  three var i ables or igina ting fro m orbit f u n ctions o f t hesim ple L ie  algebras $A_ 2 $ , $C _2 $, $G_2 $,$A_3$,$B_ 3$,  an d $ C_ 3$ are ob tained i nth efo rth comin g  paper [ @N est er enk oPate r a2009] .

Ex pone nt ia l  su bstitut i on - ---- -- -- ---- --- -- ----- -

T h ere  is ano ther appr oac h  tomu lt ivariat e orthogonalpo lynomials,  w hic h is a l s o basedon orbit functions. Suc h  polyno mia ls ca n be  construc ted  by th e c o ntinuo us and  inve rt ibl e chang e of  va ri ables $$\b e g in{ gathe re d}
\ label{s ubst_exp}
y_j=e^{2 \ pii x_j}, \quad  x_ j\in \ R ,\ qua d  j = 1,2 ,\ d ots , n .\end{gathered} $$ Conside ra n$A_n$ orbi t  fu nc tion $C _\lambd a(x)$ ,  $S_\la mbda(x)$or $E_\la mb da(x ) $ , w hen $\lamb da$ is g iven in t h e $\w $ -b asisand  $x$ i sgiv en in  the $ \ alp ha$-b asis.Ea ch
 has_$\left(\begin{smallmatrix}n+1\\ j\\\end{smallmatrix}\right)+1$_terms, where $\left(\begin{smallmatrix}n+1\\j\\\end{smallmatrix}\right)$ is_the size_of_the orbit_of_$\w_j$.

An efficient way_to find the_decompositions is to work_with products of_Weyl_group orbits, rather than with orbit functions. Their decomposition has been studied, and many_examples_have been_described_in_[@HLP]. It is useful to_be aware of the congruence_class of_each product, because all of the orbits in_its_decomposition necessarily belong_to that class. The *congruence number* $\#$ of an_orbit $\lambda$ of $A_n$, which is_also the congruence_number_of_the orbit functions $C_\lambda$_and $S_\lambda$, specifies the class. It_is calculated as follows, $$\begin{gathered}
\#(C_{(a_1,a_2,\dots,a_n)}(x))=\#(S_{(a_1,a_2,\dots,a_n)}(x))=\sum_{k=1}^n ka_k\mod_(n+1).\end{gathered}$$ In particular, each $X_j$, where $j=1,2,\dots,n$,_is in its own congruence class._During the multiplication, congruence numbers_add up_$\mod n+1$.

Polynomials in two and_three variables originating_from orbit_functions of the_simple Lie algebras $A_2$, $C_2$, $G_2$,_$A_3$, $B_3$, and_$C_3$ are obtained in the forthcoming_paper [@NesterenkoPatera2009].

Exponential_substitution
------------------------

There is another_approach_to_multivariate orthogonal_polynomials, which is_also_based on_orbit_functions. Such polynomials can be constructed_by_the continuous and invertible change of variables_$$\begin{gathered}
\label{subst_exp}
y_j=e^{2\pi i x_j}, \quad_x_j\in\R,\quad_j=1,2,\dots,n.\end{gathered}$$ Consider an $A_n$_orbit function $C_\lambda(x)$, $S_\lambda(x)$ or_$E_\lambda(x)$, when $\lambda$ is given in_the $\w$-basis_and $x$_is given in the $\alpha$-basis. Each