`TH_{1}:x=y=0`, ta thấy thỏa mãn điều kiện.
`TH_{2}:x,y\ne=0`, ta có:
`3(x^2+xy+y^2)=x+8y`
`⇔3x^2+3xy+3y^2-x-8y=0`
`⇔6.(3x^2+3xy+3y^2-x-8y)=0.6`
`⇔18x^2+18xy+18y^2-6x-48y=0`
`⇔(9x^2+18xy+9y^2)+(9x^2-6x+1)+(9y^2-48y+64)-65=0`
`⇔(3x+3y)^2+(3x-1)^2+(3y-8)^2=65`
Vì `x,y ∈Z`, mà` (3x+3y)^2+(3x-1)^2+(3y-8)^2` là tổng của 3 số chính phương, nên ta có thể tách là:
`65=0+1+64=0+16+49=36+4+25.`
Ta thấy vai trò của số nguyên âm và số nguyên dương khi bình phương là như nhau nên ta giả sử `x\gey\ge0⇒(3x+3y)^2\ge(3x-1)^2\ge(3y-8)^2`, nên ta xét 3 trường hợp:
`1) 3x+3y=64, 3x-1=8, 3y-8=0 ⇒` không có `x,y` nguyên thỏa mãn.
`2) 3x+3y=49, 3x-1=16, 3y-8=0 ⇒` không có `x,y` nguyên thỏa mãn.
`3) 3x+3y=36, 3x-1=25, 3y-8=4 ⇒ x=y=1.`
Thử lại ta thấy thỏa mãn, vì ta đang giả sử cặp số tự nhiên nên cặp số nguyên âm ta cũng lấy vì khi bình phương đều cùng giá trị nên ta có thêm cặp `(x,y)` là `(-1;-1).`
Vậy ta có cặp `(x,y)` là `(0;0),(1;1),(-1;-1).`