compute.sigma.star function

Computes intermediate correlation matrix