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 - All rights reserved.
Disclaimer