function junk(x)
  % From bandpass_filter:
  
  NP = N/2+1; % Number of distinct frequencies.
  fnyq = fsmp/2; % Nyquist frequency.
  dfnyq = NP-1;   % Discrete Nyquist frequncy.
  assert(fmin_lo <= fmin_hi)
  assert(fmax_lo <= fmax_hi)
  
  % Compute low and high frequencies of each shoulder:
  dfmin_lo = max([0,floor(fmin_lo*N/fsmp)])
  dfmin_hi = max([0,floor(fmin_hi*N/fsmp)])
  
  dfmax_lo = min([dfnyq, ceil(fmax_lo*N/fsmp)])
  dfmax_hi = min([dfnyq, ceil(fmax_hi*N/fsmp)])
  
  Y = ones(N,1);
  for k = 1:NP
    df = k - 1; % Discrete frequency
    % Low-end window weight:
    if dfmin_hi > 0
      vlo = shoulder(df,dfmin_lo,dfmin_hi);
    else
      vlo = 1.0;
    end
    if dfmax_lo < dfnyq
      vhi = shoulder(df,dfmax_hi,dfmax_lo);
    else
      vhi = 1.0;
    end
    Y(k,1)  = >SAME< * vlo*vhi;
    if k > 1
      Y(N+2-k,1)  = >SAME< * vlo*vhi;
    end
  end