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