There is, if you have the Symbolic Math Toolbox:
syms x y
Eqn = y*x == 1 + log(y/x);
y = solve(Eqn, y);
figure(1)
fplot(real(y), [0 1], '-b')
hold on
fplot(imag(y), [0 1], '-r')
hold off
grid
set(gca, 'YScale','log')
xlabel('X')
ylabel('Y')
legend('Re\{y\}', 'Im\{y\}')
y =
-wrightOmega(- log(-1/x^2) - 1)/x
