def moma(self, wt_fluxes):
reactions = set(self._adjustment_reactions())<EOL>v = self._v<EOL>obj_expr = <NUM_LIT:0><EOL>for f_reaction, f_value in iteritems(wt_fluxes):<EOL><INDENT>if f_reaction in reactions:<EOL><INDENT>obj_expr += (f_value - v[f_reaction])**<NUM_LIT:2><EOL><DEDENT><DEDENT>self._prob.set_objective(obj_expr)<EOL>self._solve(lp.ObjectiveSense.Minimize)<EOL>
Minimize the redistribution of fluxes using Euclidean distance.

        Minimizing the redistribution of fluxes using a quadratic objective
        function. The distance is minimized by minimizing the sum of
        (wild type - knockout)^2.

        Args:
            wt_fluxes: Dictionary of all the wild type fluxes that will be
                used to find a close MOMA solution. Fluxes can be expiremental
                or calculated using :meth: get_fba_flux(objective).