Sigma = {
if (klcloseparam){
return [[sigma1**2,rho*sigma1*sigma2],[rho*sigma1*sigma2,sigma2**2]]
}else{
return [[paramlogitn.sigma1**2, paramlogitn.rho*paramlogitn.sigma1*paramlogitn.sigma2],
[paramlogitn.rho*paramlogitn.sigma1*paramlogitn.sigma2, paramlogitn.sigma2**2]]
}
}