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.  
  375.  

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