dgmrf2 function

Reparametrised multivariate Gaussian distribution