function [F,X]=Func_TFD(t,y)

// Transformée de Fourier discrète
// voir commande dft()
y=y';

N=length(y);

deltat=abs(t(2)-t(1));  // intervalle Te
tmax=N*deltat; 
deltanu=1/tmax;
F=deltanu.*(0:N-1);



k=0:N-1;
for p=1:N
   X(p)=sum(y.*exp(-2.*%pi.*%i.*p.*k./N));
end
X=abs(X)/N;


endfunction