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