In order to obtain the stable stationary mass distribution, we apply a new iterative method, inspired by a previous one used in 2D turbulence, which ensures entropy increase and, hence, convergence towards an equilibrium state. Applying new iterative method, we analyze the phase transition and the difference between microcanonical and canonical ensemble in an intermediate energy region.
