This is the system of equations: $$\sqrt { x } +y=7$$ $$\sqrt { y } +x=11$$
Its pretty visible that the solution is $(x,y)=(9,4)$
For this, I put $x={ p }^{ 2 }$ and $y={ q }^{ 2 }$. Then I subtracted one equation from the another such that I got $4$ on RHS and factorized LHS to get two factors in terms of $p$ and $q$.
Then $4$ can be represented as $2*2$, $4*1$ or $1*4$. Comparing the two factors on both sides, I got the solution.
As you can see, the major drawback here is that I assumed this system has only integral solutions and then went further. Is there any way I can prove that this system indeed has only integral solutions or is there any other elegant way to solve this question?




You want \begin{cases} x=(7-y)^2\\ y=(11-x)^2\\ 0\le x\le 11\\ 0\le y\le 7 \end{cases}
The equation becomes $$ x=49-14(121-22x+x^2)+(121-22x+x^2)^2 $$ which reduces to $$ (x-9)(x^3-35x^2+397x-1444)=0 $$ (courtesy of WolframAlpha). The polynomial $f(x)=x^3-35x^2+397x-1444$ has at least one real root. It has indeed three. One of them satisfies the condition $0\le x\le 11$ and it's approximately $7.87$. With this value of $x$ we get $y\approx 9.79$ that doesn't satisfy $0\le y\le 7$.
One can be more precise: call $\alpha$ the least root of $f$. Then $7<\alpha<8$, so that $3<11-\alpha<4$ and so $9<(11-\alpha)^2<16$, which shows that the limitation on $y$ is not fulfilled.
How can you know this? Compute
Thus you know that the three roots of $f$ are $7<\alpha<8$, $12<\beta<13$ and $13<\gamma<15$.