tel/TPmatlab/integration/TP3g/polyval_symb.m

   1 function res=polyval_symb(P,x)\r
   2 \r
   3 %       polyval_symb : évaluation formelle de polynôme\r
   4 %\r
   5 % *********************************************************\r
   6 %  ATTENTION, ni polyval ni poly2sym ne marchent sur un tableau \r
   7 %             à coefficients formels !!\r
   8 %\r
   9 % *********************************************************\r
  10 %\r
  11 %\r
  12 %       res=polyval_sym(P,x) : renvoie la valeur de P en x (en formel) :\r
  13 %  si n est la longueur de P : \r
  14 %     res=P(1)*x^(n-1)*+P(2)*x^(n-2)+...+P(n-1)*x+P(n).\r
  15 %\r
  16 %       variables d'entrées : \r
  17 %    *  P : polynoômes à coefficients formels\r
  18 %    *  x : réel ou symbolique.\r
  19 %       variables de sortie : \r
  20 %    * res : valeur de P en x\r
  21 %\r
  22 %\r
  23 %\r
  24 % ************ Fonctions auxiliaires utilisées ************\r
  25 %\r
  26 %       aucune\r
  27 %\r
  28 % *********************************************************\r
  29 %\r
  30 \r
  31 % Contrôles d'entrée\r
  32 % nombre d'arguments\r
  33 if nargin~=2\r
  34    error('nombre d''arguments de la fonction incorrect');\r
  35 end\r
  36 \r
  37 \r
  38 % Corps d'algorithme.\r
  39 n=length(P);\r
  40 if (isempty(P))\r
  41     res=0;\r
  42 else\r
  43    if x==0 \r
  44       res=P(n);\r
  45    else\r
  46       res=sum(P.*x.^(n-1:-1:0));\r
  47   end\r
  48 end