  The solutions of the Balitsky-Kovchegov evolution equations are studied
numerically and compared with known analytical estimations. The rapidity and
nuclear size dependences of the saturation scale are obtained for the cases of
fixed and running coupling constant. These same dependences are studied in
experimental data, on lepton-nucleus, deuteron-nucleus and nucleus-nucleus
collisions, through geometric scaling and compared with the theoretical
calculations.
