function x=IFTform(X,fs)
%   x=IFTform(X,fs)
% Computes a 2-sided rep for x(t)
% from the Fourier Transform rep X
    N=length(X);                    %this better be a power of 2!
    y=fs*real(ifft([X,0,conj(flipud(X(2:N)')')]));
    x=[y(1:N);[y(1),flipud(y(N+2:2*N)')']];