function [F,X,Xc]=Func_TF(t,y)

// S.Ayrinhac (2012)
// Transformée de Fourier discrète utilisant la TF rapide (FFT)
// Xc : coefficients de Fourier complexes

// Tests
if size(t,2)==1 // on ne veut que des vecteurs lignes
   t=t';
end
if size(y,2)==1 // on ne veut que des vecteurs lignes
   y=y';
end

if length(t)~=length(y)
    disp('Attention ! Les dimensions des vecteurs en entrée ne sont pas identiques.');
end

N=length(t);

// Création du vecteur temps
deltat=abs(t(2)-t(1));   // intervalle de temps 
t=(0:1:N-1)*deltat;     // on recrée le vecteur temps
tmax=max(t);             // temps maximum

// Création du vecteur fréquences
deltanu=1/tmax;
F=deltanu.*(0:1:N-1);

// Création du vecteur amplitude de Fourier
Xc=fft(y)./(N/2); // coefficients de Fourier complexes (pour Func_TFI.sci)
              // le facteur 1/N : normalisation à l'amplitude
X=abs(Xc); // valeur absolue des coeffs de Fourier complexes
//X=fftshift(X); // symétrise le spectre autour de 0


    
endfunction