It can be shown that the dynamics of magnetic perturbation at the nonmonoto-nous profile of a current in tokamak is qualitatively different from that of normal profiles. Behaviour of the perturbation depends on the helical magnetic currents localized mainly between the two singular surfaces. In this paper the numerical method developed to study the non-linear evolution of tearing modes is described in detail.