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