clear all;
clf;

f=@(x) 1./(1+25*x.^2);
xx=linspace(-1,1);
plot(xx,f(xx),';target fun;')
hold on;
n=20;
pts1=linspace(-1,1,n);
p1=polyfit(pts1,f(pts1),n-1);
plot(xx,polyval(p1,xx),';equidistant;');
pts2=cos(pi*(2*(1:n)-1)/(2*n));
p2=polyfit(pts2,f(pts2),n-1);
plot(xx,polyval(p2,xx),';Csebisev;');

axis([-1 1 -2 2]);
