function c = metodo3 (f,x)

n1 = length(x);
y = f(x);

% prealloca la tabella delle differenze divise
ctab = zeros(n1,n1);

% inizializza la prima colonna
ctab(:,1) = y;

% costruisce la tabella per colonne
for j = 2:n1
    for i = 1:n1-j+1
        ctab(i,j) = (ctab(i,j-1)-ctab(i+1,j-1))/(x(i)-x(i+j-1));
    end
end
ctab;
% definisce il vettore dei coefficienti
c = ctab(1,:);
c = c(:);