function y=logth(eps,x)
% logth - Natural log with lower threshold
%
% logth :: real~'lower threshold', real -> real.

x(x<eps)=eps; y=log(x);
