Homotopy solution to the nonlinear evolution for heavy nuclei

In the paper we suggest the homotopy method for solving of the nonlinear evolution equation. This method consists of two steps. First is the analytical solution for the linearized version of the nonlinear evolution deep in the saturation region. Second, the perturbative procedure is suggested to take into account the remaining parts of the nonlinear corrections. It turns out that these corrections are rather small and can be estimated in the regular iterative procedure.