code: def getSbus(self, buses=None):
          
          bs = self.buses if buses is None else buses
          s = array([self.s_surplus(v) / self.base_mva for v in bs])
          return s
docs: Returns the net complex bus power injection vector in p.u.