function errmax=orthhestest(A)
siz=size(A);
n=siz(1);
la1=zeros(n,1);
la2=zeros(n,1);
la1=eig(A);
[H,U,beta]=orthhes(A);
la2=eig(H);
error=0;
for i=1:n
    v=la1(i)*ones(n,1);
    error=max(error,min(abs(la2-v)));
end
disp('orthhestest: maximale relative abweichung in den eigenwerten von A und H');
error=error/norm(A);
disp(num2str(error,14));
errmax=error;