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