Abstract: A stable recursive algorithm is presented for the transient simulation of interconnect systems in the high-speed VLSI and multichip modules (MCMs). The Taylor approximant are employed in the frequency domain, then a set of recursive formulas is derived by the inverse Laplace transform. The recursive convolution includes no unknowns, it is only related with the computed quantities. The algorithm concerns only the quantities at two ends of the lines, so it consumes less computer time and less memory in comparison with other methods. The algorithm is useful for both delay and crosstalk estimation and can be used to evaluate transient responses of high-speed interconnect circuits. Applications show that the method gives high accuracy.