function [c]=interp2(n,f)
% fornisce i coefficienti del polinomio di interpolazione
% della funzione f sull'intervallo -1,1 con n punti di chebichef
for k=1:20
   x(k)=cos(pi*(k-1)/19); 
end

gg=f(x);
V=ones(n,1);
for k=1:n-1
    V=[V,(x').^k];
end
c=inv(V)*gg';
end