  We consider the gl(N)-invariant Calogero-Sutherland Models with N=1,2,3,...
in a unified framework, which is the framework of Symmetric Polynomials. By the
framework we mean an isomorphism between the space of states of the
gl(N)-invariant Calogero-Sutherland Model and the space of Symmetric Laurent
Polynomials. In this framework it becomes apparent that all gl(N)-invariant
Calogero-Sutherland Models are manifestations of the same entity, which is the
commuting family of Macdonald Operators. Macdonald Operators depend on two
parameters $q$ and $t$. The Hamiltonian of gl(N)-invariant Calogero-Sutherland
Model belongs to a degeneration of this family in the limit when both $q$ and
$t$ approach the N-th elementary root of unity. This is a generalization of the
well-known situation in the case of Scalar Calogero-Sutherland Model (N=1). In
the limit the commuting family of Macdonald Operators is identified with the
maximal commutative sub-algebra in the Yangian action on the space of states of
the gl(N)-invariant Calogero-Sutherland Model. The limits of Macdonald
Polynomials which we call gl(N)-Jack Polynomials are eigenvectors of this
sub-algebra and form Yangian Gelfand-Zetlin bases in irreducible components of
the Yangian action. The gl(N)-Jack Polynomials describe the orthogonal
eigenbasis of gl(N)-invariant Calogero-Sutherland Model in exactly the same way
as Jack Polynomials describe the orthogonal eigenbasis of the Scalar Model
(N=1). For each known property of Macdonald Polynomials there is a
corresponding property of gl(N)-Jack Polynomials. As a simplest application of
these properties we compute two-point Dynamical Spin-Density and Density
Correlation Functions in the gl(2)-invariant Calogero-Sutherland Model at
integer values of the coupling constant.
