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