  // n=60, 70,71,72,73,  80, ...
x=linspace(-1,1,n);
y=x; 
for k=1:n,
   y(k)=1/(1+x(k)^2); 
 end;   // rand(x)
[ab,or]=TraceLagrange(x,y);
xbasc(),
plot2d(ab,or,2),
plot2d(x,y,-2,'000'),


//fin 