function shat = damp(a,algor)
%
f = ff2(-a,1,algor);
shat = -log(abs(f));
return
