function d=hilbertmetric(x,y)
% function d=hilbertmetric(x,y)

eigxinvy=eig(inv(x)*y);
d=log(max(eigxinvy)/min(eigxinvy));
end