multinom

Computes coefficients powers for a multinomial expansion of the form (a_1+...+a_M)^N

Usage: [c, p] = multinom(M,N)

Computes coefficients 'c' & powers 'p' for a multinomial expansion such that

a_1 + ... + a_M)^N = \sum_{k=1}^K c(k) a_1^p(k,1) ... a_M^p(k,M),

where K = nchoosek(N+M-1,N).

'c' is a Kx1 vector of coeficients
'p' is a KxM matrix of powers, with k^th row containing the powers p(k,1) through p(k,M).