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