First-order methods with extended stability regions for solving electric circuit problems

Stability control of Runge-Kutta numerical schemes is studied to increase efficiency of integrating stiff problems. The implementation of the algorithm to determine coefficients of stability polynomials with the use of the GMP library is presented. Shape and size of the stability region of a method can be preassigned using proposed algorithm. Sets of first-order methods with extended stability domains are built. The results of electrical circuits simulation show the increase of the efficiency of the constructed first-order methods in comparison with methods of higher order.