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