function [f]=MVN(x,Ginv,factor)
  f = exp(-x'*Ginv*x/2)/ factor;