function [c,x]=interpp2(n,f)
% INTERP,2: c=interp2,(n,f)
% fornisce i coefficienti del polinomio di interpolazione
% della funzione f sull'intervallo -1,1 con n punti di chebichef
% c(1) e' il coefficiente della potenza di grado massimo.
for k=1:n
   x(k)=cos(pi*(k-1)/(n-1)); 
end

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