docs: def PartialDynamicSystem(self, ieq, variable):
          """
          returns dynamical system blocks associated to output variable
          """
          if ieq == 0:
              # U1=0
              if variable == self.physical_nodes[0].variable:
                  v = Step('Ground', 0)
                  return[Gain(v, variable, 1)]