- 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.25GeV^2≤Q^2≤2.5\times10^ 3GeV^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 ourselves 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.