Two multi-steps methods: Adams and BDF was applied for calculation of phase portrait of incommensurate superstructure with the Lifshitz’s invariant. The Python language was used to create the correspond application. The result of numerical experiments was gives an excellent result in the calculation of complex dynamic systems.