DSolve[-(x*y[x]) + y[x]^2 + (x^2 + x*y[x] + 2*x*y[x]^3)*Derivative[1][y][x] == 0,y[x],x]