DSolve[Derivative[1][y][x] == (-1/2 + (3*x*y[x])/4 + x*y[x]^2 - (3*x^2*y[x]^2)/8 + x*y[x]^3 - (x^2*y[x]^3)/2 + (x^3*y[x]^3)/16)/(x*(-2 - 2*y[x] + x*y[x])),y[x],x]