Asumiendo $x\ne0$$y\ne0$, dividiendo ( $1$ ) $x$ e ( $2$ ) $y$ da
\begin{equation}
\dfrac{x^\prime}{x}=\dfrac{y^\prime}{y}
\end{equation}
así
\begin{equation}
\dfrac{d}{dt}\ln(x)=\dfrac{d}{dt}\ln(y)
\end{equation}
así
\begin{equation}
x(t)=c_1y(t)
\end{equation}
Así que la ecuación ($2$) se convierte en
\begin{equation}
y^\prime-y=-(1+c_1)y^2
\end{equation}
por lo que es fácil de Bernoulli ecuación.
La solución de la ecuación de Bernoulli da
\begin{align}
x(t)&=\dfrac{c_1}{c_2e^{-t}+c_1+1}\\
y(t)&=\dfrac{1}{c_2e^{-t}+c_1+1}
\end{align}
Por lo $x(0)=\dfrac{c_1}{c_1+c_2+1}>0$$y(0)=\dfrac{1}{c_1+c_2+1}>0$.
Así conocemos $1+c_1+c_2>0$ $c_1>0$ y $x(t)>0,y(t)>0$ todos los $t$.
Procedimiento en el supuesto de que ni $x(t)$ ni $y(t)$ nunca $0$, obtenemos soluciones para que ese es el verdadero.
Esto deja abierta la cuestión de si hay alguna otra solución para el sistema en el que lo que es falso (aparte de los obvios $x(t)=y(t)=0$).
Tenga en cuenta que
\begin{equation}
y^\prime=f(x,y)=-y^2-xy+y
\end{equation}
es continua en un barrio de todos los $(t_0,y_0)$ y que
\begin{equation}
\dfrac{\partial f}{\partial y}=-2y-x+1
\end{equation}
también es continua en un barrio de todos los $(t_0,y_0)$, por lo que no hay una única solución para la ecuación ($2$) que contiene cada una de las $(t_0,y_0)$, en este caso en particular en $(0,y_0)$. El mismo argumento vale para la ecuación ($1$). No hay una única ecuación de ($1$) que contenga $(0,x_0)$.