function gauss_value = gauss2(mu,cov,x_n)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% function gauss_value = gauss2(mu,cov,x_n)          %%
%%						      %%
%% mu = mean (column vector)			      %%	
%% cov = covariance matrix			      %%		
%% x_n = data = column vector			      %%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

dim=size(mu,1);
gauss_value=1/((2*pi)^(dim/2)*sqrt(det(cov)))...
   *exp(-0.5*(x_n-mu)'*inv(cov)*(x_n-mu));