function ret = gaussian(x,mu,sigma2),

ret = 1/sqrt(2*pi*sigma2)*exp(-(x-mu).^2/(2*sigma2));

