From Help, the procedure
StableMarriage was an element of the
Combinatorica, but it is available in the built-in
FindIndependentEdgeSet since Mathematica 10.
But it is not clear from the documentation how to use
FindIndependentEdgeSet as a Gale-Shapley algorithm straightforwardly. So, the procedure
FindIndependentEdgeSet can find an independent edge set and, for weighted graphs, maximal total-valued independent edge set. But this is no the same as a stable matching. How to solve it?
Explicitly: The input of the
StableMarriage is a pair of permutations describing male and female preferences. But that is not acceptable input for