clear
n=input('grado polinomio = ');
a=-5; 
b=5;

% nodi equispaziati
h=(b-a)/n;
xi=a:h:b;
yi=runge(xi);

% nodi di Chebyshev
xic=(a+b)/2+(b-a)/2*cos([0:n]*pi/n);
yic=runge(xic);

% diff. divise su nodi equispaziati
de=diffDivise(xi,yi);
% diff. divise su nodi di Chebyshev
dc=diffDivise(xic,yic);

x=linspace(a,b,100);
g=runge(x);

for i=1:length(x),
    pe(i)=horner(xi,x(i),de); %equispaziati
    pc(i)=horner(xic,x(i),dc); %chebyshev
end

subplot(1,2,1)
plot(xi,yi,'o',x,pe,'b-',x,g,'r-');
legend('nodi equisp.','polinomio','funzione')

subplot(1,2,2)
plot(xic,yic,'o',x,pc,'b-',x,g,'r-');
legend('nodi di Cheby.','polinomio','funzione')
