Modify logGaussian fct to accept dense covariance matrices

Now COVAR is a square array corresponding to the covariance
matrix.
MEAN and X must be column vectors with dimension compatible with
the covariance matrix.
Define priors for alpha, sigma and H

Still missing information to define prior for Cd.