I am solving this Diophantine equation
$$xy(x^2-y^2)=pq(p^2-q^2)=rs(r^2-s^2)$$
where $(x^2+y^2)^2+(p^2+q^2)^2=2(r^2+s^2)^2$,
but I am not able to proceed.
I use Euler method $pq(p^2-q^2) = rs(r^2-s^2)$.
If $s = p$, the equation reduces to solving $p^2 = q^2-qr+r^2$, two solutions of which are,
$$\{p,q,r\}=\{u^2+3v^2, u^2+2uv-3v^2, 4uv\}$$
$$\{p,q,r\}=\{u^2+3v^2, u^2-2uv-3v^2,u^2+2uv-3v^2\}.$$
But L.J.Lander gives the solution
$$\{p,q,r,s\} = \{v^5-2v, v^5+v, -2v^4+1, v^4+1\}.$$
But I am to find integer solutions for the above three equivalent rather than two with specific conditions $$(x^2+y^2)^2+(p^2+q^2)^2=2(r^2+s^2)^2$$ which will help me to know more about arithmetic progression between square. So, please provide solution or guidance to solve this type of diophantine equation or write a computer program to get integer solutions.