The new effective algorithm of the numerical solution of Dirac equation is offered for modeling the self-coordinated electronic density of ions The Gunnarson-Lundkvist formula was used for exchange -correlation potential. In NaF crystal energy of compressed ions Na⁺ and F^- is determined and their radiuses are evaluated.