y'(x)=\frac{1}{x y(x) \left(x y(x)^2+x+1\right)}