function tweaksum(n)
H = fillHn(n);
G = sqrt(H)'*sqrt(H);
G1 = fillGn(n);
totdif(G,G1)
return

function G = fillGn(n)
G = zeros(n);
for k=1:n
  for l=1:n
    s = 0;
    for p=1:min(k,l)-1
      s = s + 1/(n-p);
    end
    g = (k==l) + (k<l)/sqrt(n-k+eps) + (l<k)/sqrt(n-l+eps) + s;
    G(k,l) = g;
  end
end
return

function H=fillHn(n)
H = zeros(n);
for i=1:n
  row = ones(1,n-i);
  row = row/sum(row);
  H(i,i+1:n) = row;
end
H = H + diag(ones(n,1),0);
return