Abstract. A class of numerical combination methods is developed where the discretization error decreases as the order of combination increases. A combination is obtained for which the discretization error is absent up to the second order of smallness. The analytical error estimates are tested in the analysis of conservative systems without losses, quartz oscillators, and high-Q systems with long transients.

Keywords: dynamic system, combinational numerical method, discrete model, discretization error, computer modeling.