Return the duplication matrix Dn which is the unique n^2 by n*(n+1)/2 matrix such that Dn vech (A) = vec (A) for all symmetric n by n matrices A.
See Magnus and Neudecker (1988), Matrix Differential Calculus with Applications in Statistics and Econometrics.
Package: octave