dsolve(diff(y(x),x) = (x-y(x))^3*(y(x)+x)^3*x/(-y(x)^2+x^2-1)/y(x),y(x))