clear all;

% Grado massimo
Nmax = 20;

% Funzione da interpolare
f = @(x) cos(x) - sin(x);

% Punti di valutazione del polinomio
s = linspace(-1,1,500);

% Inizializzo i vettori degli errori
errE = zeros(1,Nmax);
errC = zeros(1,Nmax);

for n = 1:Nmax

    % Genero n+1 punti equidistanti da usare come nodi
    xE = linspace(-1,1,n+1);
    xC = -cos([0:n]*pi/n);

    % Valutazione della funzione sui nodi
    yE = f(xE);
    yC = f(xC);

    % Valutazione del polinomio nei punti s attraverso la funzione interpol
    tE = interpol_bary(xE,yE,s')';
    tC = interpol_bary(xC,yC,s')';

    % Errore massimo di interpolazione
    errE(n) = max(abs(tE-f(s)));
    errC(n) = max(abs(tC-f(s)));

end

figure(1)
semilogy(1:Nmax,errE,'b')
hold on;
semilogy(1:Nmax,errC,'m')
title('Errori');
legend('Pol Equi', 'Pol Cheby');