function[base]=base_pn(x,X)
% engendre la base de Pn utilisée dans l'écriture de Newton
% à partir des symboliques contenus dans x.
%
% variables d'entrée
% x est un vecteur de deg_max+1 symboliques repréentant les points
% du support d'interpolation.
% X représente la variable symbolique
% générique dans les écritures polynomiales.
%
% variables de sortie
% base est un vecteur de symboliques représentant la base
% 1,(X-x0),...,(X-x0)...(X-xn-1) de Pn.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                       Fonctions connexes appelées
%
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% tests d'entrée
if nargin~=2
    error('Nombre d''arguments incorrect.');
end

% corps d'algorithme proprement dit

% évite de passer une globale
deg_max=size(x,2)-1;

base=sym(ones(1,deg_max+1));gen=X*ones(1,deg_max)-x(1:deg_max);
for i=2:deg_max+1,
    base(i)=prod(gen(1:i-1));
end



