function F=interpolaABSchev(n)
% Primeramente se definen los puntos donde se va a interpolar,que seran los
% ceros de los polinomios de Tchebychev. Como queremos n+1 puntos tenemos
% que utilizar el polinomio de grado n+1.
m=n+1;
X=zeros(n+1,1);
X(1)=cos(pi/(2*m));
for i=1:n
    X(i+1)= cos(pi*(1+2*i)/(2*m));
end
Y=zeros(n+1,1);
for j=1:n+1
    Y(j)= abs(X(j));
end
% Se halla el polinomio de grado n que interpola los n+1 puntos
P=zeros(n+1,1);
P=polyfit (X,Y,n);
% Construyo una malla densa para graficar. 10000 puntos andaran OK.
W=zeros(10001,1);
W(1)=-1;
for i=1:10000
    W(i+1)= -1 + 2*i/10000;
end
% Evaluo el polinomio en la malla
Q=zeros(10001,1);
for i=1:10001
    Q(i)=polyval (P,W(i));
end
% Evaluo la funcion en la malla.
Z=zeros(10001,1);
for j=1:10001
    Z(j)= abs(W(j));
end
% Graficamos todo junto para comparar.
close
hold on
plot (W,Q,'b-')
plot (W,Z,'r-' )
hold off