function k=wavenumber(j,n)
% WAVENUMBER - Return the wavenumber corresponding to index j
  if j<=n/2+1
     k=j-1;
  else
     k=j-1-n;
  end