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