The numerical simulation of 2D simple square lattice of Ising superspins, that interact via the Ruderman–Kittel–Kasuya–Yosida (RKKY) exchange interaction is performed. It is found, that at low temperatures in the researched systems predominantly realized labyrinth domain structure. The modeled image of magnetization distribution does not coincide with the image of energy distribution. Hysteresis loops of the samples and labyrinth domain structures similarly obtained by Monte Carlo method, are also observed in physical experiments.