function U=InterpolTrigo1(f,p); 

% U=InterpolTrigo1(f,p); 
%
%Calcule la suite des coefficients du polynôme trigonométrique
%(balancé par rapport à l'origine) de degré 2^p interpolant les 
%valeurs d'une fonction f définie sur [0,1] aux 2^p points équidistribués 
%j/2^p, j=0,...,2^p-1. 

J=[0:2^p-1]';
X=f(J/2^p).*(-1).^J;  
U=fft(X,2^p)/2^p;

% ou bien aussi, à la place de la dernière ligne :  
% U= CooleyTukey(X,p)/2^p;  
