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