dmultinorm function, see page 393.
dmultinorm(xval,yval,mu.vector,sigma.matrix)
xval
yval
mu.vector
sigma.matrix