function shat = damp3(a,algor)
%
f = ff3(-a,1,algor);
shat = -log(max(abs(eig(f))));
return
