def maximize(self, reaction):
self._prob.set_objective(self.flux_expr(reaction))<EOL>self._solve()<EOL>
Solve the model by maximizing the given reaction.

        If reaction is a dictionary object, each entry is interpreted as a
        weight on the objective for that reaction (non-existent reaction will
        have zero weight).