code: def find_zero_constrained_reactions(model):
      
      return [rxn for rxn in model.reactions if
              rxn.lower_bound == 0 and
              rxn.upper_bound == 0]
docs: Return list of reactions that are constrained to zero flux.