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