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