dan@0: function [is, ie] = zpzh_idxs(N)
dan@0: md2 = mod(N,2);
dan@0: if md2 == 1
dan@0:   is =  (N + md2)/2:N;
dan@0:   ie = 1:(N + md2 - 1)/2;
dan@0: else
dan@0:   is =  N/2+1:N;
dan@0:   ie = [1:N /2];
dan@0: end
dan@0: 
dan@0: end