Solution:
We seek solutions (x,y,z) which are in arithmetic progression. Let us put y−x=z−y=d>0 so that the equation reduces to the form
3y2+2d2=2d3
Thus we get 3y2=2(d−1)d2. We conclude that 2(d−1) is 3 times a square. This is satisfied if d−1=6n2 for some n. Thus d=6n2+1 and 3y2=d2⋅2(6n2) giving us y2=4d2n2. Thus we can take y=2dn=2n(6n2+1). From this we obtain x=y−d=(2n−1)(6n2+1),z=y+d=(2n+1)(6n2+1). It is easily verified that
(x,y,z)=((2n−1)(6n2+1),2n(6n2+1),(2n+1)(6n2+1))
is indeed a solution for a fixed n and this gives an infinite set of solutions as n varies over natural numbers.