function y=myIFFT(x)
% MYIFFT inverse fast Fourier transform
%   y=myIFFT(x); computes recursively the inverse Fourier tranform of
%   the vector x whose length must be a power of 2.

n=length(x);
if n==1,
  y=x;
else
  w=exp(2i*pi/n*(0:n/2-1)');
  ze=myIFFT(x(1:2:n-1));
  zo=myIFFT(x(2:2:n));
  y=[ze+w.*zo; ze-w.*zo];
end;