Magnus and Neudecker's elimination matrix
[d]=elimination(n)
n = the size of the underlying var_cov matrix
d = the elimination matrix
e=elimination(3)