covW function

calculate the pooled within groups covariance matrix