function w=gauss_win(n,sigma)
% gauss_win - Gaussian window of given length and std deviation
%
% gauss_win :: N:natural, nonneg~'std dev as multiple of N/2' -> [[N]].

s=(1:n)'-(n+1)/2;
w=expquadf(2*s,sigma*n);
