Télécharger sissib_cov2.dgibi

Retour à la liste

Numérotation des lignes :

  1. OPTI DIME 3 MODE TRID ELEM QUA4;
  2. *
  3. **********************************************************
  4. * Test sur la procedure SISSIB, suite à la modification
  5. * pour le calcul des efforts maximaux et des covariances sur
  6. * base modale des efforts generalisés pour les elements frame
  7. *
  8. * On considere un portique de 6m x 5m
  9. * soumis à trois mouvements en direction x, y et z (representés
  10. * par les spectres de reponse)
  11. * Les poteaux sont encastrés à la base
  12.  
  13. *
  14. * On compare les covariances calculées avec SISSIB avec les
  15. * valeurs
  16. * analytiques
  17. *
  18. * Develloppé par Alberto FRAU /DEN/DANS/DM2S/SEMT/EMSI
  19. * et Nicolas ILE /DEN/DANS/DM2S/SEMT/EMSI
  20. *
  21. **********************************************************
  22. *
  23. * option graphiques
  24. GRAP1 = 0;
  25. *
  26. * Proprietés geometriques Plancher
  27. H1 = 0.15;
  28. *
  29. * Proprietes poteaux
  30. BPOT1 = 0.30;
  31. INR1 = ((1.0/12.0)*(0.30))*((0.30)**(3.0));
  32. SECT1 = ((0.30)*(0.30));
  33. TORS1 = 10.0;
  34. *
  35. * Definition du chargement sismique x, y et z
  36. LIS_FREQ = PROG 0.1 0.25 2.5 8.0 30.0 40.0 50.0 400.0 ;
  37. * liste des acc?l?rations pour seisme suivant x et y [m.s-2]
  38. LIS_ACCX = PROG 0.0175 0.0175 0.184 0.208 0.1 0.1 0.1 0.1 ;
  39. * liste des acc?l?rations pour s?isme suivant z [m.s-2]
  40. LIS_ACCZ = (LIS_ACCX)*0.6666667 ;
  41. *
  42. LIS_ACCX = (LIS_ACCX * 9.810)*3.0 ;
  43. LIS_ACCZ = (LIS_ACCZ * 9.810)*3.0 ;
  44. *
  45. SEISME_X = EVOL MANU '[Hz]' LIS_FREQ '[m.s-2]' LIS_ACCX ;
  46. SEISME_Y = EVOL MANU '[Hz]' LIS_FREQ '[m.s-2]' LIS_ACCX ;
  47. SEISME_Z = EVOL ROUG MANU '[Hz]' LIS_FREQ '[m.s-2]' LIS_ACCZ ;
  48. *
  49. * -------------------------------------------------------------
  50. *
  51. P1 = 0. 0. 0.;
  52. P2 = 6. 0. 0.;
  53. P3 = 6. 5. 0.;
  54. P4 = 0. 5. 0.;
  55. *
  56. L1 = D 30 P1 P2;
  57. L2 = D 20 P2 P3;
  58. L3 = D 30 P3 P4;
  59. L4 = D 20 P4 P1;
  60. *
  61. SUR1 = DALL L1 L2 L3 L4;
  62. SUR1 = SUR1 PLUS (0. 0. 3.0);
  63. POT1 = (P1 LIGNE TRANS 10 (0. 0. 3.0));
  64. POT2 = (P2 LIGNE TRANS 10 (0. 0. 3.0)) ET
  65. (P3 LIGNE TRANS 10 (0. 0. 3.0)) ET
  66. (P4 LIGNE TRANS 10 (0. 0. 3.0));
  67. ELIM 0.0001 (SUR1 ET POT1 ET POT2);
  68. *
  69. * modele et materiau
  70. MOD1 = MODE SUR1 MECANIQUE ELASTIQUE ISOTROPE COQ4;
  71. MAT1 = MATE MOD1 YOUNG 30000000.E6 NU 0.2 RHO 2500.0 EPAI H1;
  72. MOD2 = MODE POT1 MECANIQUE ELASTIQUE ISOTROPE POUT;
  73. MAT2 = MATE MOD2 YOUNG 30000.E6 NU 0.2 RHO 2500.0
  74. SECT SECT1 INRZ INR1 INRY INR1 TORS TORS1 VECT (0. 1.0 0.);
  75. MOD3 = MODE POT2 MECANIQUE ELASTIQUE ISOTROPE POUT;
  76. MAT3 = MATE MOD3 YOUNG 30000.E6 NU 0.2 RHO 2500.0
  77. SECT SECT1 INRZ INR1 INRY INR1 TORS TORS1 VECT (0. 1.0 0.);
  78. *
  79. * matrices de masse et rigidite
  80. RIG1 = (RIGI MOD1 MAT1) ET (RIGI MOD2 MAT2) ET (RIGI MOD3 MAT3);
  81. MAS1 = (MASS MOD1 MAT1) ET (MASS MOD2 MAT2) ET (MASS MOD3 MAT3);
  82. *
  83. * blocages
  84. BL1 = (BLOQUER 'DEPLA' 'ROTA' (P1 ET P2 ET P3 ET P4));
  85. * (BLOQUER 'UY' L3) ET (BLOQUER 'UY' L1);
  86. *
  87. * definition de la base modale
  88. OPTI EPSI LINEAIRE;
  89. TB1 = VIBR INTERVALLE 0. 1000. BASSE 4
  90. (RIG1 ET BL1) (MAS1) IMPR TBAS;
  91. CH1 = TB1.MODES. 1 .DEFORMEE_MODALE;
  92. CH2 = TB1.MODES. 2 .DEFORMEE_MODALE;
  93. CH3 = TB1.MODES. 3 .DEFORMEE_MODALE;
  94. CH4 = TB1.MODES. 4 .DEFORMEE_MODALE;
  95. DEF1 = DEFO CH1 (SUR1 ET POT1 ET POT2) VERT;
  96. DEF2 = DEFO CH2 (SUR1 ET POT1 ET POT2) VERT;
  97. DEF3 = DEFO CH3 (SUR1 ET POT1 ET POT2) VERT;
  98. DEF4 = DEFO CH4 (SUR1 ET POT1 ET POT2) VERT;
  99. DEF0 = DEFO CH1 (SUR1 ET POT1 ET POT2) 0. ROUG;
  100. *TRAC (DEF0 ET DEF1);
  101. *TRAC (DEF0 ET DEF2);
  102. *TRAC (DEF0 ET DEF3);
  103. *TRAC (DEF0 ET DEF4);
  104. *
  105. * calcul contraintes modales
  106. TB1 = SIGSOL (MOD1 ET MOD2 ET MOD3) (MAT1 ET MAT2 ET MAT3) TB1 ;
  107. SIG_M1 = (TB1.MODES. 1 .CONTRAINTE_MODALE);
  108. SIG_M2 = (TB1.MODES. 2 .CONTRAINTE_MODALE);
  109. SIG_M3 = (TB1.MODES. 3 .CONTRAINTE_MODALE);
  110. SIG_M4 = (TB1.MODES. 4 .CONTRAINTE_MODALE);
  111. SIG_M1 = REDU (CHAN GRAVITE SIG_M1 (REDU (MOD1 ET MOD2 ET MOD3) POT1))
  112. (POT1 ELEM 10);
  113. SIG_M2 = REDU (CHAN GRAVITE SIG_M2 (REDU (MOD1 ET MOD2 ET MOD3) POT1))
  114. (POT1 ELEM 10);
  115. SIG_M3 = REDU (CHAN GRAVITE SIG_M3 (REDU (MOD1 ET MOD2 ET MOD3) POT1))
  116. (POT1 ELEM 10);
  117. SIG_M4 = REDU (CHAN GRAVITE SIG_M4 (REDU (MOD1 ET MOD2 ET MOD3) POT1))
  118. (POT1 ELEM 10);
  119. *
  120. * extraction des efforts globaux pour l'element 1
  121. * mode 1
  122. NN_M1 = EXTR SIG_M1 'EFFX' 1 1 1;
  123. VY_M1 = EXTR SIG_M1 'EFFY' 1 1 1;
  124. VZ_M1 = EXTR SIG_M1 'EFFZ' 1 1 1;
  125. MT_M1 = EXTR SIG_M1 'MOMX' 1 1 1;
  126. MY_M1 = EXTR SIG_M1 'MOMY' 1 1 1;
  127. MZ_M1 = EXTR SIG_M1 'MOMZ' 1 1 1;
  128. *
  129. * extraction des efforts globaux pour l'element 1
  130. * mode 2
  131. NN_M2 = EXTR SIG_M2 'EFFX' 1 1 1;
  132. VY_M2 = EXTR SIG_M2 'EFFY' 1 1 1;
  133. VZ_M2 = EXTR SIG_M2 'EFFZ' 1 1 1;
  134. MT_M2 = EXTR SIG_M2 'MOMX' 1 1 1;
  135. MY_M2 = EXTR SIG_M2 'MOMY' 1 1 1;
  136. MZ_M2 = EXTR SIG_M2 'MOMZ' 1 1 1;
  137. *
  138. * extraction des efforts globaux pour l'element 1
  139. * mode 3
  140. NN_M3 = EXTR SIG_M3 'EFFX' 1 1 1;
  141. VY_M3 = EXTR SIG_M3 'EFFY' 1 1 1;
  142. VZ_M3 = EXTR SIG_M3 'EFFZ' 1 1 1;
  143. MT_M3 = EXTR SIG_M3 'MOMX' 1 1 1;
  144. MY_M3 = EXTR SIG_M3 'MOMY' 1 1 1;
  145. MZ_M3 = EXTR SIG_M3 'MOMZ' 1 1 1;
  146. *
  147. * extraction des efforts globaux pour l'element 1
  148. * mode 4
  149. NN_M4 = EXTR SIG_M4 'EFFX' 1 1 1;
  150. VY_M4 = EXTR SIG_M4 'EFFY' 1 1 1;
  151. VZ_M4 = EXTR SIG_M4 'EFFZ' 1 1 1;
  152. MT_M4 = EXTR SIG_M4 'MOMX' 1 1 1;
  153. MY_M4 = EXTR SIG_M4 'MOMY' 1 1 1;
  154. MZ_M4 = EXTR SIG_M4 'MOMZ' 1 1 1;
  155. *
  156. * Calcul des reactions modales
  157. TB1 = REAC BL1 TB1 ;
  158. *
  159. SI (GRAP1 EGA 1);
  160. NN1 = (DIME TB1.MODES) - 2;
  161. I = 1;
  162. REPE IND1 (NN1);
  163. CH1 = TB1.MODES. I .DEFORMEE_MODALE;
  164. FREQ1 = TB1.MODES. I .FREQUENCE;
  165. TIT1 = CHAINE 'Mode ' I ' - Frequence ' FREQ1;
  166. DEF0 = DEFO CH1 (SUR1 ET POT1 ET POT2) 0. ROUG;
  167. DEF1 = DEFO CH1 (SUR1 ET POT1 ET POT2) VERT;
  168. TRAC (1.E6 -1.E6 1.E6) (DEF1 ET DEF0) TITR TIT1;
  169. I = I + 1;
  170. FIN IND1;
  171. FINSI;
  172.  
  173.  
  174.  
  175.  
  176. *
  177. * Reponse base modale
  178. *
  179. NB_MODE = (DIME TB1.MODES) - 2;
  180. * liste des amortissements
  181. LIS_AMOR = PROG 7.e-2 ;
  182. BAS_AMOR = PROG NB_MODE * 7.e-2 ;
  183. *
  184. TAB1 = 'TABLE' ;
  185. TAB3 = 'TABLE' 'EXCITATION' ;
  186. TAB1 . 'STRUCTURE' = TB1 ;
  187. *
  188. TAB1 . 'AMORTISSEMENT' = BAS_AMOR ;
  189. TAB1 . 'EXCITATION' = TAB3 ;
  190. *TAB1 . 'TRONCATURE' = VRAI;
  191. * Excitation suivant x
  192. TAB3 . 1 = 'TABLE' ;
  193. TAB3 . 1 . 'DIRECTION' = 'X' ;
  194. TAB3 . 1 . 'SPECTRE' = SEISME_X ;
  195. TAB3 . 1 . 'AMORTISSEMENT' = LIS_AMOR ;
  196. *
  197. * Excitation suivant y
  198. TAB3 . 2 = 'TABLE' ;
  199. TAB3 . 2 . 'DIRECTION' = 'Y' ;
  200. TAB3 . 2 . 'SPECTRE' = (SEISME_Y) ;
  201. TAB3 . 2 . 'AMORTISSEMENT' = LIS_AMOR ;
  202. *
  203. * Excitation suivant z
  204. TAB3 . 3 = 'TABLE' ;
  205. TAB3 . 3 . 'DIRECTION' = 'Z' ;
  206. TAB3 . 3 . 'SPECTRE' = (SEISME_Z);
  207. TAB3 . 3 . 'AMORTISSEMENT' = LIS_AMOR ;
  208. *
  209. * recombinaison quadratique complete
  210. TAB1 . 'RECOMBINAISON_MODES' = 'CQC';
  211. * Recombinaison des directions de seisme
  212. TAB1 . 'RECOMBINAISON_DIRECTIONS' = 'QUADRATIQUE';
  213. TAB1 . 'IMPRESSION' = FAUX ;
  214. *
  215. TAB4 = 'TABLE' 'SORTIES' ;
  216. TAB1 . 'SORTIES' = TAB4 ;
  217. *
  218. * domaine sur lequel porte les resultats
  219. TAB4 . 'DOMAINE' = MOD1;
  220. TAB4 . 'DEPLACEMENTS' = VRAI ;
  221. TAB4 . 'REACTIONS' = VRAI ;
  222. TAB4 . 'CONTRAINTES' = VRAI ;
  223. TAB4 . 'ACCELERATIONS' = VRAI;
  224. *
  225. TAB5 = TABLE;
  226. TAB5. 1 = TABLE;
  227. TAB5. 1 .'TYPE' = CHAINE 'POUTR';
  228. TAB5. 1 .'MAIL' = POT1;
  229. TAB5. 1 .'MODE' = MOD1 ET MOD2 ET MOD3;
  230. *
  231. * Calcul SISSIB
  232. TAB2 = SISSIB TAB1 TAB5;
  233. *
  234. * Determination de la Pseudo acceleration pour les Modes 1 et 2
  235. * pour les direction x, y et z
  236. PSE_A1X = TAB2.'REPRISE'.'SPECTRE'. 1 . 1;
  237. PSE_A1Y = TAB2.'REPRISE'.'SPECTRE'. 1 . 2;
  238. PSE_A1Z = TAB2.'REPRISE'.'SPECTRE'. 1 . 3;
  239. PSE_A2X = TAB2.'REPRISE'.'SPECTRE'. 2 . 1;
  240. PSE_A2Y = TAB2.'REPRISE'.'SPECTRE'. 2 . 2;
  241. PSE_A2Z = TAB2.'REPRISE'.'SPECTRE'. 2 . 3;
  242. PSE_A3X = TAB2.'REPRISE'.'SPECTRE'. 3 . 1;
  243. PSE_A3Y = TAB2.'REPRISE'.'SPECTRE'. 3 . 2;
  244. PSE_A3Z = TAB2.'REPRISE'.'SPECTRE'. 3 . 3;
  245. PSE_A4X = TAB2.'REPRISE'.'SPECTRE'. 4 . 1;
  246. PSE_A4Y = TAB2.'REPRISE'.'SPECTRE'. 4 . 2;
  247. PSE_A4Z = TAB2.'REPRISE'.'SPECTRE'. 4 . 3;
  248. *
  249. * Determination de l'amplification pour les Modes 1 et 2
  250. * pour les direction x, y et z
  251. AMP_1X = (TB1.MODES. 1 .FACT_PART_NORM. 1)*(PSE_A1X);
  252. AMP_1Y = (TB1.MODES. 1 .FACT_PART_NORM. 2)*(PSE_A1Y);
  253. AMP_1Z = (TB1.MODES. 1 .FACT_PART_NORM. 3)*(PSE_A1Z);
  254. AMP_2X = (TB1.MODES. 2 .FACT_PART_NORM. 1)*(PSE_A2X);
  255. AMP_2Y = (TB1.MODES. 2 .FACT_PART_NORM. 2)*(PSE_A2Y);
  256. AMP_2Z = (TB1.MODES. 2 .FACT_PART_NORM. 3)*(PSE_A2Z);
  257. AMP_3X = (TB1.MODES. 3 .FACT_PART_NORM. 1)*(PSE_A3X);
  258. AMP_3Y = (TB1.MODES. 3 .FACT_PART_NORM. 2)*(PSE_A3Y);
  259. AMP_3Z = (TB1.MODES. 3 .FACT_PART_NORM. 3)*(PSE_A3Z);
  260. AMP_4X = (TB1.MODES. 4 .FACT_PART_NORM. 1)*(PSE_A4X);
  261. AMP_4Y = (TB1.MODES. 4 .FACT_PART_NORM. 2)*(PSE_A4Y);
  262. AMP_4Z = (TB1.MODES. 4 .FACT_PART_NORM. 3)*(PSE_A4Z);
  263. *
  264. * Efforts projetés amplifiés (mode 1 ; directions x, y, z)
  265. NN_M1X = NN_M1*AMP_1X;
  266. VY_M1X = VY_M1*AMP_1X;
  267. VZ_M1X = VZ_M1*AMP_1X;
  268. MT_M1X = MT_M1*AMP_1X;
  269. MY_M1X = MY_M1*AMP_1X;
  270. MZ_M1X = MZ_M1*AMP_1X;
  271. NN_M1Y = NN_M1*AMP_1Y;
  272. VY_M1Y = VY_M1*AMP_1Y;
  273. VZ_M1Y = VZ_M1*AMP_1Y;
  274. MT_M1Y = MT_M1*AMP_1Y;
  275. MY_M1Y = MY_M1*AMP_1Y;
  276. MZ_M1Y = MZ_M1*AMP_1Y;
  277. NN_M1Z = NN_M1*AMP_1Z;
  278. VY_M1Z = VY_M1*AMP_1Z;
  279. VZ_M1Z = VZ_M1*AMP_1Z;
  280. MT_M1Z = MT_M1*AMP_1Z;
  281. MY_M1Z = MY_M1*AMP_1Z;
  282. MZ_M1Z = MZ_M1*AMP_1Z;
  283. * Efforts projetés amplifiés (mode 2 ; directions x, y, z)
  284. NN_M2X = NN_M2*AMP_2X;
  285. VY_M2X = VY_M2*AMP_2X;
  286. VZ_M2X = VZ_M2*AMP_2X;
  287. MT_M2X = MT_M2*AMP_2X;
  288. MY_M2X = MY_M2*AMP_2X;
  289. MZ_M2X = MZ_M2*AMP_2X;
  290. NN_M2Y = NN_M2*AMP_2Y;
  291. VY_M2Y = VY_M2*AMP_2Y;
  292. VZ_M2Y = VZ_M2*AMP_2Y;
  293. MT_M2Y = MT_M2*AMP_2Y;
  294. MY_M2Y = MY_M2*AMP_2Y;
  295. MZ_M2Y = MZ_M2*AMP_2Y;
  296. NN_M2Z = NN_M2*AMP_2Z;
  297. VY_M2Z = VY_M2*AMP_2Z;
  298. VZ_M2Z = VZ_M2*AMP_2Z;
  299. MT_M2Z = MT_M2*AMP_2Z;
  300. MY_M2Z = MY_M2*AMP_2Z;
  301. MZ_M2Z = MZ_M2*AMP_2Z;
  302. * Efforts projetés amplifiés (mode 3 ; directions x, y, z)
  303. NN_M3X = NN_M3*AMP_3X;
  304. VY_M3X = VY_M3*AMP_3X;
  305. VZ_M3X = VZ_M3*AMP_3X;
  306. MT_M3X = MT_M3*AMP_3X;
  307. MY_M3X = MY_M3*AMP_3X;
  308. MZ_M3X = MZ_M3*AMP_3X;
  309. NN_M3Y = NN_M3*AMP_3Y;
  310. VY_M3Y = VY_M3*AMP_3Y;
  311. VZ_M3Y = VZ_M3*AMP_3Y;
  312. MT_M3Y = MT_M3*AMP_3Y;
  313. MY_M3Y = MY_M3*AMP_3Y;
  314. MZ_M3Y = MZ_M3*AMP_3Y;
  315. NN_M3Z = NN_M3*AMP_3Z;
  316. VY_M3Z = VY_M3*AMP_3Z;
  317. VZ_M3Z = VZ_M3*AMP_3Z;
  318. MT_M3Z = MT_M3*AMP_3Z;
  319. MY_M3Z = MY_M3*AMP_3Z;
  320. MZ_M3Z = MZ_M3*AMP_3Z;
  321. * Efforts projetés amplifiés (mode 3 ; directions x, y, z)
  322. NN_M4X = NN_M4*AMP_4X;
  323. VY_M4X = VY_M4*AMP_4X;
  324. VZ_M4X = VZ_M4*AMP_4X;
  325. MT_M4X = MT_M4*AMP_4X;
  326. MY_M4X = MY_M4*AMP_4X;
  327. MZ_M4X = MZ_M4*AMP_4X;
  328. NN_M4Y = NN_M4*AMP_4Y;
  329. VY_M4Y = VY_M4*AMP_4Y;
  330. VZ_M4Y = VZ_M4*AMP_4Y;
  331. MT_M4Y = MT_M4*AMP_4Y;
  332. MY_M4Y = MY_M4*AMP_4Y;
  333. MZ_M4Y = MZ_M4*AMP_4Y;
  334. NN_M4Z = NN_M4*AMP_4Z;
  335. VY_M4Z = VY_M4*AMP_4Z;
  336. VZ_M4Z = VZ_M4*AMP_4Z;
  337. MT_M4Z = MT_M4*AMP_4Z;
  338. MY_M4Z = MY_M4*AMP_4Z;
  339. MZ_M4Z = MZ_M4*AMP_4Z;
  340. *
  341. * determination du terme NNNN
  342. CNNNN_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1X)*(NN_M1X)));
  343. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2X)*(NN_M1X)));
  344. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3X)*(NN_M1X)));
  345. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4X)*(NN_M1X)));
  346. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1X)*(NN_M2X)));
  347. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2X)*(NN_M2X)));
  348. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3X)*(NN_M2X)));
  349. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4X)*(NN_M2X)));
  350. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1X)*(NN_M3X)));
  351. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2X)*(NN_M3X)));
  352. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3X)*(NN_M3X)));
  353. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4X)*(NN_M3X)));
  354. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1X)*(NN_M4X)));
  355. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2X)*(NN_M4X)));
  356. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3X)*(NN_M4X)));
  357. CNNNN_X = CNNNN_X + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4X)*(NN_M4X)));
  358. CNNNN_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1Y)*(NN_M1Y)));
  359. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2Y)*(NN_M1Y)));
  360. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3Y)*(NN_M1Y)));
  361. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4Y)*(NN_M1Y)));
  362. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1Y)*(NN_M2Y)));
  363. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2Y)*(NN_M2Y)));
  364. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3Y)*(NN_M2Y)));
  365. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4Y)*(NN_M2Y)));
  366. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1Y)*(NN_M3Y)));
  367. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2Y)*(NN_M3Y)));
  368. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3Y)*(NN_M3Y)));
  369. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4Y)*(NN_M3Y)));
  370. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1Y)*(NN_M4Y)));
  371. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2Y)*(NN_M4Y)));
  372. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3Y)*(NN_M4Y)));
  373. CNNNN_Y = CNNNN_Y + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4Y)*(NN_M4Y)));
  374. CNNNN_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1Z)*(NN_M1Z)));
  375. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2Z)*(NN_M1Z)));
  376. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3Z)*(NN_M1Z)));
  377. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4Z)*(NN_M1Z)));
  378. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1Z)*(NN_M2Z)));
  379. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2Z)*(NN_M2Z)));
  380. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3Z)*(NN_M2Z)));
  381. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4Z)*(NN_M2Z)));
  382. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1Z)*(NN_M3Z)));
  383. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2Z)*(NN_M3Z)));
  384. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3Z)*(NN_M3Z)));
  385. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4Z)*(NN_M3Z)));
  386. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1Z)*(NN_M4Z)));
  387. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2Z)*(NN_M4Z)));
  388. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3Z)*(NN_M4Z)));
  389. CNNNN_Z = CNNNN_Z + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4Z)*(NN_M4Z)));
  390. CNNNN_T = CNNNN_X + CNNNN_Y + CNNNN_Z;
  391. CNNNN_C = EXTR
  392. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.POU1) (POT1 ELEM 10))
  393. 'C11' 1 1 1;
  394. *
  395. * determination du terme NNVY
  396. CNNVY_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1X)*(VY_M1X)));
  397. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2X)*(VY_M1X)));
  398. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3X)*(VY_M1X)));
  399. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4X)*(VY_M1X)));
  400. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1X)*(VY_M2X)));
  401. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2X)*(VY_M2X)));
  402. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3X)*(VY_M2X)));
  403. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4X)*(VY_M2X)));
  404. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1X)*(VY_M3X)));
  405. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2X)*(VY_M3X)));
  406. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3X)*(VY_M3X)));
  407. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4X)*(VY_M3X)));
  408. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1X)*(VY_M4X)));
  409. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2X)*(VY_M4X)));
  410. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3X)*(VY_M4X)));
  411. CNNVY_X = CNNVY_X + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4X)*(VY_M4X)));
  412. CNNVY_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1Y)*(VY_M1Y)));
  413. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2Y)*(VY_M1Y)));
  414. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3Y)*(VY_M1Y)));
  415. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4Y)*(VY_M1Y)));
  416. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1Y)*(VY_M2Y)));
  417. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2Y)*(VY_M2Y)));
  418. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3Y)*(VY_M2Y)));
  419. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4Y)*(VY_M2Y)));
  420. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1Y)*(VY_M3Y)));
  421. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2Y)*(VY_M3Y)));
  422. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3Y)*(VY_M3Y)));
  423. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4Y)*(VY_M3Y)));
  424. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1Y)*(VY_M4Y)));
  425. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2Y)*(VY_M4Y)));
  426. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3Y)*(VY_M4Y)));
  427. CNNVY_Y = CNNVY_Y + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4Y)*(VY_M4Y)));
  428. CNNVY_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1Z)*(VY_M1Z)));
  429. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2Z)*(VY_M1Z)));
  430. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3Z)*(VY_M1Z)));
  431. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4Z)*(VY_M1Z)));
  432. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1Z)*(VY_M2Z)));
  433. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2Z)*(VY_M2Z)));
  434. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3Z)*(VY_M2Z)));
  435. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4Z)*(VY_M2Z)));
  436. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1Z)*(VY_M3Z)));
  437. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2Z)*(VY_M3Z)));
  438. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3Z)*(VY_M3Z)));
  439. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4Z)*(VY_M3Z)));
  440. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1Z)*(VY_M4Z)));
  441. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2Z)*(VY_M4Z)));
  442. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3Z)*(VY_M4Z)));
  443. CNNVY_Z = CNNVY_Z + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4Z)*(VY_M4Z)));
  444. CNNVY_T = CNNVY_X + CNNVY_Y + CNNVY_Z;
  445. CNNVY_C = EXTR
  446. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.POU1) (POT1 ELEM 10))
  447. 'C12' 1 1 1;
  448. *
  449. * determination du terme NNVZ
  450. CNNVZ_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1X)*(VZ_M1X)));
  451. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2X)*(VZ_M1X)));
  452. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3X)*(VZ_M1X)));
  453. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4X)*(VZ_M1X)));
  454. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1X)*(VZ_M2X)));
  455. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2X)*(VZ_M2X)));
  456. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3X)*(VZ_M2X)));
  457. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4X)*(VZ_M2X)));
  458. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1X)*(VZ_M3X)));
  459. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2X)*(VZ_M3X)));
  460. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3X)*(VZ_M3X)));
  461. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4X)*(VZ_M3X)));
  462. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1X)*(VZ_M4X)));
  463. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2X)*(VZ_M4X)));
  464. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3X)*(VZ_M4X)));
  465. CNNVZ_X = CNNVZ_X + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4X)*(VZ_M4X)));
  466. CNNVZ_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1Y)*(VZ_M1Y)));
  467. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2Y)*(VZ_M1Y)));
  468. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3Y)*(VZ_M1Y)));
  469. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4Y)*(VZ_M1Y)));
  470. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1Y)*(VZ_M2Y)));
  471. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2Y)*(VZ_M2Y)));
  472. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3Y)*(VZ_M2Y)));
  473. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4Y)*(VZ_M2Y)));
  474. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1Y)*(VZ_M3Y)));
  475. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2Y)*(VZ_M3Y)));
  476. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3Y)*(VZ_M3Y)));
  477. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4Y)*(VZ_M3Y)));
  478. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1Y)*(VZ_M4Y)));
  479. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2Y)*(VZ_M4Y)));
  480. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3Y)*(VZ_M4Y)));
  481. CNNVZ_Y = CNNVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4Y)*(VZ_M4Y)));
  482. CNNVZ_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1Z)*(VZ_M1Z)));
  483. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2Z)*(VZ_M1Z)));
  484. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3Z)*(VZ_M1Z)));
  485. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4Z)*(VZ_M1Z)));
  486. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1Z)*(VZ_M2Z)));
  487. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2Z)*(VZ_M2Z)));
  488. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3Z)*(VZ_M2Z)));
  489. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4Z)*(VZ_M2Z)));
  490. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1Z)*(VZ_M3Z)));
  491. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2Z)*(VZ_M3Z)));
  492. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3Z)*(VZ_M3Z)));
  493. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4Z)*(VZ_M3Z)));
  494. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1Z)*(VZ_M4Z)));
  495. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2Z)*(VZ_M4Z)));
  496. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3Z)*(VZ_M4Z)));
  497. CNNVZ_Z = CNNVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4Z)*(VZ_M4Z)));
  498. CNNVZ_T = CNNVZ_X + CNNVZ_Y + CNNVZ_Z;
  499. CNNVZ_C = EXTR
  500. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.POU1) (POT1 ELEM 10))
  501. 'C13' 1 1 1;
  502. *
  503. * determination du terme NNMT
  504. CNNMT_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1X)*(MT_M1X)));
  505. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2X)*(MT_M1X)));
  506. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3X)*(MT_M1X)));
  507. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4X)*(MT_M1X)));
  508. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1X)*(MT_M2X)));
  509. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2X)*(MT_M2X)));
  510. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3X)*(MT_M2X)));
  511. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4X)*(MT_M2X)));
  512. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1X)*(MT_M3X)));
  513. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2X)*(MT_M3X)));
  514. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3X)*(MT_M3X)));
  515. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4X)*(MT_M3X)));
  516. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1X)*(MT_M4X)));
  517. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2X)*(MT_M4X)));
  518. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3X)*(MT_M4X)));
  519. CNNMT_X = CNNMT_X + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4X)*(MT_M4X)));
  520. CNNMT_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1Y)*(MT_M1Y)));
  521. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2Y)*(MT_M1Y)));
  522. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3Y)*(MT_M1Y)));
  523. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4Y)*(MT_M1Y)));
  524. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1Y)*(MT_M2Y)));
  525. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2Y)*(MT_M2Y)));
  526. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3Y)*(MT_M2Y)));
  527. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4Y)*(MT_M2Y)));
  528. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1Y)*(MT_M3Y)));
  529. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2Y)*(MT_M3Y)));
  530. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3Y)*(MT_M3Y)));
  531. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4Y)*(MT_M3Y)));
  532. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1Y)*(MT_M4Y)));
  533. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2Y)*(MT_M4Y)));
  534. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3Y)*(MT_M4Y)));
  535. CNNMT_Y = CNNMT_Y + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4Y)*(MT_M4Y)));
  536. CNNMT_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1Z)*(MT_M1Z)));
  537. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2Z)*(MT_M1Z)));
  538. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3Z)*(MT_M1Z)));
  539. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4Z)*(MT_M1Z)));
  540. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1Z)*(MT_M2Z)));
  541. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2Z)*(MT_M2Z)));
  542. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3Z)*(MT_M2Z)));
  543. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4Z)*(MT_M2Z)));
  544. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1Z)*(MT_M3Z)));
  545. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2Z)*(MT_M3Z)));
  546. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3Z)*(MT_M3Z)));
  547. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4Z)*(MT_M3Z)));
  548. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1Z)*(MT_M4Z)));
  549. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2Z)*(MT_M4Z)));
  550. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3Z)*(MT_M4Z)));
  551. CNNMT_Z = CNNMT_Z + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4Z)*(MT_M4Z)));
  552. CNNMT_T = CNNMT_X + CNNMT_Y + CNNMT_Z;
  553. CNNMT_C = EXTR
  554. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.POU1) (POT1 ELEM 10))
  555. 'C14' 1 1 1;
  556. *
  557. * determination du terme NNMY
  558. CNNMY_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1X)*(MY_M1X)));
  559. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2X)*(MY_M1X)));
  560. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3X)*(MY_M1X)));
  561. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4X)*(MY_M1X)));
  562. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1X)*(MY_M2X)));
  563. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2X)*(MY_M2X)));
  564. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3X)*(MY_M2X)));
  565. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4X)*(MY_M2X)));
  566. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1X)*(MY_M3X)));
  567. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2X)*(MY_M3X)));
  568. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3X)*(MY_M3X)));
  569. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4X)*(MY_M3X)));
  570. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1X)*(MY_M4X)));
  571. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2X)*(MY_M4X)));
  572. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3X)*(MY_M4X)));
  573. CNNMY_X = CNNMY_X + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4X)*(MY_M4X)));
  574. CNNMY_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1Y)*(MY_M1Y)));
  575. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2Y)*(MY_M1Y)));
  576. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3Y)*(MY_M1Y)));
  577. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4Y)*(MY_M1Y)));
  578. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1Y)*(MY_M2Y)));
  579. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2Y)*(MY_M2Y)));
  580. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3Y)*(MY_M2Y)));
  581. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4Y)*(MY_M2Y)));
  582. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1Y)*(MY_M3Y)));
  583. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2Y)*(MY_M3Y)));
  584. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3Y)*(MY_M3Y)));
  585. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4Y)*(MY_M3Y)));
  586. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1Y)*(MY_M4Y)));
  587. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2Y)*(MY_M4Y)));
  588. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3Y)*(MY_M4Y)));
  589. CNNMY_Y = CNNMY_Y + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4Y)*(MY_M4Y)));
  590. CNNMY_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1Z)*(MY_M1Z)));
  591. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2Z)*(MY_M1Z)));
  592. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3Z)*(MY_M1Z)));
  593. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4Z)*(MY_M1Z)));
  594. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1Z)*(MY_M2Z)));
  595. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2Z)*(MY_M2Z)));
  596. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3Z)*(MY_M2Z)));
  597. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4Z)*(MY_M2Z)));
  598. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1Z)*(MY_M3Z)));
  599. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2Z)*(MY_M3Z)));
  600. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3Z)*(MY_M3Z)));
  601. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4Z)*(MY_M3Z)));
  602. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1Z)*(MY_M4Z)));
  603. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2Z)*(MY_M4Z)));
  604. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3Z)*(MY_M4Z)));
  605. CNNMY_Z = CNNMY_Z + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4Z)*(MY_M4Z)));
  606. CNNMY_T = CNNMY_X + CNNMY_Y + CNNMY_Z;
  607. CNNMY_C = EXTR
  608. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.POU1) (POT1 ELEM 10))
  609. 'C15' 1 1 1;
  610. *
  611. * determination du terme NNMZ
  612. CNNMZ_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1X)*(MZ_M1X)));
  613. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2X)*(MZ_M1X)));
  614. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3X)*(MZ_M1X)));
  615. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4X)*(MZ_M1X)));
  616. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1X)*(MZ_M2X)));
  617. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2X)*(MZ_M2X)));
  618. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3X)*(MZ_M2X)));
  619. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4X)*(MZ_M2X)));
  620. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1X)*(MZ_M3X)));
  621. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2X)*(MZ_M3X)));
  622. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3X)*(MZ_M3X)));
  623. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4X)*(MZ_M3X)));
  624. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1X)*(MZ_M4X)));
  625. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2X)*(MZ_M4X)));
  626. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3X)*(MZ_M4X)));
  627. CNNMZ_X = CNNMZ_X + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4X)*(MZ_M4X)));
  628. CNNMZ_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1Y)*(MZ_M1Y)));
  629. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2Y)*(MZ_M1Y)));
  630. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3Y)*(MZ_M1Y)));
  631. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4Y)*(MZ_M1Y)));
  632. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1Y)*(MZ_M2Y)));
  633. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2Y)*(MZ_M2Y)));
  634. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3Y)*(MZ_M2Y)));
  635. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4Y)*(MZ_M2Y)));
  636. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1Y)*(MZ_M3Y)));
  637. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2Y)*(MZ_M3Y)));
  638. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3Y)*(MZ_M3Y)));
  639. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4Y)*(MZ_M3Y)));
  640. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1Y)*(MZ_M4Y)));
  641. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2Y)*(MZ_M4Y)));
  642. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3Y)*(MZ_M4Y)));
  643. CNNMZ_Y = CNNMZ_Y + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4Y)*(MZ_M4Y)));
  644. CNNMZ_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((NN_M1Z)*(MZ_M1Z)));
  645. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M2Z)*(MZ_M1Z)));
  646. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M3Z)*(MZ_M1Z)));
  647. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M4Z)*(MZ_M1Z)));
  648. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((NN_M1Z)*(MZ_M2Z)));
  649. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 2 . 2)*((NN_M2Z)*(MZ_M2Z)));
  650. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M3Z)*(MZ_M2Z)));
  651. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M4Z)*(MZ_M2Z)));
  652. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((NN_M1Z)*(MZ_M3Z)));
  653. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((NN_M2Z)*(MZ_M3Z)));
  654. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 3 . 3)*((NN_M3Z)*(MZ_M3Z)));
  655. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M4Z)*(MZ_M3Z)));
  656. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((NN_M1Z)*(MZ_M4Z)));
  657. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((NN_M2Z)*(MZ_M4Z)));
  658. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((NN_M3Z)*(MZ_M4Z)));
  659. CNNMZ_Z = CNNMZ_Z + ((TAB2.REPRISE.PIJ. 4 . 4)*((NN_M4Z)*(MZ_M4Z)));
  660. CNNMZ_T = CNNMZ_X + CNNMZ_Y + CNNMZ_Z;
  661. CNNMZ_C = EXTR
  662. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.POU1) (POT1 ELEM 10))
  663. 'C16' 1 1 1;
  664. *
  665. * determination du terme VYVY
  666. CVYVY_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1X)*(VY_M1X)));
  667. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2X)*(VY_M1X)));
  668. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3X)*(VY_M1X)));
  669. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4X)*(VY_M1X)));
  670. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1X)*(VY_M2X)));
  671. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2X)*(VY_M2X)));
  672. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3X)*(VY_M2X)));
  673. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4X)*(VY_M2X)));
  674. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1X)*(VY_M3X)));
  675. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2X)*(VY_M3X)));
  676. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3X)*(VY_M3X)));
  677. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4X)*(VY_M3X)));
  678. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1X)*(VY_M4X)));
  679. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2X)*(VY_M4X)));
  680. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3X)*(VY_M4X)));
  681. CVYVY_X = CVYVY_X + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4X)*(VY_M4X)));
  682. CVYVY_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1Y)*(VY_M1Y)));
  683. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2Y)*(VY_M1Y)));
  684. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3Y)*(VY_M1Y)));
  685. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4Y)*(VY_M1Y)));
  686. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1Y)*(VY_M2Y)));
  687. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2Y)*(VY_M2Y)));
  688. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3Y)*(VY_M2Y)));
  689. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4Y)*(VY_M2Y)));
  690. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1Y)*(VY_M3Y)));
  691. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2Y)*(VY_M3Y)));
  692. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3Y)*(VY_M3Y)));
  693. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4Y)*(VY_M3Y)));
  694. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1Y)*(VY_M4Y)));
  695. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2Y)*(VY_M4Y)));
  696. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3Y)*(VY_M4Y)));
  697. CVYVY_Y = CVYVY_Y + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4Y)*(VY_M4Y)));
  698. CVYVY_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1Z)*(VY_M1Z)));
  699. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2Z)*(VY_M1Z)));
  700. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3Z)*(VY_M1Z)));
  701. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4Z)*(VY_M1Z)));
  702. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1Z)*(VY_M2Z)));
  703. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2Z)*(VY_M2Z)));
  704. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3Z)*(VY_M2Z)));
  705. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4Z)*(VY_M2Z)));
  706. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1Z)*(VY_M3Z)));
  707. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2Z)*(VY_M3Z)));
  708. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3Z)*(VY_M3Z)));
  709. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4Z)*(VY_M3Z)));
  710. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1Z)*(VY_M4Z)));
  711. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2Z)*(VY_M4Z)));
  712. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3Z)*(VY_M4Z)));
  713. CVYVY_Z = CVYVY_Z + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4Z)*(VY_M4Z)));
  714. CVYVY_T = CVYVY_X + CVYVY_Y + CVYVY_Z;
  715. CVYVY_C = EXTR
  716. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.POU1) (POT1 ELEM 10))
  717. 'C22' 1 1 1;
  718. *
  719. * determination du terme VYVZ
  720. CVYVZ_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1X)*(VZ_M1X)));
  721. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2X)*(VZ_M1X)));
  722. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3X)*(VZ_M1X)));
  723. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4X)*(VZ_M1X)));
  724. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1X)*(VZ_M2X)));
  725. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2X)*(VZ_M2X)));
  726. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3X)*(VZ_M2X)));
  727. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4X)*(VZ_M2X)));
  728. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1X)*(VZ_M3X)));
  729. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2X)*(VZ_M3X)));
  730. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3X)*(VZ_M3X)));
  731. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4X)*(VZ_M3X)));
  732. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1X)*(VZ_M4X)));
  733. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2X)*(VZ_M4X)));
  734. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3X)*(VZ_M4X)));
  735. CVYVZ_X = CVYVZ_X + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4X)*(VZ_M4X)));
  736. CVYVZ_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1Y)*(VZ_M1Y)));
  737. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2Y)*(VZ_M1Y)));
  738. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3Y)*(VZ_M1Y)));
  739. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4Y)*(VZ_M1Y)));
  740. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1Y)*(VZ_M2Y)));
  741. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2Y)*(VZ_M2Y)));
  742. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3Y)*(VZ_M2Y)));
  743. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4Y)*(VZ_M2Y)));
  744. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1Y)*(VZ_M3Y)));
  745. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2Y)*(VZ_M3Y)));
  746. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3Y)*(VZ_M3Y)));
  747. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4Y)*(VZ_M3Y)));
  748. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1Y)*(VZ_M4Y)));
  749. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2Y)*(VZ_M4Y)));
  750. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3Y)*(VZ_M4Y)));
  751. CVYVZ_Y = CVYVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4Y)*(VZ_M4Y)));
  752. CVYVZ_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1Z)*(VZ_M1Z)));
  753. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2Z)*(VZ_M1Z)));
  754. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3Z)*(VZ_M1Z)));
  755. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4Z)*(VZ_M1Z)));
  756. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1Z)*(VZ_M2Z)));
  757. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2Z)*(VZ_M2Z)));
  758. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3Z)*(VZ_M2Z)));
  759. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4Z)*(VZ_M2Z)));
  760. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1Z)*(VZ_M3Z)));
  761. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2Z)*(VZ_M3Z)));
  762. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3Z)*(VZ_M3Z)));
  763. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4Z)*(VZ_M3Z)));
  764. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1Z)*(VZ_M4Z)));
  765. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2Z)*(VZ_M4Z)));
  766. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3Z)*(VZ_M4Z)));
  767. CVYVZ_Z = CVYVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4Z)*(VZ_M4Z)));
  768. CVYVZ_T = CVYVZ_X + CVYVZ_Y + CVYVZ_Z;
  769. CVYVZ_C = EXTR
  770. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.POU1) (POT1 ELEM 10))
  771. 'C23' 1 1 1;
  772. *
  773. * determination du terme VYMT
  774. CVYMT_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1X)*(MT_M1X)));
  775. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2X)*(MT_M1X)));
  776. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3X)*(MT_M1X)));
  777. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4X)*(MT_M1X)));
  778. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1X)*(MT_M2X)));
  779. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2X)*(MT_M2X)));
  780. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3X)*(MT_M2X)));
  781. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4X)*(MT_M2X)));
  782. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1X)*(MT_M3X)));
  783. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2X)*(MT_M3X)));
  784. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3X)*(MT_M3X)));
  785. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4X)*(MT_M3X)));
  786. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1X)*(MT_M4X)));
  787. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2X)*(MT_M4X)));
  788. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3X)*(MT_M4X)));
  789. CVYMT_X = CVYMT_X + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4X)*(MT_M4X)));
  790. CVYMT_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1Y)*(MT_M1Y)));
  791. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2Y)*(MT_M1Y)));
  792. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3Y)*(MT_M1Y)));
  793. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4Y)*(MT_M1Y)));
  794. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1Y)*(MT_M2Y)));
  795. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2Y)*(MT_M2Y)));
  796. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3Y)*(MT_M2Y)));
  797. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4Y)*(MT_M2Y)));
  798. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1Y)*(MT_M3Y)));
  799. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2Y)*(MT_M3Y)));
  800. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3Y)*(MT_M3Y)));
  801. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4Y)*(MT_M3Y)));
  802. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1Y)*(MT_M4Y)));
  803. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2Y)*(MT_M4Y)));
  804. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3Y)*(MT_M4Y)));
  805. CVYMT_Y = CVYMT_Y + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4Y)*(MT_M4Y)));
  806. CVYMT_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1Z)*(MT_M1Z)));
  807. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2Z)*(MT_M1Z)));
  808. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3Z)*(MT_M1Z)));
  809. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4Z)*(MT_M1Z)));
  810. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1Z)*(MT_M2Z)));
  811. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2Z)*(MT_M2Z)));
  812. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3Z)*(MT_M2Z)));
  813. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4Z)*(MT_M2Z)));
  814. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1Z)*(MT_M3Z)));
  815. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2Z)*(MT_M3Z)));
  816. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3Z)*(MT_M3Z)));
  817. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4Z)*(MT_M3Z)));
  818. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1Z)*(MT_M4Z)));
  819. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2Z)*(MT_M4Z)));
  820. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3Z)*(MT_M4Z)));
  821. CVYMT_Z = CVYMT_Z + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4Z)*(MT_M4Z)));
  822. CVYMT_T = CVYMT_X + CVYMT_Y + CVYMT_Z;
  823. CVYMT_C = EXTR
  824. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.POU1) (POT1 ELEM 10))
  825. 'C24' 1 1 1;
  826. *
  827. * determination du terme VYMY
  828. CVYMY_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1X)*(MY_M1X)));
  829. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2X)*(MY_M1X)));
  830. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3X)*(MY_M1X)));
  831. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4X)*(MY_M1X)));
  832. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1X)*(MY_M2X)));
  833. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2X)*(MY_M2X)));
  834. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3X)*(MY_M2X)));
  835. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4X)*(MY_M2X)));
  836. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1X)*(MY_M3X)));
  837. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2X)*(MY_M3X)));
  838. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3X)*(MY_M3X)));
  839. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4X)*(MY_M3X)));
  840. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1X)*(MY_M4X)));
  841. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2X)*(MY_M4X)));
  842. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3X)*(MY_M4X)));
  843. CVYMY_X = CVYMY_X + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4X)*(MY_M4X)));
  844. CVYMY_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1Y)*(MY_M1Y)));
  845. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2Y)*(MY_M1Y)));
  846. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3Y)*(MY_M1Y)));
  847. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4Y)*(MY_M1Y)));
  848. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1Y)*(MY_M2Y)));
  849. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2Y)*(MY_M2Y)));
  850. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3Y)*(MY_M2Y)));
  851. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4Y)*(MY_M2Y)));
  852. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1Y)*(MY_M3Y)));
  853. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2Y)*(MY_M3Y)));
  854. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3Y)*(MY_M3Y)));
  855. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4Y)*(MY_M3Y)));
  856. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1Y)*(MY_M4Y)));
  857. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2Y)*(MY_M4Y)));
  858. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3Y)*(MY_M4Y)));
  859. CVYMY_Y = CVYMY_Y + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4Y)*(MY_M4Y)));
  860. CVYMY_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1Z)*(MY_M1Z)));
  861. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2Z)*(MY_M1Z)));
  862. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3Z)*(MY_M1Z)));
  863. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4Z)*(MY_M1Z)));
  864. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1Z)*(MY_M2Z)));
  865. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2Z)*(MY_M2Z)));
  866. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3Z)*(MY_M2Z)));
  867. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4Z)*(MY_M2Z)));
  868. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1Z)*(MY_M3Z)));
  869. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2Z)*(MY_M3Z)));
  870. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3Z)*(MY_M3Z)));
  871. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4Z)*(MY_M3Z)));
  872. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1Z)*(MY_M4Z)));
  873. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2Z)*(MY_M4Z)));
  874. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3Z)*(MY_M4Z)));
  875. CVYMY_Z = CVYMY_Z + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4Z)*(MY_M4Z)));
  876. CVYMY_T = CVYMY_X + CVYMY_Y + CVYMY_Z;
  877. CVYMY_C = EXTR
  878. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.POU1) (POT1 ELEM 10))
  879. 'C25' 1 1 1;
  880. *
  881. * determination du terme VYMZ
  882. CVYMZ_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1X)*(MZ_M1X)));
  883. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2X)*(MZ_M1X)));
  884. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3X)*(MZ_M1X)));
  885. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4X)*(MZ_M1X)));
  886. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1X)*(MZ_M2X)));
  887. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2X)*(MZ_M2X)));
  888. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3X)*(MZ_M2X)));
  889. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4X)*(MZ_M2X)));
  890. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1X)*(MZ_M3X)));
  891. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2X)*(MZ_M3X)));
  892. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3X)*(MZ_M3X)));
  893. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4X)*(MZ_M3X)));
  894. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1X)*(MZ_M4X)));
  895. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2X)*(MZ_M4X)));
  896. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3X)*(MZ_M4X)));
  897. CVYMZ_X = CVYMZ_X + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4X)*(MZ_M4X)));
  898. CVYMZ_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1Y)*(MZ_M1Y)));
  899. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2Y)*(MZ_M1Y)));
  900. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3Y)*(MZ_M1Y)));
  901. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4Y)*(MZ_M1Y)));
  902. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1Y)*(MZ_M2Y)));
  903. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2Y)*(MZ_M2Y)));
  904. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3Y)*(MZ_M2Y)));
  905. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4Y)*(MZ_M2Y)));
  906. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1Y)*(MZ_M3Y)));
  907. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2Y)*(MZ_M3Y)));
  908. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3Y)*(MZ_M3Y)));
  909. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4Y)*(MZ_M3Y)));
  910. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1Y)*(MZ_M4Y)));
  911. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2Y)*(MZ_M4Y)));
  912. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3Y)*(MZ_M4Y)));
  913. CVYMZ_Y = CVYMZ_Y + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4Y)*(MZ_M4Y)));
  914. CVYMZ_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((VY_M1Z)*(MZ_M1Z)));
  915. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M2Z)*(MZ_M1Z)));
  916. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M3Z)*(MZ_M1Z)));
  917. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M4Z)*(MZ_M1Z)));
  918. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VY_M1Z)*(MZ_M2Z)));
  919. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 2 . 2)*((VY_M2Z)*(MZ_M2Z)));
  920. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M3Z)*(MZ_M2Z)));
  921. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M4Z)*(MZ_M2Z)));
  922. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VY_M1Z)*(MZ_M3Z)));
  923. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VY_M2Z)*(MZ_M3Z)));
  924. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 3 . 3)*((VY_M3Z)*(MZ_M3Z)));
  925. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M4Z)*(MZ_M3Z)));
  926. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VY_M1Z)*(MZ_M4Z)));
  927. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VY_M2Z)*(MZ_M4Z)));
  928. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VY_M3Z)*(MZ_M4Z)));
  929. CVYMZ_Z = CVYMZ_Z + ((TAB2.REPRISE.PIJ. 4 . 4)*((VY_M4Z)*(MZ_M4Z)));
  930. CVYMZ_T = CVYMZ_X + CVYMZ_Y + CVYMZ_Z;
  931. CVYMZ_C = EXTR
  932. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.POU1) (POT1 ELEM 10))
  933. 'C26' 1 1 1;
  934. *
  935. * determination du terme VZVZ
  936. CVZVZ_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((VZ_M1X)*(VZ_M1X)));
  937. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M2X)*(VZ_M1X)));
  938. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M3X)*(VZ_M1X)));
  939. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M4X)*(VZ_M1X)));
  940. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M1X)*(VZ_M2X)));
  941. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((VZ_M2X)*(VZ_M2X)));
  942. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M3X)*(VZ_M2X)));
  943. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M4X)*(VZ_M2X)));
  944. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M1X)*(VZ_M3X)));
  945. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M2X)*(VZ_M3X)));
  946. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((VZ_M3X)*(VZ_M3X)));
  947. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M4X)*(VZ_M3X)));
  948. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M1X)*(VZ_M4X)));
  949. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M2X)*(VZ_M4X)));
  950. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M3X)*(VZ_M4X)));
  951. CVZVZ_X = CVZVZ_X + ((TAB2.REPRISE.PIJ. 4 . 4)*((VZ_M4X)*(VZ_M4X)));
  952. CVZVZ_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((VZ_M1Y)*(VZ_M1Y)));
  953. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M2Y)*(VZ_M1Y)));
  954. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M3Y)*(VZ_M1Y)));
  955. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M4Y)*(VZ_M1Y)));
  956. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M1Y)*(VZ_M2Y)));
  957. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 2 . 2)*((VZ_M2Y)*(VZ_M2Y)));
  958. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M3Y)*(VZ_M2Y)));
  959. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M4Y)*(VZ_M2Y)));
  960. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M1Y)*(VZ_M3Y)));
  961. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M2Y)*(VZ_M3Y)));
  962. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 3 . 3)*((VZ_M3Y)*(VZ_M3Y)));
  963. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M4Y)*(VZ_M3Y)));
  964. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M1Y)*(VZ_M4Y)));
  965. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M2Y)*(VZ_M4Y)));
  966. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M3Y)*(VZ_M4Y)));
  967. CVZVZ_Y = CVZVZ_Y + ((TAB2.REPRISE.PIJ. 4 . 4)*((VZ_M4Y)*(VZ_M4Y)));
  968. CVZVZ_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((VZ_M1Z)*(VZ_M1Z)));
  969. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M2Z)*(VZ_M1Z)));
  970. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M3Z)*(VZ_M1Z)));
  971. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M4Z)*(VZ_M1Z)));
  972. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M1Z)*(VZ_M2Z)));
  973. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 2 . 2)*((VZ_M2Z)*(VZ_M2Z)));
  974. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M3Z)*(VZ_M2Z)));
  975. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M4Z)*(VZ_M2Z)));
  976. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M1Z)*(VZ_M3Z)));
  977. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M2Z)*(VZ_M3Z)));
  978. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 3 . 3)*((VZ_M3Z)*(VZ_M3Z)));
  979. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M4Z)*(VZ_M3Z)));
  980. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M1Z)*(VZ_M4Z)));
  981. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M2Z)*(VZ_M4Z)));
  982. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M3Z)*(VZ_M4Z)));
  983. CVZVZ_Z = CVZVZ_Z + ((TAB2.REPRISE.PIJ. 4 . 4)*((VZ_M4Z)*(VZ_M4Z)));
  984. CVZVZ_T = CVZVZ_X + CVZVZ_Y + CVZVZ_Z;
  985. CVZVZ_C = EXTR
  986. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.POU1) (POT1 ELEM 10))
  987. 'C33' 1 1 1;
  988. *
  989. * determination du terme VZMT
  990. CVZMT_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((VZ_M1X)*(MT_M1X)));
  991. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M2X)*(MT_M1X)));
  992. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M3X)*(MT_M1X)));
  993. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M4X)*(MT_M1X)));
  994. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M1X)*(MT_M2X)));
  995. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((VZ_M2X)*(MT_M2X)));
  996. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M3X)*(MT_M2X)));
  997. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M4X)*(MT_M2X)));
  998. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M1X)*(MT_M3X)));
  999. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M2X)*(MT_M3X)));
  1000. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((VZ_M3X)*(MT_M3X)));
  1001. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M4X)*(MT_M3X)));
  1002. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M1X)*(MT_M4X)));
  1003. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M2X)*(MT_M4X)));
  1004. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M3X)*(MT_M4X)));
  1005. CVZMT_X = CVZMT_X + ((TAB2.REPRISE.PIJ. 4 . 4)*((VZ_M4X)*(MT_M4X)));
  1006. CVZMT_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((VZ_M1Y)*(MT_M1Y)));
  1007. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M2Y)*(MT_M1Y)));
  1008. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M3Y)*(MT_M1Y)));
  1009. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M4Y)*(MT_M1Y)));
  1010. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M1Y)*(MT_M2Y)));
  1011. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 2 . 2)*((VZ_M2Y)*(MT_M2Y)));
  1012. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M3Y)*(MT_M2Y)));
  1013. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M4Y)*(MT_M2Y)));
  1014. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M1Y)*(MT_M3Y)));
  1015. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M2Y)*(MT_M3Y)));
  1016. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 3 . 3)*((VZ_M3Y)*(MT_M3Y)));
  1017. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M4Y)*(MT_M3Y)));
  1018. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M1Y)*(MT_M4Y)));
  1019. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M2Y)*(MT_M4Y)));
  1020. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M3Y)*(MT_M4Y)));
  1021. CVZMT_Y = CVZMT_Y + ((TAB2.REPRISE.PIJ. 4 . 4)*((VZ_M4Y)*(MT_M4Y)));
  1022. CVZMT_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((VZ_M1Z)*(MT_M1Z)));
  1023. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M2Z)*(MT_M1Z)));
  1024. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M3Z)*(MT_M1Z)));
  1025. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M4Z)*(MT_M1Z)));
  1026. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M1Z)*(MT_M2Z)));
  1027. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 2 . 2)*((VZ_M2Z)*(MT_M2Z)));
  1028. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M3Z)*(MT_M2Z)));
  1029. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M4Z)*(MT_M2Z)));
  1030. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M1Z)*(MT_M3Z)));
  1031. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M2Z)*(MT_M3Z)));
  1032. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 3 . 3)*((VZ_M3Z)*(MT_M3Z)));
  1033. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M4Z)*(MT_M3Z)));
  1034. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M1Z)*(MT_M4Z)));
  1035. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M2Z)*(MT_M4Z)));
  1036. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M3Z)*(MT_M4Z)));
  1037. CVZMT_Z = CVZMT_Z + ((TAB2.REPRISE.PIJ. 4 . 4)*((VZ_M4Z)*(MT_M4Z)));
  1038. CVZMT_T = CVZMT_X + CVZMT_Y + CVZMT_Z;
  1039. CVZMT_C = EXTR
  1040. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.POU1) (POT1 ELEM 10))
  1041. 'C34' 1 1 1;
  1042. *
  1043. * determination du terme VZMY
  1044. CVZMY_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((VZ_M1X)*(MY_M1X)));
  1045. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M2X)*(MY_M1X)));
  1046. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M3X)*(MY_M1X)));
  1047. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M4X)*(MY_M1X)));
  1048. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M1X)*(MY_M2X)));
  1049. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((VZ_M2X)*(MY_M2X)));
  1050. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M3X)*(MY_M2X)));
  1051. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M4X)*(MY_M2X)));
  1052. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M1X)*(MY_M3X)));
  1053. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M2X)*(MY_M3X)));
  1054. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((VZ_M3X)*(MY_M3X)));
  1055. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M4X)*(MY_M3X)));
  1056. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M1X)*(MY_M4X)));
  1057. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M2X)*(MY_M4X)));
  1058. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M3X)*(MY_M4X)));
  1059. CVZMY_X = CVZMY_X + ((TAB2.REPRISE.PIJ. 4 . 4)*((VZ_M4X)*(MY_M4X)));
  1060. CVZMY_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((VZ_M1Y)*(MY_M1Y)));
  1061. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M2Y)*(MY_M1Y)));
  1062. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M3Y)*(MY_M1Y)));
  1063. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M4Y)*(MY_M1Y)));
  1064. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M1Y)*(MY_M2Y)));
  1065. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 2 . 2)*((VZ_M2Y)*(MY_M2Y)));
  1066. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M3Y)*(MY_M2Y)));
  1067. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M4Y)*(MY_M2Y)));
  1068. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M1Y)*(MY_M3Y)));
  1069. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M2Y)*(MY_M3Y)));
  1070. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 3 . 3)*((VZ_M3Y)*(MY_M3Y)));
  1071. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M4Y)*(MY_M3Y)));
  1072. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M1Y)*(MY_M4Y)));
  1073. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M2Y)*(MY_M4Y)));
  1074. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M3Y)*(MY_M4Y)));
  1075. CVZMY_Y = CVZMY_Y + ((TAB2.REPRISE.PIJ. 4 . 4)*((VZ_M4Y)*(MY_M4Y)));
  1076. CVZMY_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((VZ_M1Z)*(MY_M1Z)));
  1077. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M2Z)*(MY_M1Z)));
  1078. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M3Z)*(MY_M1Z)));
  1079. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M4Z)*(MY_M1Z)));
  1080. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M1Z)*(MY_M2Z)));
  1081. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 2 . 2)*((VZ_M2Z)*(MY_M2Z)));
  1082. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M3Z)*(MY_M2Z)));
  1083. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M4Z)*(MY_M2Z)));
  1084. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M1Z)*(MY_M3Z)));
  1085. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M2Z)*(MY_M3Z)));
  1086. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 3 . 3)*((VZ_M3Z)*(MY_M3Z)));
  1087. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M4Z)*(MY_M3Z)));
  1088. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M1Z)*(MY_M4Z)));
  1089. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M2Z)*(MY_M4Z)));
  1090. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 4 . 3)*((VZ_M3Z)*(MY_M4Z)));
  1091. CVZMY_Z = CVZMY_Z + ((TAB2.REPRISE.PIJ. 4 . 4)*((VZ_M4Z)*(MY_M4Z)));
  1092. CVZMY_T = CVZMY_X + CVZMY_Y + CVZMY_Z;
  1093. CVZMY_C = EXTR
  1094. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.POU1) (POT1 ELEM 10))
  1095. 'C35' 1 1 1;
  1096. *
  1097. * determination du terme VZMZ
  1098. CVZMZ_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((VZ_M1X)*(MZ_M1X)));
  1099. CVZMZ_X = CVZMZ_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M2X)*(MZ_M1X)));
  1100. CVZMZ_X = CVZMZ_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M3X)*(MZ_M1X)));
  1101. CVZMZ_X = CVZMZ_X + ((TAB2.REPRISE.PIJ. 4 . 1)*((VZ_M4X)*(MZ_M1X)));
  1102. CVZMZ_X = CVZMZ_X + ((TAB2.REPRISE.PIJ. 2 . 1)*((VZ_M1X)*(MZ_M2X)));
  1103. CVZMZ_X = CVZMZ_X + ((TAB2.REPRISE.PIJ. 2 . 2)*((VZ_M2X)*(MZ_M2X)));
  1104. CVZMZ_X = CVZMZ_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M3X)*(MZ_M2X)));
  1105. CVZMZ_X = CVZMZ_X + ((TAB2.REPRISE.PIJ. 4 . 2)*((VZ_M4X)*(MZ_M2X)));
  1106. CVZMZ_X = CVZMZ_X + ((TAB2.REPRISE.PIJ. 3 . 1)*((VZ_M1X)*(MZ_M3X)));
  1107. CVZMZ_X = CVZMZ_X + ((TAB2.REPRISE.PIJ. 3 . 2)*((VZ_M2X)*(MZ_M3X)));
  1108. CVZMZ_X = CVZMZ_X + ((TAB2.REPRISE.PIJ. 3 . 3)*((VZ_M3X)*(MZ_M3X)));
  1109. CVZMZ_X = CVZMZ_X + ((TAB2.REPRISE.PIJ. 4 . 3&