K. Phillips, “R Functions to Symbolically Compute the Central Moments of the Multivariate Normal Distribution”, J. Stat. Soft., vol. 33, no. 1, pp. 1–14, Feb. 2010.