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