construct_covmat function

Constructing a covariance matrix for a variable number of phenotypes