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