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