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