code: def find_unconstrained_reactions(model):
      
      lower_bound, upper_bound = helpers.find_bounds(model)
      return [rxn for rxn in model.reactions if
              rxn.lower_bound <= lower_bound and
              rxn.upper_bound >= upper_bound]
docs: Return list of reactions that are not constrained at all.