  It is shown that the $gl(3)$ polynomial integrable system introduced by
Sokolov-Turbiner is equivalent to the $gl(3)$ quantum Euler-Arnold top in a
constant magnetic field. Their Hamiltonian and their 3rd order Integral can be
rewritten in terms of $gl(3)$ algebra generators. All $gl(3)$ generators can be
represented by the non-linear elements of the universal enveloping algebra of
the 5-dimensional Heisenberg algebra $h_5(\hat{p}_{1,2},\hat{q}_{1,2}, I)$.
Four different representations of the $h_5$ Heisenberg algebra are used: (I) by
differential operators in two real (complex) variables, (II) by
finite-difference operators on uniform or exponential lattices. We discovered
the existence of two 2-parametric bilinear and trilinear elements (denoted $H$
and $I$, respectively) of the universal enveloping algebra $U(gl(3))$ such that
their Lie bracket (commutator) can be written as a linear superposition of {\it
nine} so-called {\it artifacts} - the special bilinear elements of $U(gl(3))$,
which vanish once the representation of the $gl(3)$-algebra generators is
written in terms of the $h_5(\hat{p}_{1,2},\hat{q}_{1,2}, I)$-algebra
generators. In this representation all nine artifacts vanish, two of the
above-mentioned elements of $U(gl(3))$ commute(!). If $(\hat{p},\hat{q})$ are
represented by finite-difference/discrete operators on uniform or exponential
lattice, the Hamiltonian and the Integral of the 3-body elliptic Calogero model
become the isospectral, finite-difference operators on uniform-uniform or
exponential-exponential lattices (or mixed) with polynomial coefficients. If
$(\hat{p},\hat{q})$ are written in complex $(z, \bar{z})$ variables the
Hamiltonian corresponds to a complexification of the 3-body elliptic Calogero
model on $C^2$.
