A new method capable of simulating transient behavior of tightly coupled circuits with high precision and speed is proposed. It is based on the tranrformation of the differential equations into integral equations. The method is semi-analytic and is of the third order. It gives a net speed advantage (about one order of magnitude) over classical methods, especially for high precision (typically analog) circuits. I3 also gives an a priori error estimate which will reduce rejected steps.