  In this paper the semi-classical approach to the solution of non-linear
evolution equation is developed. We found the solution in the entire kinematic
region to the non-linear evolution equation that governs the dynamics in the
high parton density QCD. The large impact parameter ($b_t$) behavior of the
solution is discussed as well as the way how to include the non-perturbative
QCD corrections in this region of $b_t$. The geometrical scaling behavior and
other properties of the solution in the saturation (Color Glass Condensate)
kinematic domain are analyzed. We obtain the asymptotic behavior for the
physical observables and found the unitarity bounds for them.
