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