function d = dish(gam)
% the "disharmony" of a column g
d = 0;

for col = 1:size(gam,2)
  g = gam(:,col);
  dij = 0;
  for i=1:length(g)
    for j=i+1:length(g)
      dij = dij + abs(g(i)-g(j));
      %dij = max(dij,abs(g(i)-g(j)));
    end
  end
  %d = max(d,dij);
  d = d + dij;
end

%d = d + max(abs(sum(gam')));
d = d + sum(abs(sum(gam')));

return