code: def run_fba_minimized(self, reaction):
          
  
          epsilon = self._args.epsilon
          solver = self._get_solver()
  
          p = fluxanalysis.FluxBalanceProblem(self._mm, solver)
  
          start_time = time.time()
  
          # Maximize reaction flux
          try:
              p.maximize(reaction)
          except fluxanalysis.FluxBalanceError as e:
              self.report_flux_balance_error(e)
  
          fluxes = {r: p.get_flux(r) for r in self._mm.reactions}
  
          # Run flux minimization
          flux_var = p.get_flux_var(reaction)
          p.prob.add_linear_constraints(flux_var == p.get_flux(reaction))
          p.minimize_l1()
  
          logger.info('Solving took {:.2f} seconds'.format(
              time.time() - start_time))
  
          count = 0
          for reaction_id in self._mm.reactions:
              flux = p.get_flux(reaction_id)
              if abs(flux - fluxes[reaction_id]) > epsilon:
                  count += 1
              yield reaction_id, flux
          logger.info('Minimized reactions: {}'.format(count))
docs: Run normal FBA and flux minimization on model.