compute_marginal_cor function

Computes non-positive-semidefinite approximation of correlation matrix