samer@4: function w=gauss_win(n,sigma)
samer@4: % gauss_win - Gaussian window of given length and std deviation
samer@4: %
samer@4: % gauss_win :: N:natural, nonneg~'std dev as multiple of N/2' -> [[N]].
samer@4: 
samer@4: s=(1:n)'-(n+1)/2;
samer@16: w=expquadf(2*s,sigma*n);