function y=wgauss(x,mi,sigma)

y=1/(sqrt(2*pi)*sigma)*exp(-((x-mi)^2)/(2*sigma^2));