function [x,y]=chebrunge(n)
f=inline('1/(1+25*x^2)');
x=zeros(n+1,1);
y=zeros(n+1,1);
for i=0:n
    x(i+1)=cos(((2*i+1)/(n+1))*(pi/2));
    y(i+1)=f(x(i+1));
end;


