DSolve[Derivative[1][y][x] == (y[x]*(x + x^2 + x*y[x] + y[x]^2))/x^2,y[x],x]