% UOFX The inverse function u(x) of x=u*log(u)
%
function y=uofx(x)
if x<=10
  u1=1.0125+.8577*x-.129013*x^2+.0208645*x^3 ...
     -.00176148*x^4+.000057941*x^5;
else
  x1=log(x);
  u1=(x/x1)/(1-log(x1)/(1+x1));
end
for n=1:4
  u0=u1;
  u1=(u0+x)/(1+log(u0));
end
y=u1;