
function g = Gaussian(x,mu,sigma)
 % Return the value of N(x;mu,sigma)

 D=length(x);
 s=0;
 for i=1:D,
   s = s + ((x(i)-mu(i))/sigma(i))^2;
end
 
g =  exp(-0.5*s);
