A modification of the trapezoidal method is suggested. The modification concludes in replace of the constant weighting coefficient by a linear function proportional to ratio of a carrying solution step to its given maximum value. The stability domain of the suggested method tends to the trapezoidal method stability domain with step decrease. The method is applied for analysis of narrow-band unstable electronic circuits.