dsolve(diff(y(x),x) = (-3*x^2*y(x)-2*x^3-2*x-x*y(x)^2-y(x)+x^3*y(x)^3+3*x^4*y(x)^2+3*x^5*y(x)+x^6)/x/(x*y(x)+x^2+1),y(x))