Called by callmultmoments to compute representation of a multivariate normal moment using recursive algorithm
multmoments(moment,current.matrix,current.cell,moment.rep,row_col)
vector c(k1,...,kn) specifying the moment X1**k1 *...* Xn**kn
upper-triangular integer matrix under consideration in recursion
cell in current matrix under consideration in recursion
current set of representations; mult.moments adds each satisfying matrix to moment.rep
matrix giving rows and columns for square matrix for each cell
moment representation, moment.rep, augmented with additional representations
Each row of the representation gives the exponents for a single product of covariance terms. For example, (1,2,0) represents S11**1 S12**2 S22**0, where the Sij are the covariances.
This function would normally only be called by callmultmoments.
K Phillips, Symbolic Computation of the Central Moments of the Multivariate Normal Distribution, Journal of Statistical Software, 2010.
callmultmoments (symmoments)