docs: def PartialDynamicSystem(self, ieq, variable):
          """
          returns dynamical system blocks associated to output variable
          """
          if ieq == 0:
              # U1=0
              if variable == self.variables[0]:
                  return[Gain(self.commands[0], variable, -self.Tmax)]