function mesamson
global m n
rDG = 0; rGD = 0;
m=2; n=5;

while 1
  P = pnormdim(rand(m^n,1),1);
  [hh,Ht,H] = gethhn(P,m,n);
  hh = hh+eye(n);
  Dnorm = max(sum(hh,2));
  
  sh = sqrt(hh);
  Gnorm = sqrt(max(eig(sh'*sh)));

  rDG = max(rDG,Dnorm/Gnorm);
  rGD = max(rGD,Gnorm/Dnorm);
  
  fprintf('rDG = %1.5f         rGD = %1.5f   \n',rDG,rGD);
  
end


return