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