A special class of large-scale nonlinear systems of equations of block bordered structure is generated in VLSI circuit simulation from the sparse and irregular nature of the circuit matrices. Several Newton-like parallel methods for solving the block bordered equations are presented. Static and dynamic load balancing problems that arise in solving the circuit equations on multiprocessors are discussed. Computational results obtained using the BBN GP1000, a shared memory multiprocessor through a distributed architecture, are presented to show the effectiveness of dynamic and static load balancing. These experiments include a simulation of the op-amp 741 amplifier and a large analog filter simulation that leads to a multilevel block bordered system of nonlinear equations.< <ETX xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">></ETX>