DifferentialRootReduce[Exp[1/(x^3 + 1)], x]
DifferentialRootReduce[Exp[1/(x^3 + 1)], x]; Integrate[%, x]
DifferentialRootReduce[Exp[1/(x^3 + 1)], x]; Integrate[%, x]; D[%, x] == Exp[1/(x^3 + 1)] // FullSimplify