#GENERATE MULTIVARIATE NORMAL SAMPLE---------------------------
simulate.multnorm _ function(mu0,Q,M=1)
{
  s _ length(mu0)
  Z _  matrix(rnorm(s*M),nrow=s)
  V <- t(chol(Q) )
  x _ (V%*%Z)+mu0
  return(x)
}