  We present the numerical solution of the non-linear evolution equation for
DIS on nuclei for $x = 10^{-2} \div 10^{-7}$. We demonstrate that the solution
to the non-linear evolution equation is quite different from the Glauber -
Mueller formula which was used as the initial condition for the equation.
  We illustrate the energy profit for performing DIS experiments on nuclei.
  However, it turns out that the gain is quite modest: $x_{Au} \simeq 5 x_{\rm
proton} $ for the same parton density.
  We find that the saturation scale $Q^2_s \propto A^{{1/3}}$. For gold the
saturation scale $Q_{s,Au} \simeq 1.5 GeV$ at $x= 10^{-3}$. Such a large value
leads to considerable contribution of the high density QCD phase to RHIC data
and reveals itself in essential damping for both $xG_A$ and $F_{2A}$.
