Télécharger la plaquette Cast3M


fr en

Télécharger quadratu.procedur

Retour à la liste

Numérotation des lignes : 

  1. * QUADRATU  PROCEDUR  AM        99/07/13    21:15:08     3632           
  2. 'DEBPROC' QUADRATU LOI*'MOT' MU*'FLOTTANT' SIGMA*'FLOTTANT' N*'ENTIER';
  3. *|=====================================================================
  4. *|
  5. *| << OBJET >> :
  6. *|
  7. *| Procedure permettant de calculer les points et poids d'intégration
  8. *| associés à une variable aléatoire.
  9. *|
  10. *|
  11. *| << ENTREE >> :
  12. *| -- ARGUMENTS OBLIGATOIRES DANS TOUS LES CAS --
  13. *| 
  14. *| LOI = objet de type MOT designant la loi suivie par la 
  15. *|                variable aléatoire ('NORM', 'LOGN',
  16. *|                'EXPO', 'UNIF')
  17. *|
  18. *| MU = objet de type REEL designant la valeur moyenne
  19. *|                    de la variable aleatoire.
  20. *|
  21. *| SIGMA = objet de type REEL designant l'ecart-type
  22. *|                       de la variable aléatoire.
  23. *|
  24. *| N = objet de type ENTIER designant le nombre de points
  25. *|              d'integration a calculer.
  26. *|
  27. *|
  28. *|
  29. *| << SORTIE >> :
  30. *|
  31. *| RESTAB = objet de type TABLE dont les indices sont des objets de 
  32. *|          type ENTIER.
  33. *|
  34. *| RESTAB.i = objet de type TABLE dont les indices sont des entiers
  35. *| RESTAB.i.'POINT' = objet de type REEL contenant le i-eme point 
  36. *|                  d'integration
  37. *| RESTAB.i.'POIDS' = objet de type REEL contenant le i-eme poids 
  38. *|                  d'integration
  39. *|
  40. *|=====================================================================
  41. *                                                                               
  42. * ______________________                                                        
  43. *                                                                               
  44. * DECLARATIONS GENERALES                                                        
  45. * ______________________                                                        
  46. *                                                                               
  47. SI = 'MOT' 'SI' ; SINON = 'MOT' 'SINON' ; FINSI = 'MOT' 'FINSI' ;               
  48. REPETER = 'MOT' 'REPETER'; QUITTER = 'MOT' 'QUITTER'; FIN = 'MOT' 'FIN';        
  49. EXISTE = 'MOT' 'EXISTE' ; EGA = 'MOT' 'EGA' ;                                   
  50. ERREUR = 'MOT' 'ERREUR' ;                                                       
  51. ENTIER   = 'MOT' 'ENTIER  ' ;                                                   
  52. FLOTTANT = 'MOT' 'FLOTTANT' ;                                                   
  53. MOT      = 'MOT' 'MOT     ' ;                                                   
  54. * _____                                                        
  55. *                                                                               
  56. * TESTS                                                        
  57. * _____                                                        
  58. *                                                                               
  59. SI (EGA LOI 'EXPO');
  60.   SI (NEG MU SIGMA);
  61.     MESS 'ERREUR : VOUS AVEZ SPECIFIE UNE VALEUR MOYENNE DIFFERENTE ';
  62.     MESS '         DE L ECART-TYPE. ';
  63.     MESS '     !!! POUR UNE LOI EXPONENTIELLE : mu = sigma !!!';
  64.   QUITTER QUAD;
  65.   FINSI;
  66. FINSI;
  67. *
  68. SI ( N '>EG' 11 );
  69.   MESS 'ERREUR : LE NOMBRE DE POINTS D INTEGRATION EST LIMITE A 10';
  70.   QUITTER QUAD;
  71. FINSI;
  72. *
  73. SI (NEG LOI 'UNIF');
  74.   SI (NEG LOI 'EXPO');
  75.     SI (NEG LOI 'LOGN');
  76.       SI (NEG LOI 'NORM');
  77. MESS 'ERREUR : VOUS AVEZ SPECIFIE UNE LOI DE PROBABILITE INCONNUE';
  78. MESS '         LA LOI ' LOI ' N EST PAS DEFINIE ';
  79. QUITTER QUAD;
  80.       FINSI;
  81.     FINSI;
  82.   FINSI;
  83. FINSI;
  84. * _____________________________________                                         
  85. *                                                                               
  86. * POIDS ET POINTS D INTEGRATION TABULES                                         
  87. * _____________________________________                                         
  88. *                ______________
  89. *                
  90. *                 LOI NORMALE
  91. *                ______________
  92. *
  93. *
  94. VaNORM = TABLE;
  95. REPE BOUC 10;
  96.   VaNORM. &BOUC = TABLE;
  97.   VaNORM.  &BOUC .'POINT' = TABLE;
  98.   VaNORM.  &BOUC .'POIDS' = TABLE;
  99. FIN BOUC;
  100. *
  101. *
  102. VaNORM. 1 .'POINT'. 1 = 0        ;VaNORM. 1 .'POIDS'.  1 = 1        ;
  103. VaNORM. 2 .'POINT'.  1 = 1        ;VaNORM. 2 .'POIDS'.  1 = .5      ;
  104. VaNORM. 2 .'POINT'.  2 = -1        ;VaNORM. 2 .'POIDS'.  2 = .5      ;
  105. VaNORM. 3 .'POINT'.  1 = -1.732051 ;VaNORM. 3 .'POIDS'.  1 = .166667 ;
  106. VaNORM. 3 .'POINT'.  2 = 1.732051 ;VaNORM. 3 .'POIDS'.  2 = .166667 ;
  107. VaNORM. 3 .'POINT'.  3 = 0        ;VaNORM. 3 .'POIDS'.  3 = .666667 ;
  108. VaNORM. 4 .'POINT'.  1 = -.741964 ;VaNORM. 4 .'POIDS'.  1 = .454124 ;
  109. VaNORM. 4 .'POINT'.  2 = .741964 ;VaNORM. 4 .'POIDS'.  2 = .454124 ;
  110. VaNORM. 4 .'POINT'.  3 = -2.334414 ;VaNORM. 4 .'POIDS'.  3 = .045876 ;
  111. VaNORM. 4 .'POINT'.  4 = 2.334414 ;VaNORM. 4 .'POIDS'.  4 = .045876 ;
  112. VaNORM. 5 .'POINT'.  1 = -1.355626 ;VaNORM. 5 .'POIDS'.  1 = .222076 ;
  113. VaNORM. 5 .'POINT'.  2 = 1.355626 ;VaNORM. 5 .'POIDS'.  2 = .222076 ;
  114. VaNORM. 5 .'POINT'.  3 = -2.85697  ;VaNORM. 5 .'POIDS'.  3 = .011257 ;
  115. VaNORM. 5 .'POINT'.  4 = 2.85697  ;VaNORM. 5 .'POIDS'.  4 = .011257 ;
  116. VaNORM. 5 .'POINT'.  5 = 0        ;VaNORM. 5 .'POIDS'.  5 = .533333 ;
  117. VaNORM. 6 .'POINT'.  1 = -.616707 ;VaNORM. 6 .'POIDS'.  1 = .408828 ;
  118. VaNORM. 6 .'POINT'.  2 = .616707 ;VaNORM. 6 .'POIDS'.  2 = .408828 ;
  119. VaNORM. 6 .'POINT'.  3 = -1.889176 ;VaNORM. 6 .'POIDS'.  3 = .088616 ;
  120. VaNORM. 6 .'POINT'.  4 = 1.889176 ;VaNORM. 6 .'POIDS'.  4 = .088616 ;
  121. VaNORM. 6 .'POINT'.  5 = -3.324257 ;VaNORM. 6 .'POIDS'.  5 = .002556 ;
  122. VaNORM. 6 .'POINT'.  6 = 3.324257 ;VaNORM. 6 .'POIDS'.  6 = .002556 ;
  123. VaNORM. 7 .'POINT'.  1 = -1.154405 ;VaNORM. 7 .'POIDS'.  1 = .240123 ;
  124. VaNORM. 7 .'POINT'.  2 = 1.154405 ;VaNORM. 7 .'POIDS'.  2 = .240123 ;
  125. VaNORM. 7 .'POINT'.  3 = -2.366759 ;VaNORM. 7 .'POIDS'.  3 = .030757 ;
  126. VaNORM. 7 .'POINT'.  4 = 2.366759 ;VaNORM. 7 .'POIDS'.  4 = .030757 ;
  127. VaNORM. 7 .'POINT'.  5 = -3.75044  ;VaNORM. 7 .'POIDS'.  5 = .000548 ;
  128. VaNORM. 7 .'POINT'.  6 = 3.75044  ;VaNORM. 7 .'POIDS'.  6 = .000548 ;
  129. VaNORM. 7 .'POINT'.  7 = 0        ;VaNORM. 7 .'POIDS'.  7 = .457143 ;
  130. VaNORM. 8 .'POINT'.  1 = -.53908  ;VaNORM. 8 .'POIDS'.  1 = .373012 ;
  131. VaNORM. 8 .'POINT'.  2 = .53908  ;VaNORM. 8 .'POIDS'.  2 = .373012 ;
  132. VaNORM. 8 .'POINT'.  3 = -1.636519 ;VaNORM. 8 .'POIDS'.  3 = .11724  ;
  133. VaNORM. 8 .'POINT'.  4 = -2.802486 ;VaNORM. 8 .'POIDS'.  4 = .009635 ;
  134. VaNORM. 8 .'POINT'.  5 = 2.802486 ;VaNORM. 8 .'POIDS'.  5 = .009635 ;
  135. VaNORM. 8 .'POINT'.  6 = -4.144547 ;VaNORM. 8 .'POIDS'.  6 = .000113 ;
  136. VaNORM. 8 .'POINT'.  7 = 4.144547 ;VaNORM. 8 .'POIDS'.  7 = .000113 ;
  137. VaNORM. 8 .'POINT'.  8 = 1.636519 ;VaNORM. 8 .'POIDS'.  8 = .11724  ;
  138. VaNORM. 9 .'POINT'.  1 = 4.512746 ;
  139. VaNORM. 9 .'POIDS'.  1 = 2.234584e-05 ;
  140. VaNORM. 9 .'POINT'.  2 = 2.076848 ;VaNORM. 9 .'POIDS'.  2 = .049916 ;
  141. VaNORM. 9 .'POINT'.  3 = -3.205429 ;VaNORM. 9 .'POIDS'.  3 = .002789 ;
  142. VaNORM. 9 .'POINT'.  4 = 3.205429 ;VaNORM. 9 .'POIDS'.  4 = .002789 ;
  143. VaNORM. 9 .'POINT'.  5 = 1.023256 ;VaNORM. 9 .'POIDS'.  5 = .244098 ;
  144. VaNORM. 9 .'POINT'.  6 = -2.076848 ;VaNORM. 9 .'POIDS'.  6 = .049916 ;
  145. VaNORM. 9 .'POINT'.  7 = -4.512746 ;
  146. VaNORM. 9 .'POIDS'.  7 = 2.234584e-05 ;
  147. VaNORM. 9 .'POINT'.  8 = 0        ;VaNORM. 9 .'POIDS'.  8 = .406349 ;
  148. VaNORM. 9 .'POINT'.  9 = -1.023256 ;VaNORM. 9 .'POIDS'.  9 = .244098 ;
  149. VaNORM. 10 .'POINT'.  1 = -3.581823 ;VaNORM. 10 .'POIDS'.  1 = .000758 ;
  150. VaNORM. 10 .'POINT'.  2 = .484936 ;VaNORM. 10 .'POIDS'.  2 = .344642 ;
  151. VaNORM. 10 .'POINT'.  3 = -2.484326 ;VaNORM. 10 .'POIDS'.  3 = .019112 ;
  152. VaNORM. 10 .'POINT'.  4 = -.484936 ;VaNORM. 10 .'POIDS'.  4 = .344642 ;
  153. VaNORM. 10 .'POINT'.  5 = 1.465989 ;VaNORM. 10 .'POIDS'.  5 = .135484 ;
  154. VaNORM. 10 .'POINT'.  6 = 2.484326 ;VaNORM. 10 .'POIDS'.  6 = .019112 ;
  155. VaNORM. 10 .'POINT'.  7 = 3.581823 ;VaNORM. 10 .'POIDS'.  7 = .000758 ;
  156. VaNORM. 10 .'POINT'.  8 = -4.859463 ;
  157. VaNORM. 10 .'POIDS'.  8 = 4.310653e-06 ;
  158. VaNORM. 10 .'POINT'.  9 = 4.859463 ;
  159. VaNORM. 10 .'POIDS'.  9 = 4.310653e-06 ;
  160. VaNORM. 10 .'POINT'.  10 = -1.465989 ;
  161. VaNORM. 10 .'POIDS'.  10 = .135484 ;
  162. *                _________________
  163. *                
  164. *                 LOI LOG-NORMALE
  165. *                _________________
  166. *
  167. * Les points et poids d integration de la loi normale sont
  168. * utilises ici. (Changement de variable)
  169. *
  170. *                ___________________
  171. *                
  172. *                 LOI EXPONENTIELLE
  173. *                ___________________
  174. VaEXPO = TABLE;
  175. REPE BOUC 10;
  176.   VaEXPO. &BOUC = TABLE;
  177.   VaEXPO. &BOUC .'POINT' = TABLE;
  178.   VaEXPO. &BOUC .'POIDS' = TABLE;
  179. FIN BOUC;
  180. VaEXPO. 1 .'POINT'. 1 = 1        ;VaEXPO. 1 .'POIDS'. 1 = 1        ;
  181. VaEXPO. 2 .'POINT'. 1 = .585786 ;VaEXPO. 2 .'POIDS'. 1 = .853553 ;
  182. VaEXPO. 2 .'POINT'. 2 = 3.414214 ;VaEXPO. 2 .'POIDS'. 2 = .146447 ;
  183. VaEXPO. 3 .'POINT'. 1 = .415775 ;VaEXPO. 3 .'POIDS'. 1 = .711093 ;
  184. VaEXPO. 3 .'POINT'. 2 = 2.29428  ;VaEXPO. 3 .'POIDS'. 2 = .278518 ;
  185. VaEXPO. 3 .'POINT'. 3 = 6.289945 ;VaEXPO. 3 .'POIDS'. 3 = .010389 ;
  186. VaEXPO. 4 .'POINT'. 1 = 9.395071 ;VaEXPO. 4 .'POIDS'. 1 = .000539 ;
  187. VaEXPO. 4 .'POINT'. 2 = 4.53662  ;VaEXPO. 4 .'POIDS'. 2 = .038888 ;
  188. VaEXPO. 4 .'POINT'. 3 = .322548 ;VaEXPO. 4 .'POIDS'. 3 = .603154 ;
  189. VaEXPO. 4 .'POINT'. 4 = 1.745761 ;VaEXPO. 4 .'POIDS'. 4 = .357419 ;
  190. VaEXPO. 5 .'POINT'. 1 = 1.413403 ;VaEXPO. 5 .'POIDS'. 1 = .398667 ;
  191. VaEXPO. 5 .'POINT'. 2 = 3.596426 ;VaEXPO. 5 .'POIDS'. 2 = .075942 ;
  192. VaEXPO. 5 .'POINT'. 3 = 7.08581  ;VaEXPO. 5 .'POIDS'. 3 = .003612 ;
  193. VaEXPO. 5 .'POINT'. 4 = 12.640801 ;
  194. VaEXPO. 5 .'POIDS'. 4 = 2.336997e-05 ;
  195. VaEXPO. 5 .'POINT'. 5 = .26356  ;VaEXPO. 5 .'POIDS'. 5 = .521756 ;
  196. VaEXPO. 6 .'POINT'. 1 = 1.188932 ;VaEXPO. 6 .'POIDS'. 1 = .417001 ;
  197. VaEXPO. 6 .'POINT'. 2 = .222847 ;VaEXPO. 6 .'POIDS'. 2 = .458965 ;
  198. VaEXPO. 6 .'POINT'. 3 = 5.775144 ;VaEXPO. 6 .'POIDS'. 3 = .010399 ;
  199. VaEXPO. 6 .'POINT'. 4 = 15.982874 ;
  200. VaEXPO. 6 .'POIDS'. 4 = 8.985479e-07 ;
  201. VaEXPO. 6 .'POINT'. 5 = 9.837467 ;VaEXPO. 6 .'POIDS'. 5 = .000261 ;
  202. VaEXPO. 6 .'POINT'. 6 = 2.992736 ;VaEXPO. 6 .'POIDS'. 6 = .113373 ;
  203. VaEXPO. 7 .'POINT'. 1 = 1.026665 ;VaEXPO. 7 .'POIDS'. 1 = .421831 ;
  204. VaEXPO. 7 .'POINT'. 2 = .193044 ;VaEXPO. 7 .'POIDS'. 2 = .409319 ;
  205. VaEXPO. 7 .'POINT'. 3 = 12.73418  ;
  206. VaEXPO. 7 .'POIDS'. 3 = 1.586546e-05 ;
  207. VaEXPO. 7 .'POINT'. 4 = 19.395728 ;
  208. VaEXPO. 7 .'POIDS'. 4 = 3.170315e-08 ;
  209. VaEXPO. 7 .'POINT'. 5 = 2.567877 ;VaEXPO. 7 .'POIDS'. 5 = .147126 ;
  210. VaEXPO. 7 .'POINT'. 6 = 4.900353 ;VaEXPO. 7 .'POIDS'. 6 = .020634 ;
  211. VaEXPO. 7 .'POINT'. 7 = 8.182153 ;VaEXPO. 7 .'POIDS'. 7 = .001074 ;
  212. VaEXPO. 8 .'POINT'. 1 = 7.045905 ;VaEXPO. 8 .'POIDS'. 1 = .002795 ;
  213. VaEXPO. 8 .'POINT'. 2 = 2.251087 ;VaEXPO. 8 .'POIDS'. 2 = .175795 ;
  214. VaEXPO. 8 .'POINT'. 3 = .17028  ;VaEXPO. 8 .'POIDS'. 3 = .369189 ;
  215. VaEXPO. 8 .'POINT'. 4 = 22.863132 ;
  216. VaEXPO. 8 .'POIDS'. 4 = 1.048001e-09 ;
  217. VaEXPO. 8 .'POINT'. 5 = .903702 ;VaEXPO. 8 .'POIDS'. 5 = .418787 ;
  218. VaEXPO. 8 .'POINT'. 6 = 4.2667   ;VaEXPO. 8 .'POIDS'. 6 = .033343 ;
  219. VaEXPO. 8 .'POINT'. 7 = 10.758516 ;
  220. VaEXPO. 8 .'POIDS'. 7 = 9.076509e-05 ;
  221. VaEXPO. 8 .'POINT'. 8 = 15.740679 ;
  222. VaEXPO. 8 .'POIDS'. 8 = 8.485747e-07 ;
  223. VaEXPO. 9 .'POINT'. 1 = .80722  ;VaEXPO. 9 .'POIDS'. 1 = .411214 ;
  224. VaEXPO. 9 .'POINT'. 2 = .152322 ;VaEXPO. 9 .'POIDS'. 2 = .336126 ;
  225. VaEXPO. 9 .'POINT'. 3 = 26.374072 ;
  226. VaEXPO. 9 .'POIDS'. 3 = 3.290874e-11 ;
  227. VaEXPO. 9 .'POINT'. 4 = 13.466237 ;
  228. VaEXPO. 9 .'POIDS'. 4 = 6.592123e-06 ;
  229. VaEXPO. 9 .'POINT'. 5 = 18.833598 ;
  230. VaEXPO. 9 .'POIDS'. 5 = 4.110769e-08 ;
  231. VaEXPO. 9 .'POINT'. 6 = 3.783474 ;VaEXPO. 9 .'POIDS'. 6 = .047461 ;
  232. VaEXPO. 9 .'POINT'. 7 = 9.372985 ;VaEXPO. 9 .'POIDS'. 7 = .000305 ;
  233. VaEXPO. 9 .'POINT'. 8 = 6.204957 ;VaEXPO. 9 .'POIDS'. 8 = .005600 ;
  234. VaEXPO. 9 .'POINT'. 9 = 2.005135 ;VaEXPO. 9 .'POIDS'. 9 = .199288 ;
  235. VaEXPO. 10 .'POINT'. 1 = 29.920697 ;
  236. VaEXPO. 10 .'POIDS'. 1 = 9.911827e-13 ;
  237. VaEXPO. 10 .'POINT'. 2 = 21.996586 ;
  238. VaEXPO. 10 .'POIDS'. 2 = 1.839565e-09 ;
  239. VaEXPO. 10 .'POINT'. 3 = 16.279258 ;
  240. VaEXPO. 10 .'POIDS'. 3 = 4.249314e-07 ;
  241. VaEXPO. 10 .'POINT'. 4 = 11.843786 ;
  242. VaEXPO. 10 .'POIDS'. 4 = 2.825923e-05 ;
  243. VaEXPO. 10 .'POINT'. 5 = 8.330153 ;VaEXPO. 10 .'POIDS'. 5 = .000753 ;
  244. VaEXPO. 10 .'POINT'. 6 = 5.552496 ;VaEXPO. 10 .'POIDS'. 6 = .009502 ;
  245. VaEXPO. 10 .'POINT'. 7 = 3.401434 ;VaEXPO. 10 .'POIDS'. 7 = .062087 ;
  246. VaEXPO. 10 .'POINT'. 8 = 1.808343 ;VaEXPO. 10 .'POIDS'. 8 = .218068 ;
  247. VaEXPO. 10 .'POINT'. 9 = .729455 ;VaEXPO. 10 .'POIDS'. 9 = .40112  ;
  248. VaEXPO. 10 .'POINT'. 10 = .137793 ;VaEXPO. 10 .'POIDS'. 10 = .308441 ;
  249. *
  250. *                ___________________
  251. *                
  252. *                   LOI UNIFORME
  253. *                ___________________
  254. VaUNIF = TABLE;
  255. REPE BOUC 10;
  256.   VaUNIF. &BOUC = TABLE;
  257.   VaUNIF. &BOUC .'POINT' = TABLE;
  258.   VaUNIF. &BOUC .'POIDS' = TABLE;
  259. FIN BOUC;
  260. VaUNIF. 1 .'POINT'. 1 = .5      ;VaUNIF. 1 .'POIDS'. 1 = 1        ;
  261. VaUNIF. 2 .'POINT'. 1 = .788675 ;VaUNIF. 2 .'POIDS'. 1 = .5      ;
  262. VaUNIF. 2 .'POINT'. 2 = .211325 ;VaUNIF. 2 .'POIDS'. 2 = .5      ;
  263. VaUNIF. 3 .'POINT'. 1 = .112702 ;VaUNIF. 3 .'POIDS'. 1 = .277778 ;
  264. VaUNIF. 3 .'POINT'. 2 = .887298 ;VaUNIF. 3 .'POIDS'. 2 = .277778 ;
  265. VaUNIF. 3 .'POINT'. 3 = .5      ;VaUNIF. 3 .'POIDS'. 3 = .444444 ;
  266. VaUNIF. 4 .'POINT'. 1 = .069432 ;VaUNIF. 4 .'POIDS'. 1 = .173927 ;
  267. VaUNIF. 4 .'POINT'. 2 = .330009 ;VaUNIF. 4 .'POIDS'. 2 = .326073 ;
  268. VaUNIF. 4 .'POINT'. 3 = .669991 ;VaUNIF. 4 .'POIDS'. 3 = .326073 ;
  269. VaUNIF. 4 .'POINT'. 4 = .930568 ;VaUNIF. 4 .'POIDS'. 4 = .173927 ;
  270. VaUNIF. 5 .'POINT'. 1 = .04691  ;VaUNIF. 5 .'POIDS'. 1 = .118463 ;
  271. VaUNIF. 5 .'POINT'. 2 = .230765 ;VaUNIF. 5 .'POIDS'. 2 = .239314 ;
  272. VaUNIF. 5 .'POINT'. 3 = .5      ;VaUNIF. 5 .'POIDS'. 3 = .284444 ;
  273. VaUNIF. 5 .'POINT'. 4 = .769235 ;VaUNIF. 5 .'POIDS'. 4 = .239314 ;
  274. VaUNIF. 5 .'POINT'. 5 = .95309  ;VaUNIF. 5 .'POIDS'. 5 = .118463 ;
  275. VaUNIF. 6 .'POINT'. 1 = .033765 ;VaUNIF. 6 .'POIDS'. 1 = .085662 ;
  276. VaUNIF. 6 .'POINT'. 2 = .966235 ;VaUNIF. 6 .'POIDS'. 2 = .085662 ;
  277. VaUNIF. 6 .'POINT'. 3 = .830605 ;VaUNIF. 6 .'POIDS'. 3 = .180381 ;
  278. VaUNIF. 6 .'POINT'. 4 = .61931  ;VaUNIF. 6 .'POIDS'. 4 = .233957 ;
  279. VaUNIF. 6 .'POINT'. 5 = .38069  ;VaUNIF. 6 .'POIDS'. 5 = .233957 ;
  280. VaUNIF. 6 .'POINT'. 6 = .169395 ;VaUNIF. 6 .'POIDS'. 6 = .180381 ;
  281. VaUNIF. 7 .'POINT'. 1 = .025446 ;VaUNIF. 7 .'POIDS'. 1 = .064742 ;
  282. VaUNIF. 7 .'POINT'. 2 = .129234 ;VaUNIF. 7 .'POIDS'. 2 = .139853 ;
  283. VaUNIF. 7 .'POINT'. 3 = .870766 ;VaUNIF. 7 .'POIDS'. 3 = .139853 ;
  284. VaUNIF. 7 .'POINT'. 4 = .974554 ;VaUNIF. 7 .'POIDS'. 4 = .064742 ;
  285. VaUNIF. 7 .'POINT'. 5 = .297077 ;VaUNIF. 7 .'POIDS'. 5 = .190915 ;
  286. VaUNIF. 7 .'POINT'. 6 = .5      ;VaUNIF. 7 .'POIDS'. 6 = .20898  ;
  287. VaUNIF. 7 .'POINT'. 7 = .702923 ;VaUNIF. 7 .'POIDS'. 7 = .190915 ;
  288. VaUNIF. 8 .'POINT'. 1 = .762766 ;VaUNIF. 8 .'POIDS'. 1 = .156853 ;
  289. VaUNIF. 8 .'POINT'. 2 = .898333 ;VaUNIF. 8 .'POIDS'. 2 = .111191 ;
  290. VaUNIF. 8 .'POINT'. 3 = .980145 ;VaUNIF. 8 .'POIDS'. 3 = .050614 ;
  291. VaUNIF. 8 .'POINT'. 4 = .591717 ;VaUNIF. 8 .'POIDS'. 4 = .181342 ;
  292. VaUNIF. 8 .'POINT'. 5 = .408283 ;VaUNIF. 8 .'POIDS'. 5 = .181342 ;
  293. VaUNIF. 8 .'POINT'. 6 = .237234 ;VaUNIF. 8 .'POIDS'. 6 = .156853 ;
  294. VaUNIF. 8 .'POINT'. 7 = .101667 ;VaUNIF. 8 .'POIDS'. 7 = .111191 ;
  295. VaUNIF. 8 .'POINT'. 8 = .019855 ;VaUNIF. 8 .'POIDS'. 8 = .050614 ;
  296. VaUNIF. 9 .'POINT'. 1 = .015920 ;VaUNIF. 9 .'POIDS'. 1 = .040637 ;
  297. VaUNIF. 9 .'POINT'. 2 = .193314 ;VaUNIF. 9 .'POIDS'. 2 = .130305 ;
  298. VaUNIF. 9 .'POINT'. 3 = .081984 ;VaUNIF. 9 .'POIDS'. 3 = .090324 ;
  299. VaUNIF. 9 .'POINT'. 4 = .806686 ;VaUNIF. 9 .'POIDS'. 4 = .130305 ;
  300. VaUNIF. 9 .'POINT'. 5 = .918016 ;VaUNIF. 9 .'POIDS'. 5 = .090324 ;
  301. VaUNIF. 9 .'POINT'. 6 = .98408  ;VaUNIF. 9 .'POIDS'. 6 = .040637 ;
  302. VaUNIF. 9 .'POINT'. 7 = .337873 ;VaUNIF. 9 .'POIDS'. 7 = .156174 ;
  303. VaUNIF. 9 .'POINT'. 8 = .5      ;VaUNIF. 9 .'POIDS'. 8 = .16512  ;
  304. VaUNIF. 9 .'POINT'. 9 = .662127 ;VaUNIF. 9 .'POIDS'. 9 = .156174 ;
  305. VaUNIF. 10 .'POINT'. 1 = .839705 ;VaUNIF. 10 .'POIDS'. 1 = .109543 ;
  306. VaUNIF. 10 .'POINT'. 2 = .574437 ;VaUNIF. 10 .'POIDS'. 2 = .147762 ;
  307. VaUNIF. 10 .'POINT'. 3 = .716698 ;VaUNIF. 10 .'POIDS'. 3 = .134633 ;
  308. VaUNIF. 10 .'POINT'. 4 = .425563 ;VaUNIF. 10 .'POIDS'. 4 = .147762 ;
  309. VaUNIF. 10 .'POINT'. 5 = .067468 ;VaUNIF. 10 .'POIDS'. 5 = .074726 ;
  310. VaUNIF. 10 .'POINT'. 6 = .160295 ;VaUNIF. 10 .'POIDS'. 6 = .109543 ;
  311. VaUNIF. 10 .'POINT'. 7 = .283302 ;VaUNIF. 10 .'POIDS'. 7 = .134633 ;
  312. VaUNIF. 10 .'POINT'. 8 = .932532 ;VaUNIF. 10 .'POIDS'. 8 = .074726 ;
  313. VaUNIF. 10 .'POINT'. 9 = .986953 ;VaUNIF. 10 .'POIDS'. 9 = .033336 ;
  314. VaUNIF. 10 .'POINT'. 10 = .013047 ;VaUNIF. 10 .'POIDS'. 10 = .033336 ;
  315. *
  316. * _________________                                                        
  317. *                                                                               
  318. * DEBUT DES CALCULS                                                        
  319. * _________________                                                        
  320. *
  321. RESTAB = TABLE;
  322. REPE BOUC N;
  323.   RESTAB. &BOUC = TABLE;
  324. FIN BOUC;
  325. *
  326. *MESS '%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%';
  327. *MESS '%%% CALCUL DES ' N ' POINTS ET POIDS D INTEGRATION';
  328. *MESS '%%% POUR UNE LOI ' LOI;
  329. *MESS '%%% DE MOYENNE   ' MU;
  330. *MESS '%%% D ECART_TYPE ' SIGMA;
  331. *MESS '%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%';
  332. *
  333. *
  334. SI (EGA LOI 'NORM');
  335. REPE BOUC N;
  336.   RESTAB .&BOUC. 'POINT' = MU
  337.      + ((VaNORM. N .'POINT'. &BOUC) * SIGMA);
  338.   RESTAB .&BOUC. 'POIDS' = VaNORM. N .'POIDS'. &BOUC;
  339. FIN BOUC;
  340. FINSI; 
  341. *
  342. *
  343. SI (EGA LOI 'LOGN');
  344. C = SIGMA / MU;
  345. REPE BOUC N;
  346.   RESTAB .&BOUC. 'POINT' = MU / ((1+(C**2))**0.5)* 
  347.      (EXP ((VaNORM. N .'POINT'. &BOUC) * ((LOG (1+(C**2)))**0.5)));
  348.   RESTAB .&BOUC. 'POIDS' = VaNORM. N .'POIDS'. &BOUC;
  349. FIN BOUC;
  350. FINSI; 
  351. *
  352. *
  353. SI (EGA LOI 'EXPO');
  354. REPE BOUC N;
  355.   RESTAB .&BOUC. 'POINT' = (VaEXPO. N .'POINT'. &BOUC) * MU;
  356.   RESTAB .&BOUC. 'POIDS' = VaEXPO. N .'POIDS'. &BOUC;
  357. FIN BOUC;
  358. FINSI; 
  359. *
  360. *
  361. *
  362. *
  363. SI (EGA LOI 'UNIF');
  364. A = MU - ((3**0.5)*SIGMA);
  365. B = MU + ((3**0.5)*SIGMA);
  366. REPE BOUC N;
  367.   RESTAB .&BOUC. 'POINT' = A + ((VaUNIF. N .'POINT'. &BOUC) * (B-A));
  368.   RESTAB .&BOUC. 'POIDS' = VaUNIF. N .'POIDS'. &BOUC;
  369. FIN BOUC;
  370. FINSI; 
  371. *
  372. *
  373. 'FINPROC' RESTAB;
  374.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales