function y = fgauss(x,mu,sigma)
y = exp(-0.5*((x-mu)/sigma).^2)/(sqrt(2*pi)*sigma);
