A numerical value of the estimate of global discretization error is used in the simulation process of transient and steady-state modes in dynamic systems. This allowed the calculation accuracy to by increased significantly with low computational expanses.