  The numerical solutions of the non-linear evolution equation are shown to
display the ``geometric'' scaling recently discovered in the experimental data.
The phenomena hold both for proton and nucleus targets for all $x$ below
$10^{-2}$ and $0.25 {\rm GeV^{2}}\le Q^2 \le 2.5\times10^3 {\rm GeV^{2}}$. The
scaling is practically exact (few percent error) in the saturation region. In
addition, an approximate scaling is found in the validity domain of the linear
evolution where it holds with about 10% accuracy.
  Basing on the scaling phenomena we determine the saturation scale $Q_s(x)$
and study both its $x$-dependence and the atomic number dependence for the
nuclei.

High Energy Physics - Phenomenology