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