compute_canonical_covs function

Compute canonical covariance matrices.