These are then in turn passed along the scipy library “optimize.root” function which solves the system, generating a list of branch currents and node voltages.