function showfft ( )
%
%% showfft() compares the ordinary DFT with the FFT.
%
%  For various reasons, this code is OBSOLETE.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'showfft():\n' );

x = randn(1024,1)+sqrt(-1)*randn(1024,1);
disp('   n       DFT Flops    FFTRecur Flops    FFT Flops');
disp('------------------------------------------------------')
for n = [2 4 8 16 32 64 128 256]
   tic();
   y = DFT(x(1:n));
   dftFlops = toc();
   tic();
   y = FFTRecur(x(1:n));
   recurFlops = toc();
   tic();
   y = fft(x(1:n));
   fftFlops = toc();  
   disp(sprintf(' %4.0f   %10.0f    %10.0f      %10.0f', n,dftFlops,recurFlops,fftFlops))
end

  return
end

