function z=kmvn(x)

% KMVN        Multivariate Normal Density Function.

[n d]=size(x);

s=diag(x*x');

z=(2*pi)^(-d/2)*exp(-0.5*s);           