function [hh,Ht,H] = gethhn(P,m,n)
hh = zeros(n);
for i=1:n
  for j=i+1:n
    %h = etaij(P,m,n,i,j);
    h = mxetaij(P,m,n,i,j);  % the new one, as in Samson
    hh(i,j) = h;
  end
end
Ht = 1+sum(hh,2)';
H = max(Ht);
return