By Shalosh B. Ekhad and Doron Zeilberger
Written: Feb. 20, 2022
Abstract: Using recurrences gotten from the Apagodu-Zeilberger Multivariate Almkvist-Zeilberger algorithm we present a linear-time, and constant-space, algorithm to compute the general mixed moments of the k-variate general normal distribution, with any covariance matrix, for any specific k. Besides their obvious importance in statistics, these numbers are also very significant in enumerative combinatorics, since they count in how many ways, in a species with k different genders, a bunch of individuals can all get married, keeping track of the different kinds of heterosexual marriages. We completely implement our algorithm (with an accompanying Maple package, MVNM.txt) for the bivariate and trivariate cases (and hence taking care of our own 2-sex society and a putative 3-sex society), but alas, the actual recurrences for larger k took too long for us to compute. We leave them as computational challenges.
1 c12 c13
c12 1 c23
c13 c23 1
for SYMBOLIC (i.e. general) c12,c13, c23, for
1 ≤ m1,m2,m3 ≤ 20
The input file generates the output file
The input file generates the output file