If we denote a=x−y and b=x+y, the equation rewrites as a2+b2=2a3, 2a−1= b2a2.
We see that 2a-1 is the square of a rational number, hence it is the square of an (odd) integer. Let 2a-1 = (2n+1)2.
It follows that a=2n2+2n+1 and b=a(2n+1)=(2n+1)(2n2+2n+1).
Finally, we obtain x=2n3+4n2+3n+1, y=2n3+2n2+n, where n is an arbitrary integer number.