function [c,x]=interpp1(n,a,b,f)
% INTERPP1   c=interpp1(n,a,b,f)
% fornisce i coefficienti del polinomio di interpolazione
% della funzione f sull'intervallo a,b con n punti equispaziati
% c(1) e' il coefficiente della potenza di grado piu' alto.
x=a:(b-a)/(n-1):b;
gg=f(x);
V=ones(n,1);
for k=1:n-1
    V=[(x').^k,V];
end
c=V\gg';
end
