function maxr = submulthh(T,hh)
[n,n] = size(hh);
maxr = 0;
for i=1:n
  for j=i+1:n
    jrents = parentseq(T,j,1);
    for k=i+1:j-1
      if ismember(i,jrents) & ismember(k,jrents)
        r = hh(i,j) / (hh(i,k) * hh(k,j) );
        if r > 1.01
          keyboard
        end
        maxr = max(maxr,r);
      end
    end
  end
end
return


function pseq = parentseq(T,u,r)
v = u;
el = 0;
pseq = [u];
while v~=r
  el = el+1;
  rents = find(T(:,v));
  v = rents;
  pseq = [v pseq];
end
return