construct_covmat_multi function

Constructing a covariance matrix for multiple phenotypes