%Funktion f, Intervallgrenzen a und b, Zahl der Tschebyschef-Stuetzstellen n
function ak=interpol(f, a, b, n)

D=zeros(1,4*n);

for j=0:n-1
  D(2*(j+1))=f( ((b-a)/2) * cos((2*j+1)*pi/(2*n)) + (a+b)/2 );
  D(2*(2*n-j))=D(2*(j+1));
endfor
aTilde=4*ffTU(D);

ak= zeros(1,n);
ak(1)=0.5*aTilde(1);
for j=2:n
  ak(j)=aTilde(j);
endfor

sn=@(x) P(ak,x, a, b);
%fplot(sn, [a-0.1*(b-a),b+0.1*(b-a)], 500)
fplot(sn, [a,b], 500)
endfunction

function y=P(akk, x, a, b)
n=size(akk,2);
y=0;
for j=0:(n-1)
  y=y+akk(j+1)*cos(j*acos(2*x/(b-a)-(b+a)/(b-a)));
endfor
endfunction
