function [c]=interp1(n,a,b,f)
% fornisce i coefficienti del polinomio di interpolazione
% della funzione f sull'intervallo a,b con n punti equispaziati
x=a:(b-a)/(n-1):b;
gg=f(x);
V=ones(n,1);
for k=1:n-1
    V=[V,(x').^k];
end
c=inv(V)*gg';
end