code: def vars(self):
          
          return self.independent_vars + self.dependent_vars + [self.sigmas[var] for var in self.dependent_vars]
docs: :return: Returns a list of dependent, independent and sigma variables, in that order.