Phillips, K. “R Functions to Symbolically Compute the Central Moments of the Multivariate Normal Distribution”. Journal of Statistical Software, Code Snippets, vol. 33, no. 1, Feb. 2010, pp. 1-14, doi:10.18637/jss.v033.c01.