tel/TPmatlab/integration/TP3e/base_pn.m

   1 function[base]=base_pn(x,X)\r
   2 % engendre la base de Pn utilisée dans l'écriture de Newton\r
   3 % à partir des symboliques contenus dans x.\r
   4 %\r
   5 % variables d'entrée\r
   6 % x est un vecteur de deg_max+1 symboliques repréentant les points\r
   7 % du support d'interpolation.\r
   8 % X représente la variable symbolique\r
   9 % générique dans les écritures polynomiales.\r
  10 %\r
  11 % variables de sortie\r
  12 % base est un vecteur de symboliques représentant la base\r
  13 % 1,(X-x0),...,(X-x0)...(X-xn-1) de Pn.\r
  14 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\r
  15 %                       Fonctions connexes appelées\r
  16 %\r
  17 %\r
  18 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\r
  19 \r
  20 % tests d'entrée\r
  21 if nargin~=2\r
  22     error('Nombre d''arguments incorrect.');\r
  23 end\r
  24 \r
  25 % corps d'algorithme proprement dit\r
  26 \r
  27 % évite de passer une globale\r
  28 deg_max=size(x,2)-1;\r
  29 \r
  30 base=sym(ones(1,deg_max+1));gen=X*ones(1,deg_max)-x(1:deg_max);\r
  31 for i=2:deg_max+1,\r
  32     base(i)=prod(gen(1:i-1));\r
  33 end\r
  34 \r
  35 \r
  36 \r