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