Homotopy solution to non-linear evolution for heavy nuclei

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