tel/TPmatlab/integration/TP3d/recherche_poids.m

   1 function [w1,w2,w3]=recherche_poids(t1,t2,t3)\r
   2 % recherche les poids d'une formule de quadrature sur [-1,1]\r
   3 % exacte pour toute les fonctions polynomes de degré au plus deux;\r
   4 % Le traitement est effectué en symbolique.\r
   5 %\r
   6 % variables d'entrée\r
   7 % t1,t2,t3 désignent trois dates distinctes de [-1,1] en lesquelles\r
   8 % sont effectuées les mesures de la fonction à intégrer;\r
   9 % elles serviront de support pour la formule de quadrature.\r
  10 %\r
  11 % variables de sortie\r
  12 % w1,w2,w3 sont les poids associés qui rendent l'intégration exacte\r
  13 % pour toute fonction polynome de degré au plus deux. Voir documentation\r
  14 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\r
  15 %                       Fonctions connexes appelées\r
  16 %\r
  17 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\r
  18 \r
  19 % tests d'entrée\r
  20 if nargin~=3\r
  21     error('Nombre d''arguments incorrect.');\r
  22 end\r
  23 \r
  24 % corps d'algorithme proprement dit\r
  25 modele=[t1,t1^2]';\r
  26 mat=sym(ones(3));\r
  27 mat(2:3,1)=modele;\r
  28 mat(2:3,2)=subs(modele,t1,t2);mat(2:3,3)=subs(modele,t1,t3);\r
  29 sol=(mat\sym([2 0 2/3]'))';\r
  30 w1=sol(1);w2=sol(2);w3=sol(3);\r
  31 \r