 function f = log10x( x );
%
  [s1,s2]=size(x);
    f  =  log10((ones(s1,s2)+x)./(ones(s1,s2)-x));
return;
