X=poly(0,'x');

// points uniformément répartis

n=10;

x = (-1:2/n:1);

l = x;

// polynomes de base l_i

for i=1:n,
  l(i) = 1;
  for j=1:n,
    if (i<>j) then
      l(i)=l(i)*(X-x(j))/(x(i)-x(j));
    end   
  end
end

// calcul du polynome d'interpolation

p=0;

for i=1:n,
  p = p+(1/(1+x(i)^2))*l(i);
end


// affichage du polynome

p

// graphe du polynome

t = (-1.2:0.1:1.2);

plot2d(t,horner(p,t),1);
plot2d(t,1 ./ ( 1+t .* t),2);



