mix_matrix function

Estimate mixture covariance matrix from within-group covariance matrices