function p=interpolazione(f,x)
%ritorna i coefficienti del polinomio che
%interpola f sui punti x(1),...x(end)

n=length(x);

%valuta f
y=zeros(n,1);
for i=1:n
    y(i)=f(x(i));
end

%genera la matrice di Vandermonde
A=zeros(n,n);
for i=1:n
    for j=1:n
        A(i,j)=x(i)^(n-j);
    end
end

cond(A)

p=A\y;

