dsolve({diff(x(t),t) = -y(t)+x(t)*(x(t)^2+y(t)^2-1), diff(y(t),t) = x(t)+y(t)*(x(t)^2+y(t)^2-1)})