% Calculate wavenumbers in radians per pixel for an FFT array
% of size N.
function k = wavenumber(N)

dk = 2*pi/N; % The stepsize for wavenumber is also the fundamental 
   % freqency on the domain.

k = dk*[0:floor(N/2), -floor((N-1)/2):-1];



