A parameter matrix decomposition (PMD) multiparameter modification formalism is presented for the resolution of nonlinear matrix equations. An efficient computational algorithm is developed for a general voltage-dependent

-terminal device model encompassing current nonlinear IC models, including diode, BJT, and MOSFET models. A PMD simulator is implemented for dc analysis of BJT analog and digital circuits. Previous problems of numerical instability and excessive computational cost using rank 1 modification are avoided.