% tchebfft.m
%
%   coefficients de Tchebycheff d'une fonction de x, -1<= x <= 1
%
nomf=input('donnez la fonction, par exemple: 3*x+4, etc. N.B. x doit pouvoir etre un vecteur    ','s');
n=1;while n>0,
 n=input('nombre de subdivisions?  (stop si <=0  ');
 if n<=0,
   break;
 end;
 % valeurs de f:
 clear x;x=cos( pi*(0:n)/n );
 clear f;f=eval(nomf);
 % passage a 2*n
 clear g;g=[f(n+1:-1:1) f(2:n)];
 coef=cos(pi*(0:2*n-1)) .* fft(g) /n;
 coef(1:n+1),
end
