The given equation can be written as
8x3+y3−6xy=4⇔(2x)3+y3+13−3⋅2x⋅y⋅1=5⇔(2x+y+1)(4x2+y2+1−2xy−2x−y)=5(1)
⇔21(2x+y+1)[(2x−y)2+(2x−1)2+(y−1)2]=5(2)
From (2), since (2x−y)2+(2x−1)2+(y−1)2>0, we have 2x+y+1>0, and hence from (1) we get
2x+y+1=1 or 2x+y+1=5⇔2x+y=0 or 2x+y=4.
From (1), for 2x+y=4, we have
4x2+y2+1−2xy−2x−y=1⇔(2x+y)2−6xy−(2x+y)=0⇔xy=2
From the system 2x+y=4, xy=2 we get the solution (x,y)=(1,2).
Also, from (1) for 2x+y=0 we have
4x2+y2+1−2xy−2x−y=5⇔(2x+y)2−6xy−(2x+y)=4⇔xy=−32.
However, from the system 2x+y=0, xy=−32 we get no solutions.
We also can work in the following way: The equation can be written as
8x3+y3−6xy=4⇔(2x+y)3−3⋅2x⋅y⋅(2x+y)−6xy=4, and so by putting 2x+y=s, 2xy=p, we get
s3−3ps−3p=4⇒p=3(s+1)s3−4.
Since 3p∈Z, we have that
s+1s3−4=s+1s3+1−5=s2−s+1−s+15∈Z.
---
Thus s+1 must be a divisor of 5, that is
s+1∈{−1,1,−5,5} or s∈{−2,0,−6,4},
And hence we find the pairs
(s,p)=(−2,4),(s,p)=(0,−34),(s,p)=(−6,344),(s,p)=(4,4),
From which only the last gives integer values for x,y, i. e. x=1,y=2.