function [f] = ifct(a)

// Output variables initialisation (not found in input variables)
f=[];

// Display mode
mode(0);

// Display warning for floating point exception
ieee(1);

//x=-cos(pi/N*(0:N)'')
N = max(size(a));



F = mtlb_fft([a(1,:);[a(2:N,:);a(N-1:-1:2,:)]/2]);
f = F(1:N,:);

f = real(f(N:-1:1,:));
f = real(f);
endfunction
