compute_sigma function

Compute covariance matrix for some reference-based methods (JR, CIR)