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