The system of equations is symmetric: if you swap x and y, for example, then the third equation stays the same and the first two equations are swapped. Hence, we can assume without loss of generality that x≥y≥z. Then the system of equations becomes:
x2−yzy2−zxz2−xy=y−z+1,=x−z+1,=x−y+1.
Subtracting the second equation from the first, we obtain x2−y2+z(x−y)=y−x, or (x−y)(x+y+z+1)=0. This yields x=y or x+y+z=−1. Subtracting the third equation from the second, we obtain y2−z2+x(y−z)=y−z, or (y−z)(y+z+x−1)=0. This yields y=z or x+y+z=1.
We now distinguish two cases: x=y and x=y. In the first case, we have y=z, as otherwise we would have x=y=z for which the first equation becomes 0=1, a contradiction. Now it follows that x+y+z=1, or 2x+z=1. Substituting y=x and z=1−2x in the first equation yields x2−x(1−2x)=x−(1−2x)+1, which can be simplified to 3x2−x=3x, or 3x2=4x. We get x=0 or x=34. With x=0, we find y=0, z=1, but does not satisfy our assumption x≥y≥z. Thus, the only remaining possibility is x=34, which gives the triple (34,34,−35). We verify that this is indeed a solution.
Now consider the case x=y. Then we have x+y+z=−1, hence we cannot have x+y+z=1, and we see that y=z. Now x+y+z=−1 yields x+2z=−1, hence x=−1−2z. Now the first equality becomes (−1−2z)2−z2=1, which can be simplified to 3z2+4z=0. From this, we conclude that z=0 or z=−34. With z=0, we find y=0, x=−1, which does not satisfy our assumption x≥y≥z. Hence, the only remaining possibility is z=−34, and this gives rise to the triple (35,−34,−34). We verify that this is indeed a solution.
By also considering the permutations of these two solutions, we find all six solutions: (34,34,−35), (34,−35,34), (−35,34,34), (35,−34,−34), (−34,35,−34), and (−34,−34,35). □