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)