This paper presents a hierarchical method for performing symbolic simulation in the time domain using a state-variable approach. The paper illustrates the development of the symbolic state-variable matrices for the subcircuits of a partitioned network, and the methodology for the hierarchical analysis. The resultant state equations are combined to produce the compact state-variable equations for the entire system. This hierarchical method exhibits a large reduction in the number of operations and in the size of the characteristic matrices needed to characterize the statevariable equations for the network. This enables the application of the method to large-scale networks. The method accommodates a model for ideal operational amplifiers.