
%                     N
%       X(k) =       sum  x(n)*exp(-j*2*pi*(k-1)*(n-1)/N), 1 <= k <= N.
%                    n=1

clear all;
format long;



% f=inline('2*exp(-i*1*x)+4*exp(-i*5*x)+8*exp(i*3*x)');
%f=inline('2*exp(-i*0*x)+2*exp(i*20*x)+2*exp(i*19*x)');
f=inline('exp(i*2*t)');

% KENDALL.
n=8;
t=(0:2*n)'*(2*pi/(2*n+1));
ft=feval(f,t);

t_start=cputime;
for l=-n:n
    exp_KA=exp(-i*l*t);
    c(l+n+1,1)=(exp_KA'*ft)/(2*n+1);
end
t_end=cputime;
tempo(1)=t_end-t_start;

% FFT.

n=2*n;
h=2*pi/n;
t=(0:n-1)'*h;
ft=feval(f,t);

t_start=cputime;
c_fft=fft(ft,n)*h/(2*pi);
t_end=cputime;
tempo(2)=t_end-t_start;

format long;
tempo