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