[tex]ab+c=x^2[/tex], [tex]bc+a=y^2[/tex], [tex]ca+b=z^2[/tex]
[tex]c=x^2-ab[/tex], [tex]b(x^2-ab)+a=y^2[/tex], [tex]a(x^2-ab)+b=z^2[/tex]
[tex]ab(x^2-ab)+a^2=y^2a[/tex], [tex]ab(x^2-ab)+b^2=z^2b[/tex]
[tex]z^2b-b^2=y^2a-a^2[/tex]
[tex]b(z^2-b)=a(y^2-a)[/tex]
[tex]b.\frac{z^2-b}{a}=a.\frac{y^2-a}{b}[/tex]
The fractions are integeres from the previous equations. If [tex]a^2|(z^2-b)[/tex] then [tex]a|x^2[/tex] and so [tex]a|c[/tex]. Otherwise [tex]a|b[/tex].
So now lets assume [tex]a\ge b \ge c[/tex]. So [tex]a[/tex] is either equal to [tex]b[/tex] or equal to [tex]c[/tex].
1. [tex]a=b[/tex]: then [tex]z=\pm y[/tex] then [tex]a^2+c=x^2[/tex]. The system of equations is
[tex]a^2+c=x^2[/tex], [tex]ac+a=y^2[/tex]. So [tex]a(x^2-a^2+1)=y^2[/tex]. If [tex]gcd(a,x^2-a+1)\ne 1[/tex] then [tex]a|(x^2+1)[/tex] and so [tex]a|(c+1)[/tex] so [tex]a=c+1[/tex]. So [tex]a^2+a-1=x^2[/tex]. Thne [tex]4a^2+4a-4=4x^2[/tex], so [tex](2a+1)^2-5=4x^2[/tex] and so [tex](2a+1+2x)(2a+1-2x)=5[/tex]. We get a system of equations. Clearly, the fisrt multiple is bigger than the second one (assuming x is positive, which is ok) so [tex]2a+1+2x=5[/tex], [tex]2a+1-2x=1[/tex] is the first option, we get [tex]2a+1=3[/tex], so [tex]a=1[/tex] Then [tex]b=1, c=0[/tex]. So [tex](1,1,0)[/tex] is a solution. The second option is [tex]2a+1+2x=-1[/tex], [tex]2a+1-2x=-5[/tex], so [tex]2a+1=-3[/tex], [tex]a=-2[/tex], [tex]b=-2[/tex], [tex]c=-3[/tex]. So [tex](-2,-2,-3)[/tex] is a solution.
2. [tex]a=c[/tex] then [tex]a=b=c[/tex] then [tex]3a^2=x^2[/tex] which has no nonzero solutions, so [tex](0,0,0)[/tex] is one solution.
NB: my solution is not complete since i have not tooak into consideration the negative options. When i wrote [tex]a|(c+1)[/tex] and assumes [tex]a \ge c[/tex] i concluded [tex]a=c+1[/tex]. But the other options are trivial. And the 2nd solution is not a solution after a check

NB2: i have forgotten the other case: [tex]gcd(a,x^2-a+1)=1[/tex].