tel/TPmatlab/integration/TP3i/demo_stocke_points_poids.m

   1 % le script demo_stocke_points_poids permet de stocker sur le fichier pointspoids_comp\r
   2 % les différents points et de poids de gauss (pour les quatre méthodes) qui sont : \r
   3 %  - calculés manuellement (voir TP3b)\r
   4 %  - calculés de façon symbolique (voir TP3g)\r
   5 %  - calculés par diagonalisation (voir TP3h). \r
   6 %\r
   7 % \r
   8 % ************ Fonctions auxiliaires utilisées ************\r
   9 %\r
  10 %  points_poids_gauss, points_poids_gauss_symb, points_poids_gauss_diago\r
  11 %\r
  12 % *********************************************************\r
  13 \r
  14 \r
  15 clear     nAmax nBmax nCmax ...\r
  16           pointsleg poidsleg ...\r
  17           pointslegsymb poidslegsymb ...\r
  18           pointslegdiag poidslegdiag ...\r
  19           pointstch poidstch ...\r
  20           pointstchsymb poidstchsymb ...\r
  21           pointstchdiag poidstchdiag ...\r
  22           pointsher poidsher ...\r
  23           pointshersymb poidshersymb ...\r
  24           pointsherdiag poidsherdiag ...\r
  25           pointslag poidslag ...\r
  26           pointslagsymb poidslagsymb ...\r
  27           pointslagdiag poidslagdiag ...\r
  28 \r
  29 disp('méthode directe : ');\r
  30 nAmax=input('entrez l''entier maximal : ') ;\r
  31 disp('méthode symbolique : ');\r
  32 nBmax=input('entrez l''entier maximal : ') ;\r
  33 disp('méthode par diagonalisation : ');\r
  34 nCmax=input('entrez l''entier maximal : ') ;\r
  35 \r
  36 % calcul des poids et points de Gauss-Legendre.\r
  37 % calcul manuel : \r
  38 [pointsleg,poidsleg]=points_poids_gauss(1,nAmax);\r
  39 % calcul symbolique : \r
  40 [pointslegsymb,poidslegsymb]=points_poids_gauss_symb(1,nBmax);\r
  41 % calcul par diagonalisation : \r
  42 [pointslegdiag,poidslegdiag]=points_poids_gauss_diago(1,nCmax);\r
  43 \r
  44 % calcul des poids et points de Gauss-Tchebychev.\r
  45 % calcul manuel : \r
  46 [pointstch,poidstch]=points_poids_gauss(2,nAmax);\r
  47 % calcul symbolique : \r
  48 [pointstchsymb,poidstchsymb]=points_poids_gauss_symb(2,nBmax);\r
  49 % calcul par diagonalisation : \r
  50 [pointstchdiag,poidstchdiag]=points_poids_gauss_diago(2,nCmax);\r
  51 \r
  52 % calcul des poids et points de Gauss-Hermite.\r
  53 % calcul manuel : \r
  54 [pointsher,poidsher]=points_poids_gauss(3,nAmax);\r
  55 % calcul symbolique : \r
  56 [pointshersymb,poidshersymb]=points_poids_gauss_symb(3,nBmax);\r
  57 % calcul par diagonalisation : \r
  58 [pointsherdiag,poidsherdiag]=points_poids_gauss_diago(3,nCmax);\r
  59 \r
  60 % calcul des poids et points de Gauss-Laguerre.\r
  61 % calcul manuel : \r
  62 [pointslag,poidslag]=points_poids_gauss(4,nAmax);\r
  63 % calcul symbolique : \r
  64 [pointslagsymb,poidslagsymb]=points_poids_gauss_symb(4,nBmax);\r
  65 % calcul par diagonalisation : \r
  66 [pointslagdiag,poidslagdiag]=points_poids_gauss_diago(4,nCmax);\r
  67 \r
  68 %Stockage sur fichier\r
  69 save pointspoids_comp ...\r
  70           nAmax nBmax nCmax ...\r
  71           pointsleg poidsleg ...\r
  72           pointslegsymb poidslegsymb ...\r
  73           pointslegdiag poidslegdiag ...\r
  74           pointstch poidstch ...\r
  75           pointstchsymb poidstchsymb ...\r
  76           pointstchdiag poidstchdiag ...\r
  77           pointsher poidsher ...\r
  78           pointshersymb poidshersymb ...\r
  79           pointsherdiag poidsherdiag ...\r
  80           pointslag poidslag ...\r
  81           pointslagsymb poidslagsymb ...\r
  82           pointslagdiag poidslagdiag ...\r
  83 \r
  84 \r