I need to solve the following equation for (x,y) $$x^y = y^x = 3$$
Everytime I run a numerical method for this problem, I get $$ (x,y) = (1.82546...,1.82546..) $$
I expect there to be a solution where $x \neq y$ but I cannot seem to find it by any means.
Does such a solution exist?
Assume $0<x<y$ and $x^y=y^x=3$. Then in particular $\frac{\ln x}x=\frac{\ln y}y$. The derivative of $f(t)=\frac{\ln t}t$ is $f'(t)=\frac{\frac1t\cdot t-1\cdot \ln t}{t^2}=\frac{1-\ln t}{t^2}>0$ for $0<t<e$, hence $f$ is injective on $(0,e]$. We conclude that $y>e$. Also note that for $x,y>0$ we have $x^y\le 1$ if $x\le 1$. Therefore $x>1$.
From $x>1$ and $y^x=3$ it follows that $y<3$. Then $$x=\sqrt[3]{x^3}\ge \sqrt[3]{x^y}=\sqrt[3]{y^x}>\sqrt[3]y>\sqrt[3]e>1+\frac13+\frac1{18}=\frac{25}{18}>\frac54.$$ Thus $$y^x>e^{\frac54}>1+\frac54+\frac{25}{32}=\frac{97}{32}>3,$$ contradiction.
Remark: The above does not require any calculator. All "numerical" calculations have been reduced to adding and multiplying fractions with manageably small denominator, and we use that for $t>0$ we have $e^t=\sum_{k=0}^\infty\frac{t^k}{k!}>1+t+\frac{t^2}2$.
Of course, if one trusts in calculator results, $\sqrt[3]e\approx 1.395612425$ and $e^{\sqrt[3]e}\approx4.037446449$ and that is of course also $>3$.