wolffd@0: function [indices, pot] = find_most_prob_entry(pot)
wolffd@0: % function [indices, pot] = find_most_prob_entry(pot)
wolffd@0: % function [indices, pot] = find_most_prob_entry(pot)
wolffd@0: % Find the indices of the argmax, and set all other enties to 0.
wolffd@0: 
wolffd@0: %indices = argmax(pot.T);
wolffd@0: [m i] = max(pot.T(:));
wolffd@0: indices = ind2subv(pot.sizes, i);
wolffd@0: pot.T = 0*myones(pot.sizes);
wolffd@0: pot.T(i) = m;