Solve the system of equations with $x=y$ Solve the system of equations: $\left\{\begin{array}{l}\sqrt{x^2+(y-2)(x-y)}+\sqrt{xy}=2y\\\sqrt{xy+x+5}-\dfrac{6x-5}{4}=\dfrac{1}{4}\left(\sqrt{2y+1}-2\right)^2\end{array}\right.$

I used wolframalpha.com and got the only solution: $(x;y)=(4;4)$
And I guess that we can get $x=y$ from first equation.
And this is my try

We have $\sqrt{x^2+(y-2)(x-y)}-y+\sqrt{xy}-y=0$
$\Leftrightarrow (x-y)\left(\dfrac{x+2y-2}{\sqrt{x^2+(y-2)(x-y)}+y}+\dfrac{y}{\sqrt{xy}+y}\right)=0$

But I can't prove that $\dfrac{x+2y-2}{\sqrt{x^2+(y-2)(x-y)}+y}+\dfrac{y}{\sqrt{xy}+y}\ne0$.
So who can help me?

If we allow the non-principal square root, you can actually have more than one solution to the system,
$$\left\{\begin{array}{l}\sqrt{x^2+(y-2)(x-y)}\color{red}\pm\sqrt{xy}=2y\\\sqrt{xy+x+5}-\dfrac{6x-5}{4}=\dfrac{1}{4}\left(\sqrt{2y+1}-2\right)^2\end{array}\right.$$
The positive case is solved by $x=y = 4$, but the negative case can be solved by,
$$x=2.771913334036360207047659037932620797468\dots\\ y=0.4518271741040497609916929241427121641268\dots$$
which are roots of $12$-deg equations and can be found by using resultants.
