samer@4: function w=hgauss_win(n,sigma)
samer@4: % w=hgauss_win(n,sigma): half-Gaussian window 
samer@4: % 	half a gaussian window of length n
samer@4: %	sigma is std dev rel to window length
samer@4: 
samer@16: w=expquadf((0:n)',sigma*n);