function l = compute_lsf(N, sigma)

x = 1:N;
l = exp( -( (x-1-N/2).^2 )/(2*sigma*sigma) );

% normalize to unit area (DC gain = 1)
area = sum(sum(l));
l = l / area;