Isserlis' theorem

formula that allows one to compute higher-order moments of the multivariate normal distribution in terms of its covariance matrix