ABSTRACT. Infinite series of boundary conditions that are consistent with even-order higher symmetries find ensure the integrability of a Burgers type equation are constructed.KEY WORDS: integrable equations, higher symmetries, boundary conditions, recursion operators.The nonlinear equations of mathematical physics contain the important class of integrable equations that can be either integrated by the inverse problem method or Linearized by differential substitutions (S-and C-integrable equations, respectively, in Calogero's terminology [1]). From the viewpoint of the symmetry approach [2], such systems axe distinguished by the fact that they have infinitely many higher symmetries (commuting flows). From the analytical viewpoint, the integrable equations are characterized by the fact that they admit a wide class of exact and asymptotic solutions. The presence of higher symmetries permits one to reduce the construction of soliton and finite-gap solutions of partial differential equations to solving ordinary differential equations.Sometimes, the solution is subjected to additional boundary conditions. In that connection, the following question arises: What boundary conditions can be satisfied by an exact solution (constructed with the help of higher symmetries). In other words, we must describe boundary conditions consistent with the higher symmetries (the precise definition will be given later). The first examples of nontrivial boundary conditions of this sort were obtained by Moser [3] for the Toda lattice and by Sklyanin [4] for nonlinear Schzfidinger type equations.In [5], an algorithm for constructing integrable boundary condition was suggested. This algorithm involves finding differential reductions for some systems of differential equations that are constructed in a unified way from the higher symmetries of the original equation. In the present paper, we consider the integrable many-component Burgers type equations constructed in [6]. These are systems of N equations of the form u,~ +2ajku u, + bjk,=uJuku '~, i = 1,... ,N, (0.1) where u ~ = u~(t, z), {a~k, b}k,~ } is a collection of constants satisfying certain polynomial identities, and summation from 1 to N ox-er repeated indices is assumed. For these systems, we find an infinite series of boundary conditions consistent with the higher symmetries and having the formwhere U = (u I , ..., u N) and Uk = cgkU/cDz ~ . The boundary conditions from this infinite series can be obtained by applying a differential operator to the simplest of these conditions (by analogy with the recursion operator for symmetries, this operator is naturally referred to as the recursion operator for integrable boundary conditions). Not that so far the authors do not know any similar recursion operators for boundary conditions for other integrable equations. Integrable many-component polynomial systems are closely related to nonassociative algebraic structures; that is why the statement of results and the proofs essentially use the language of nonassociative algebras.