The software is presented for simulation of anisoplanatic effect and its influence on performance of adaptive optical phase conjugation system in inhomogeneous turbulent atmosphere. Atmospheric turbulence was simulated with the help of a set of moving random phase screens with arbitrary statistics. Both reference and target are supposed to be the point light sources. To simulate atmospheric turbulence we applied the concept of a number of moving random phase screens with Kolmogorov spectrum. In our investigations we used the model of Shack-Hartmann wavefront sensor and the ideal model of wavefront adaptive mirror that is assumed to reproduce a given number of Zernike polynomials without time delays. The designed software allows to calculate instantaneous and average values of phase correction errors at a different angle between a reference beacon and target source. Simulations can be made with a broad range of parameters of adaptive system and atmospheric turbulence. The approach enables us to estimate residual aberrations as well as to calculate instant parameters of system performance - point spread function (PSF), optical transfer function (OTF) - and system isoplanatic angle. The model of system allows to change the control algorithm of phase corresction. Both common phase conjugation and weighted phase conjugation algorithm have been tested.