CONSEJO:
$$\left(x-xy(x)\right)+\left(y(x)+x^2\right)y'(x)=0\Longleftrightarrow$$ $$\left(x^2+y(x)\right)y'(x)-y(x)x+x=0\Longleftrightarrow$$
Deje $\text{Q}(x,y)=x-xy$ y $\text{T}(x,y)=x^2+y$.
Esta no es una ecuación exacta, porque $\frac{\partial\text{Q}(x,y)}{\partial y}=-x\ne2x=\frac{\partial\text{T}(x,y)}{\partial x}$.
Encuentre un factor integrante $\mu(y)$ tal que sea exacto:
$$\mu(y)\text{Q}(x,y)+\mu\cdot\frac{\text{d}y(x)}{\text{d}x}\cdot\text{T}(x,y)=0$$
Esto significa que:
$$\frac{\partial}{\partial y}\left(\mu(y)\text{Q}(x,y)\right)=\frac{\partial}{\partial y}\left(\mu(y)\text{T}(x,y)\right)$$
$$\frac{\text{d}\mu(y)}{\text{d}y}\left(x-yx\right)-\mu(y)x=2x\mu(y)\Longleftrightarrow$$ $$\frac{\mu'(y)}{\mu(y)}=-\frac{3}{y-1}\Longleftrightarrow$$ $$\int\frac{\mu'(y)}{\mu(y)}\space\text{d}y=\int-\frac{3}{y-1}\space\text{d}y\Longleftrightarrow$$ $$\ln\left(\mu(y)\right)=-3\ln\left(y-1\right)\Longleftrightarrow$$ $$\mu(y)=\frac{1}{(y-1)^3}$$
Deje $\text{P}(x,y)=-\frac{x}{(y-1)^2}$ y $\text{R}(x,y)=\frac{x^2+y}{(y-1)^3}$.
Esta es una ecuación exacta, porque $\frac{\partial\text{P}(x,y)}{\partial y}=\frac{2x}{(y-1)^3}=\frac{\partial\text{R}(x,y)}{\partial x}$.
Defina $f(x,y)$ tal que $\frac{\partial f(x,y)}{\partial x}=\text{P}(x,y)$ y $\frac{\partial f(x,y)}{\partial y}=\text{R}(x,y)$:
Entonces, la solución se dará por $f(x,y)=\text{C}$, donde $\text{C}$ es una constante arbitraria.
Integre $\frac{\partial f(x,y)}{\partial x}$ para encontrar $f(x,y)$ y donde $g(y)$ es una función arbitraria de $y$:
$$f(x,y)=\int-\frac{x}{(y-1)^2}\space\text{d}x=-\frac{x^2}{2(y-1)^2}+g(y)$$
Ahora, diferencie $f(x,y)$ para encontrar $g(y)$:
$$g'(y)=\frac{y}{(y-1)^3}$$
Entonces, al final tienes que resolver:
$$-\frac{x^2}{2(y-1)^2}+\frac{1-2y}{2(y-1)^2}=\text{C}$$