My approach so far:
Let $x^2+y^2+xy=n^2,$ where $n\in\mathbb Z$.
$\implies (x+y) ^2-xy=n^2$
$\implies (x+y) ^2-n^2=xy$
$\implies (x+y+n) (x+y-n) =xy$
Only one case is possible:
When $x+y+n=xy$ and $x+y-n=1$. On adding, we get: $2(x+y) =xy+1$. Since LHS is an even number so is RHS. That's why, $xy$ must be odd$\implies x$ and $y$ both are odds. So I checked for few odd numbers and found out that for $x=3$ and $y=5,\;x^2+y^2+xy$ becomes $49$ which is a perfect square.. I didn't dare to check for more odds as there could be many...
Recently I found out that for $x=5$ and $y=16$(an even number), $x^2+y^2+xy=361=19^2.$ (Surprising!!)
So now I can say, I am stuck very badly.. All of my observations miserably failed...
Please suggest something!