DSolve[Derivative[1][y][x] == ((1 + x*y[x])*(1 + x + x^2 + 2*x*y[x] + x^2*y[x] + x^2*y[x]^2))/x^5,y[x],x]