Suppose that the equation has a solution (x1,y1) with x1≤y1
Clearing denominators, we can write the equation under the form x2−(4y−1)x+y2+y=0
that is, a quadratic in x. One of the roots is x1, therefore. by Vieta's theorem, the second one is 4y1−1−x1.
Observe that 4y1−1−x1≥4x1−1−x1=3x1−1≥2x1>0
Hence, 4y1−1−x1 is a positive integer. It follows that (4y1−1−x1,y1)
is another solution of the system, Because the equation is symmetric, we obtain that (x2+y2)=(y1,4y1−1−x1)
is also a solution.To end the proof, observe that x2+y2=5y1−1−x1+y1
and that (1,1) is a solution. Thus , we can generate infinitely many solutions: (1,1)→(1,2)→(2,6)→