In this paper, we present a method for obtaining stable time marching equations of distributed multiport networks with nonlinear loads. This method allows to highly decrease the number of frequency samples required by other stable TM methods. Stability is analyzed and ensured in the Z‐domain.