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