function f=tdf(g)
% transformada discreta de Fourier

N=length(g);
D=zeros(1,N);
tol=1e-12;
% cálculo da DFT
for ij=1:N
  j=ij-1;
  for ik=1:N
    k=ik-1;
    D(ij)=D(ij)+g(ik)*exp(-1i*j*2*pi*k/N);
  end
  if abs(D(ij)) < tol
    D(ij)=0;
  end
end
f=D;