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