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