Factor the left-hand side to obtain (x√a−y√b)(x√a+y√b)=1 Cubing both sides yields ((x3a+3xy2b)√a−(3x2ya+y3b)√b)((x3+3xy2)√a+(3x2ya+y3b)√b=1. Multiplying out, we obtain a(x3a+3xy2b)2−b(3x2ya+y3b)2=1 Therefore, if (x1,y1) is a solution of the equation, so is (x2,y2) with x2=x31a+3x1y21 and y2=3x21y1a+y31b. Clearly ,x2>x1 and y2>y1.Continuing in this way we obtain infinitely many solutions in positive integers.