construct_covmat_single function

Constructing a covariance matrix for a single phenotype