tel/TPmatlab/erreur_algorithmique/TP1d/sommeB.m

   1 function res=sommeB(p,n)\r
   2 \r
   3 %       sommeB : calcul de somme (version B).\r
   4 %\r
   5 % *********************************************************\r
   6 %\r
   7 %       res=sommeB(p,n)\r
   8 %       \r
   9 %  calcul de la somme \sum_{i=1}^n  1/(i^p) en utilisant un calcul \r
  10 %  matriciel et en sommant dans l'ordre des indices croissants.\r
  11 %\r
  12 %       variables d'entrées :\r
  13 %   * p,n : deux entiers \r
  14 %\r
  15 %       variables de sortie\r
  16 %   * res: résultats de la somme\r
  17 % \r
  18 % \r
  19 % ************ Fonctions auxiliaires utilisées ************\r
  20 %\r
  21 %       aucune\r
  22 %\r
  23 % *********************************************************\r
  24 \r
  25 \r
  26 \r
  27 % Contrôles d'entrée\r
  28 % nombre d'arguments\r
  29 if nargin~=2\r
  30    error('nombre d''arguments de la fonction incorrect');\r
  31 end\r
  32 % autres tests éventuels \r
  33    \r
  34 \r
  35 % Corps d'algorithme\r
  36 res=sum(1./((1:n).^p));\r