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

