function Y = gauss(X,mu,sigma)
Y = exp(-0.5*((X-mu)/sigma).^2)./sqrt(2*pi*sigma^2);
% Y = Y/exp(-0.5*(mu/sigma).^2)./sqrt(2*pi*sigma^2);