Không mất tính tổng quát, giả sử \(x\ge y\ge z\ge1\).
Khi đó ta có: \(13=xyz+x^2+y^2+z^2\ge z^3+3z^2\)
suy ra \(z=1\).
\(12=xy+x^2+y^2\ge y^2+y^2+y^2=3y^2\)
\(\Rightarrow y=1\)hoặc \(y=2\).
Với \(y=1\): \(x^2+1+1+x=13\Leftrightarrow x^2+x-11=0\)không có nghiệm nguyên dương.
Với \(y=2\): \(x^2+2^2+1^2+1.2.x=13\Leftrightarrow x^2+2x-8=0\Leftrightarrow\left(x-2\right)\left(x+4\right)=0\)
\(\Rightarrow x=2\)thỏa mãn.
Vậy phương trình có nghiệm là \(\left(1,2,2\right)\)và các hoán vị.