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