function [y] = log10(x)
%y=log10(x)
%compute the logarithm, base 10, of x

y = log(x+eps) ./ log(10);
