[Proof]Let
⎩⎨⎧x=n(4n2−1)+1,y=1−n(4n2−1),z=1−4n2n=1,2,⋯
Then, we have
x2(x−1)===[n(4n2−1)+1]2⋅n(4n2−1)[n(4n2−1)−1]2⋅n(4n2−1)+(4n2−1)2.4n2y2(1−y)+z2(1−z)
That is, x3+y3+z3=x2+y2+z2.
Therefore, the original equation has infinitely many integer solutions.