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