function f=spectre(S,T);

%                  f=spectre(S,T);
%
% Calcule, étant donné un signal analogique discrétisé t --> S(t) donné comme 
% un vecteur colonne avec
%                       t=-T:T/M:T-T/M 
% les valeurs du spectre approché de S (en colonne) aux points 
%                    omega=k pi/T, k=-M:1:M-1

N=length(S);
M=fix(N/2);
w=-ones(N,1);
for k=1:2:N
  w(k)=-w(k);
end
u=fft(S);
u=w.*u;
u1=u(1:M);
u2=u(M+1:N);
f=[u2;u1];
f=(2*T/N)*f;