Expectation of Stratonovich iterated integrals of Wiener processes

The solution of a (stochastic) differential equation (SDE) can be locally approximated by a stochastic expansion, a linear combination of iterated integrals. Quantities of interest, like moments, can then be approximated with the expansion. We present a formula for the case where the drivers of the equation are time and Wiener processes. We also present a Mathematica implementation of the result.