% SUOFX Symbolic counterpart of uofx
%
function y=suofx(dig,x)
syms u0 u1 
digits(dig)
x0=double(x);
if x0<=10
  du1=1.0125+.8577*x0-.129013*x0^2+.0208645*x0^3 ...
     -.00176148*x0^4+.000057941*x0^5;
else
  x1=log(x0);
  du1=(x0/x1)/(1-log(x1)/(1+x1));
end
u1=vpa(du1);
x=vpa(x);
%for n=1:5
for n=1:6
  u0=u1;
  u1=vpa((u0+x)/(1+log(u0)));
end
y=u1;