Télécharger kbmmsp_impl2D.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : kbmmsp_impl2D.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ***********************************************************
  5. ***********************************************************
  6. **** APPROCHE VF "Cell-Centred Formulation" pour la ****
  7. **** solution des ****
  8. **** Equations d'Euler pour un gaz parfait multiespes ****
  9. **** Implicit: calcul du jacobien du residu ****
  10. **** ****
  11. **** Cas gaz multiespes, "calorically perfect" ****
  12. **** ****
  13. **** Methodes: AUSM+ ****
  14. **** ****
  15. **** S. Kudriakov DRN/DMT/SEMT/LTMF DEC 2002 ****
  16. ***********************************************************
  17. ***********************************************************
  18.  
  19. 'OPTION' 'DIME' 2 ;
  20. 'OPTION' 'ELEM' QUA4 ;
  21. 'OPTION' 'ECHO' 1 ;
  22. 'OPTION' 'TRAC' 'X' ;
  23.  
  24. *
  25. *** GRAPH
  26. *
  27.  
  28. GRAPH = FAUX ;
  29. GRAPH = VRAI ;
  30.  
  31. ERRTOL = 5.0D-3 ;
  32. DELTA = 0.000001 ;
  33. METO = 'AUSMPLM ' ;
  34. ***************************
  35. ***** DOMAINE SPATIAL ****
  36. ***************************
  37.  
  38.  
  39. A0 = 0.0D0 0.0D0;
  40. A1 = 1.0D0 0.0D0;
  41. A2 = 2.0D0 0.0D0;
  42. A3 = 3.0D0 0.0D0;
  43. ************************
  44. A0A1 = A0 'DROIT' 1 A1;
  45. A1A2 = A1 'DROIT' 1 A2;
  46. A2A3 = A2 'DROIT' 1 A3;
  47. ************************
  48. DOM1 = 'TRANSLATION' A0A1 1 (0.0 1.0) ;
  49. DOM2 = DOM1 'PLUS' (0.0 1.0) ;
  50. DOM3 = DOM2 'PLUS' (0.0 1.0) ;
  51. DOM4 = 'TRANSLATION' A1A2 1 (0.0 1.0) ;
  52. DOM5 = DOM4 'PLUS' (0.0 1.0) ;
  53. DOM6 = DOM5 'PLUS' (0.0 1.0) ;
  54. DOM7 = 'TRANSLATION' A2A3 1 (0.0 1.0) ;
  55. DOM8 = DOM7 'PLUS' (0.0 1.0) ;
  56. DOM9 = DOM8 'PLUS' (0.0 1.0) ;
  57.  
  58.  
  59. ANGLE = 30.0D0 ;
  60. ORIG1 = 0.0D0 0.0D0 ;
  61.  
  62. 'MESSAGE';
  63. 'MESSAGE' (CHAIN 'ANGLE = ' ANGLE);
  64. 'MESSAGE';
  65.  
  66. DOM1 = DOM1 'TOURNER' ANGLE ORIG1 ;
  67. DOM2 = DOM2 'TOURNER' ANGLE ORIG1 ;
  68. DOM3 = DOM3 'TOURNER' ANGLE ORIG1 ;
  69. DOM4 = DOM4 'TOURNER' ANGLE ORIG1 ;
  70. DOM5 = DOM5 'TOURNER' ANGLE ORIG1 ;
  71. DOM6 = DOM6 'TOURNER' ANGLE ORIG1 ;
  72. DOM7 = DOM7 'TOURNER' ANGLE ORIG1 ;
  73. DOM8 = DOM8 'TOURNER' ANGLE ORIG1 ;
  74. DOM9 = DOM9 'TOURNER' ANGLE ORIG1 ;
  75.  
  76.  
  77. DOMTOT = DOM1 'ET' DOM2 'ET' DOM3 'ET' DOM4 'ET' DOM5 'ET'
  78. DOM6 'ET' DOM7 'ET' DOM8 'ET' DOM9 'ELIMINATION' 0.0001 ;
  79.  
  80.  
  81. ***********************************************************
  82. *** These things are necessary in order to work with ***
  83. *** each square separately ***
  84. ***********************************************************
  85.  
  86. $DOMTOT = 'MODELISER' DOMTOT 'EULER';
  87.  
  88. $DOM1 = 'MODELISER' DOM1 'EULER';
  89. $DOM2 = 'MODELISER' DOM2 'EULER';
  90. $DOM3 = 'MODELISER' DOM3 'EULER';
  91. $DOM4 = 'MODELISER' DOM4 'EULER';
  92. $DOM5 = 'MODELISER' DOM5 'EULER';
  93. $DOM6 = 'MODELISER' DOM6 'EULER';
  94. $DOM7 = 'MODELISER' DOM7 'EULER';
  95. $DOM8 = 'MODELISER' DOM8 'EULER';
  96. $DOM9 = 'MODELISER' DOM9 'EULER';
  97.  
  98. TDOMTOT = 'DOMA' $DOMTOT 'VF';
  99.  
  100. TDOM1 = 'DOMA' $DOM1 'VF';
  101. TDOM2 = 'DOMA' $DOM2 'VF';
  102. TDOM3 = 'DOMA' $DOM3 'VF';
  103. TDOM4 = 'DOMA' $DOM4 'VF';
  104. TDOM5 = 'DOMA' $DOM5 'VF';
  105. TDOM6 = 'DOMA' $DOM6 'VF';
  106. TDOM7 = 'DOMA' $DOM7 'VF';
  107. TDOM8 = 'DOMA' $DOM8 'VF';
  108. TDOM9 = 'DOMA' $DOM9 'VF';
  109.  
  110. MDOM1 = TDOM1 . 'QUAF' ;
  111. MDOM2 = TDOM2 . 'QUAF' ;
  112. MDOM3 = TDOM3 . 'QUAF' ;
  113. MDOM4 = TDOM4 . 'QUAF' ;
  114. MDOM5 = TDOM5 . 'QUAF' ;
  115. MDOM6 = TDOM6 . 'QUAF' ;
  116. MDOM7 = TDOM7 . 'QUAF' ;
  117. MDOM8 = TDOM8 . 'QUAF' ;
  118. MDOM9 = TDOM9 . 'QUAF' ;
  119.  
  120. **** old stuff $DOMTOT = 'DOMA' DOMTOT ;
  121.  
  122. MDOMTOT = TDOMTOT . 'QUAF' ;
  123.  
  124. 'ELIMINATION' (MDOMTOT ET MDOM1) 0.0001 ;
  125. 'ELIMINATION' (MDOMTOT ET MDOM2) 0.0001 ;
  126. 'ELIMINATION' (MDOMTOT ET MDOM3) 0.0001 ;
  127. 'ELIMINATION' (MDOMTOT ET MDOM4) 0.0001 ;
  128. 'ELIMINATION' (MDOMTOT ET MDOM5) 0.0001 ;
  129. 'ELIMINATION' (MDOMTOT ET MDOM6) 0.0001 ;
  130. 'ELIMINATION' (MDOMTOT ET MDOM7) 0.0001 ;
  131. 'ELIMINATION' (MDOMTOT ET MDOM8) 0.0001 ;
  132. 'ELIMINATION' (MDOMTOT ET MDOM9) 0.0001 ;
  133.  
  134.  
  135.  
  136.  
  137. ****************************************************
  138. * Creating a mesh WITHOUT mesh
  139. ****************************************************
  140. P1 = 0.5 0.5 ;
  141. ELP1 = 'MANU' 'POI1' P1 ;
  142. MELLIM = 'DIFF' ELP1 ELP1 ;
  143.  
  144. ********************************************************
  145. ******* MODELE DU GAS *****************************
  146. ********************************************************
  147.  
  148. NESP = 4;
  149.  
  150. *
  151. *** GAS: H_2, O_2, H_2O, N_2
  152. *
  153. * CP, CV en J/Kg/K @ T = 3000
  154. *
  155.  
  156. PGAS = 'TABLE' ;
  157.  
  158. PGAS . 'CP' = 'TABLE' ;
  159. PGAS . 'CP' . 'H2 ' = .18729066D+05 ;
  160. PGAS . 'CP' . 'O2 ' = .11886820D+04 ;
  161. PGAS . 'CP' . 'H2O ' = .31209047D+04 ;
  162. PGAS . 'CP' . 'N2 ' = .12993995D+04 ;
  163.  
  164.  
  165.  
  166. PGAS . 'CV' = 'TABLE' ;
  167. PGAS . 'CV' . 'H2 ' = .14571861D+05 ;
  168. PGAS . 'CV' . 'O2 ' = .92885670D+03 ;
  169. PGAS . 'CV' . 'H2O ' = .26589930D+04 ;
  170. PGAS . 'CV' . 'N2 ' = .10024563D+04;
  171.  
  172.  
  173.  
  174. *
  175. **** Especes qui sont dans les equations d'Euler
  176. *
  177.  
  178. PGAS . 'ESPEULE' = 'MOTS' 'H2 ' 'O2 ' 'H2O ' ;
  179.  
  180. *
  181. **** Espece qui n'y est pas
  182. *
  183.  
  184.  
  185. PGAS . 'ESPNEULE' = 'N2 ';
  186.  
  187. ******************************************************
  188. ************ Les CHPOINTs *************************
  189. ******************************************************
  190.  
  191. MCEN = 'DOMA' $DOMTOT CENTRE ;
  192.  
  193. *** density **
  194. RN = 'BRUI' 'BLAN' 'UNIF' 1.11 0.5 MCEN ;
  195. *** pressure **
  196. PN = 'BRUI' 'BLAN' 'UNIF' 1234.1 500.2 MCEN ;
  197. *** CP ****
  198. CP = ('MANUEL' 'CHPO' MCEN 1 'H2'
  199. (PGAS . 'CP' . 'H2 ') 'NATU' 'DISCRET') 'ET'
  200. ('MANUEL' 'CHPO' MCEN 1 'O2'
  201. (PGAS . 'CP' . 'O2 ') 'NATU' 'DISCRET') 'ET'
  202. ('MANUEL' 'CHPO' MCEN 1 'H2O'
  203. (PGAS . 'CP' . 'H2O ') 'NATU' 'DISCRET') 'ET'
  204. ('MANUEL' 'CHPO' MCEN 1 'N2'
  205. (PGAS . 'CP' . 'N2 ') 'NATU' 'DISCRET') ;
  206. *** CV ****
  207. CV = ('MANUEL' 'CHPO' MCEN 1 'H2'
  208. (PGAS . 'CV' . 'H2 ') 'NATU' 'DISCRET') 'ET'
  209. ('MANUEL' 'CHPO' MCEN 1 'O2'
  210. (PGAS . 'CV' . 'O2 ') 'NATU' 'DISCRET') 'ET'
  211. ('MANUEL' 'CHPO' MCEN 1 'H2O'
  212. (PGAS . 'CV' . 'H2O ') 'NATU' 'DISCRET') 'ET'
  213. ('MANUEL' 'CHPO' MCEN 1 'N2'
  214. (PGAS . 'CV' . 'N2 ') 'NATU' 'DISCRET') ;
  215. *** mass fractions ****
  216. FF = 'BRUI' 'BLAN' 'UNIF' 0.2 0.05 MCEN ;
  217. CHF = 'NOMC' 'H2' FF 'NATU' 'DISCRET' ;
  218.  
  219. FS = 'BRUI' 'BLAN' 'UNIF' 0.3 0.04 MCEN ;
  220. CHS = 'NOMC' 'O2' FS 'NATU' 'DISCRET' ;
  221.  
  222. FT = 'BRUI' 'BLAN' 'UNIF' 0.15 0.045 MCEN ;
  223. CHT = 'NOMC' 'H2O' FT 'NATU' 'DISCRET' ;
  224.  
  225. YNB = CHF 'ET' CHS ;
  226. YN = YNB 'ET' CHT ;
  227.  
  228. UNCH = 'MANUEL' 'CHPO' MCEN 1 'SCAL' 1.0 ;
  229.  
  230. YNESP = UNCH '-' ('PSCAL' UNCH YN ('MOTS' 'SCAL' 'SCAL' 'SCAL')
  231. (PGAS . 'ESPEULE')) ;
  232. YNESP = 'NOMC' 'N2 ' YNESP 'NATU' 'DISCRET' ;
  233.  
  234. YNTOT = YN 'ET' YNESP ;
  235.  
  236. RYN = YN '*' RN (PGAS . 'ESPEULE') ('MOTS' 'SCAL' 'SCAL' 'SCAL')
  237. (PGAS . 'ESPEULE') ;
  238. *** gamma ***************
  239. LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
  240. CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
  241. CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
  242. GAMN = CPTOT '/' CVTOT ;
  243. **** velocities **********
  244. CSONN = (GAMN '*' PN) '/' RN ;
  245. CSONN = 'KOPS' CSONN '**' 0.5 ;
  246. UXN = -0.005 * CSONN ;
  247. UYN = 0.002 * CSONN ;
  248. ***************************************************
  249. UINF = 'MANUEL' 'CHPO' MCEN 1 'SCAL' 0.001
  250. 'NATU' 'DISCRET' ;
  251. UPRI = 'MANUEL' 'CHPO' MCEN 1 'SCAL' 0.0001
  252. 'NATU' 'DISCRET' ;
  253.  
  254. ****************************************************************
  255. ****************************************************************
  256. ****************************************************************
  257.  
  258. GN = ('NOMC' (RN '*' UXN) 'UX' 'NATU' 'DISCRET') 'ET'
  259. <$iv sTyld="font: normal normal 1em/1.2em monospace; margin:0; padding:0; background:none; vertical-align:top;"> ('NOMC' (RN '*' UYN) 'UY' 'NATU' 'DISCRET') ;
  •  
  •  
  •  
  • ECIN = 0.5D0 '*' RN '*' ((UXN '*' UXN) '+' (UYN '*' UYN)) ;
  • RETN = 'NOMC' 'SCAL' ((PN '/' (GAMN '-' 1.0)) '+' ECIN)
  • 'NATU' 'DISCRET' ;
  •  
  • VIT PRES TTEMP FRAM GAMM = 'PRIM' 'PERFMULT' PGAS RN GN RETN RYN ;
  • ****************************************************************
  • ********* Application of the Operator 'KONV' ******************
  • ****************************************************************
  • NOMDEN = 'RN ' ;
  • NOMMOX = 'RUXN' ;
  • NOMMOY = 'RUYN' ;
  • NOMRET = 'RETN' ;
  • NOMYH2 = 'RYH2' ;
  • NOMYO2 = 'RYO2' ;
  • NOMYHO = 'RYHO' ;
  •  
  • LINCO = 'MOTS' NOMDEN NOMMOX NOMMOY NOMRET
  • NOMYH2 NOMYO2 NOMYHO ;
  • ******************************************************************
  • MRN MUN MPN MYN MGAMN = 'PRET' 'PERFMULT' 1 1 $DOMTOT
  • RN VIT PN YN GAMN ;
  •  
  • RRESID RFLOT = 'KONV' 'VF' 'PERFMULT' 'RESI' METO LINCO
  • $DOMTOT MRN MUN MPN MGAMN MYN PGAS UINF UPRI MELLIM ;
  •  
  • ******************************************************************
  • ******************************************************************
  •  
  • RMAT = 'KONV' 'VF' 'PERFMULT' 'JACOCONS' METO $DOMTOT PGAS LINCO
  • MELLIM RN VIT PN YN UINF UPRI ;
  •  
  • ******************************************************************
  •  
  • PCEN9 = 'POIN' 1 ('DOMA' $DOM9 CENTRE) ;
  •  
  • *****************************************************************
  • **** Beggining of the procedure **********
  • *****************************************************************
  •  
  • 'DEBPROC' JAEXAC RMAT*'MATRIK' $DOMA*'MMODEL' LISTINCO*'LISTMOTS'
  • PPRIM*'POINT' PDUAL*'POINT' MOTPRI*'MOT' MOTDUA*'MOT' ;
  •  
  • MCEN = 'DOMA' $DOMA CENTRE ;
  • ELT1 = 'MANUEL' 'POI1' PPRIM ;
  • NDIM = 'DIME' LISTINCO ;
  • CHPUN = 'MANUEL' 'CHPO' ELT1 1 MOTPRI 1.0 'NATURE' 'DISCRET' ;
  • **********************
  • 'REPETER' BL1 NDIM ;
  • MOTCEL = 'EXTRAIRE' LISTINCO &BL1 ;
  • CHPUN = CHPUN 'ET' ('MANUEL' 'CHPO' MCEN 1 MOTCEL
  • 0.0 'NATURE' 'DISCRET') ;
  • 'FIN' BL1 ;
  • *********************
  • CHCH = 'KOPS' RMAT 'MULT' CHPUN ;
  • SCAL = 'EXTRAIRE' CHCH PDUAL MOTDUA ;
  • 'FINPROC' SCAL ;
  •  
  • ****************************************************************
  • ***** End of the procedure JAEXAC ********
  • ****************************************************************
  •  
  • ****************************************************************
  • **************** Procedure JACNUM ************************
  • ****************************************************************
  • 'DEBPROC' JACNUM $DOMA*'MMODEL' NOMMET*'MOT' RN*'CHPOINT'
  • GN*'CHPOINT' RETN*'CHPOINT' RYN*'CHPOINT' GAMN*'CHPOINT'
  • PGAS*'TABLE ' LISTINCO*'LISTMOTS' MELLIM*'MAILLAGE'
  • UINF*'CHPOINT' UPRI*'CHPOINT'
  • PPRIM*'POINT' PDUAL*'POINT'
  • MOTPRI*'MOT' MOTDUA*'MOT' EPSILON*'FLOTTANT';
  •  
  •  
  • RNCEL = 'COPIER' RN ;
  • GNCEL = 'COPIER' GN ;
  • RETNCEL = 'COPIER' RETN ;
  • RYNCEL = 'COPIER' RYN ;
  •  
  • VIT PRES TTEMP FRAM GAMM = 'PRIM' 'PERFMULT'
  • PGAS RNCEL GNCEL RETNCEL RYNCEL ;
  •  
  • MRN MUN MPN MYN MGAMN = 'PRET' 'PERFMULT' 1 1 $DOMTOT
  • RNCEL VIT PRES FRAM GAMN ;
  •  
  • CHPRESO RFLOT = 'KONV' 'VF' 'PERFMULT' 'RESI' NOMMET
  • LISTINCO $DOMA MRN MUN MPN MGAMN MYN PGAS UINF UPRI ;
  •  
  • VAL0 = 'EXTRAIRE' CHPRESO PDUAL MOTDUA ;
  • ******************************************
  • dens0 = 'EXTRAIRE' RN PPRIM 'SCAL' ;
  •  
  • CN2 = GAMN '*' (PRES '/' RN) ;
  •  
  • cson0 = ('EXTRAIRE' CN2 PPRIM 'SCAL') '**' 0.5 ;
  •  
  • ret0 = ('EXTRAIRE' RETN PPRIM 'SCAL') ;
  • ******************************************
  • NDIM = 'DIME' LISTINCO ;
  • 'REPETER' BL1 NDIM ;
  • MOTCEL = 'EXTRAIRE' LISTINCO &BL1 ;
  • 'SI' ('EGA' MOTCEL MOTPRI) ;
  • ICEL = &BL1 ;
  • 'QUITTER' BL1 ;
  • 'FINSI' ;
  • 'FIN' BL1 ;
  •  
  • 'SI' (ICEL > NDIM) ;
  • 'MESSAGE' 'Procedure JACNUM' ;
  • 'MESSAGE' 'MOTPRI = ??? ';
  • 'ERREUR' 21 ;
  • 'FINSI' ;
  •  
  • ELT1 = 'MANUEL' 'POI1' PPRIM ;
  • ****************************************************
  • * ICEL = 1 -> On perturbe la densité
  • ****************************************************
  • 'SI' ('EGA' ICEL 1) ;
  • DELTATOT = (EPSILON * dens0) ;
  • RNCEL = ('MANUEL' 'CHPO' ELT1 1 'SCAL' DELTATOT
  • 'NATURE' 'DISCRET') 'ET' RN ;
  • GNCEL = 'COPIER' GN ;
  • RETNCEL = 'COPIER' RETN ;
  • RYNCEL = 'COPIER' RYN ;
  • 'FINSI' ;
  • ****************************************************
  • * ICEL = 2 -> On perturbe la q.d.m. long l'ax x
  • ****************************************************
  • 'SI' ('EGA' ICEL 2) ;
  • DELTATOT = (EPSILON * dens0 * cson0) ;
  • GNCEL = ('MANUEL' 'CHPO' ELT1 1 'UX' DELTATOT
  • 'NATURE' 'DISCRET') 'ET' GN ;
  • RNCEL = 'COPIER' RN ;
  • RETNCEL = 'COPIER' RETN ;
  • RYNCEL = 'COPIER' RYN ;
  • 'FINSI' ;
  • ****************************************************
  • * ICEL = 3 -> On perturbe la q.d.m. long l'ax y
  • ****************************************************
  • 'SI' ('EGA' ICEL 3) ;
  • DELTATOT = (EPSILON * dens0 * cson0) ;
  • GNCEL = ('MANUEL' 'CHPO' ELT1 1 'UY' DELTATOT
  • 'NATURE' 'DISCRET') 'ET' GN ;
  • RNCEL = 'COPIER' RN ;
  • RETNCEL = 'COPIER' RETN ;
  • RYNCEL = 'COPIER' RYN ;
  • 'FINSI' ;
  • ****************************************************
  • * ICEL = 4 -> On perturbe l'énergie totale
  • ****************************************************
  • 'SI' ('EGA' ICEL 4) ;
  • DELTATOT = (EPSILON '*' ret0) ;
  • RETNCEL = ('MANUEL' 'CHPO' ELT1 1 'SCAL' DELTATOT
  • 'NATURE' 'DISCRET') 'ET' RETN ;
  • RNCEL = 'COPIER' RN ;
  • GNCEL = 'COPIER' GN ;
  • RYNCEL = 'COPIER' RYN ;
  • 'FINSI' ;
  • ****************************************************
  • * ICEL = 5 -> On perturbe (\rho Y_{H2})
  • ****************************************************
  • 'SI' ('EGA' ICEL 5) ;
  • DELTATOT = (EPSILON * dens0 ) ;
  • RYNCEL = ('MANUEL' 'CHPO' ELT1 1 'H2' DELTATOT
  • 'NATURE' 'DISCRET') 'ET' RYN ;
  • RNCEL = 'COPIER' RN ;
  • GNCEL = 'COPIER' GN ;
  • RETNCEL = 'COPIER' RETN ;
  • 'FINSI' ;
  • ****************************************************
  • * ICEL = 6 -> On perturbe (\rho Y_{O2})
  • ****************************************************
  • 'SI' ('EGA' ICEL 6) ;
  • DELTATOT = (EPSILON * dens0) ;
  • RYNCEL = ('MANUEL' 'CHPO' ELT1 1 'O2' DELTATOT
  • 'NATURE' 'DISCRET') 'ET' RYN ;
  • RNCEL = 'COPIER' RN ;
  • GNCEL = 'COPIER' GN ;
  • RETNCEL = 'COPIER' RETN ;
  • 'FINSI' ;
  • ****************************************************
  • * ICEL = 7 -> On perturbe (\rho Y_{HO})
  • ****************************************************
  • 'SI' ('EGA' ICEL 7) ;
  • DELTATOT = (EPSILON * dens0) ;
  • RYNCEL = ('MANUEL' 'CHPO' ELT1 1 'H2O' DELTATOT
  • 'NATURE' 'DISCRET') 'ET' RYN ;
  • RNCEL = 'COPIER' RN ;
  • GNCEL = 'COPIER' GN ;
  • RETNCEL = 'COPIER' RETN ;
  • 'FINSI' ;
  • ****************************************************
  • VIT PRES TTEMP FRAM GAMM = 'PRIM' 'PERFMULT'
  • PGAS RNCEL GNCEL RETNCEL RYNCEL ;
  •  
  • MRNF MUNF MPNF MYNF MGAMNF = 'PRET' 'PERFMULT' 1 1 $DOMTOT
  • RNCEL VIT PRES FRAM GAMM ;
  •  
  • CHPRES1 RFLOT = 'KONV' 'VF' 'PERFMULT' 'RESI' NOMMET
  • LISTINCO $DOMA MRNF MUNF MPNF MGAMNF MYNF PGAS UINF UPRI MELLIM ;
  •  
  • VAL1 = 'EXTRAIRE' CHPRES1 PDUAL MOTDUA ;
  • aaa = (((VAL1 '-' VAL0) '/' 100000.0) '/' DELTATOT ) ;
  • * 'LISTE' aaa ;
  • * 'LISTE' ((VAL1 '-' VAL0) '/' DELTATOT) ;
  • * 'MESSAGE' (CHAINE 'VAL1=' FORMAT '(F17.9)' VAL1) ;
  • * 'MESSAGE' (CHAINE 'VAL0=' FORMAT '(F17.9)' VAL0) ;
  • * 'LISTE' ret0 ;
  • * 'LISTE' DELTATOT ;
  • * 'LISTE' ICEL ;
  • 'FINPROC' ((VAL1 '-' VAL0) '/' DELTATOT) ;
  • ****************************************************************
  • ****************************************************************
  • DRR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMDEN ;
  • DGXR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMMOX ;
  • DGYR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMMOY ;
  • DRETR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMRET ;
  • DYH2R = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMYH2 ;
  • DYO2R = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMYO2 ;
  • DYHOR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMYHO ;
  • ***************************************
  • DRGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMDEN ;
  • DGXGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMMOX ;
  • DGYGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMMOY ;
  • DRETGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMRET ;
  • DYH2GX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMYH2 ;
  • DYO2GX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMYO2 ;
  • DYHOGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMYHO ;
  • ***************************************
  • DRGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMDEN ;
  • DGXGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMMOX ;
  • DGYGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMMOY ;
  • DRETGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMRET ;
  • DYH2GY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMYH2 ;
  • DYO2GY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMYO2 ;
  • DYHOGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMYHO ;
  • ***************************************
  • DRRET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMDEN ;
  • DGXRET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMMOX ;
  • DGYRET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMMOY ;
  • DRETRET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMRET ;
  • DYH2RET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMYH2 ;
  • DYO2RET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMYO2 ;
  • DYHORET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMYHO ;
  • ***************************************
  • DRYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMDEN ;
  • DGXYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMMOX ;
  • DGYYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMMOY ;
  • DRETYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMRET ;
  • DYH2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMYH2 ;
  • DYO2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMYO2 ;
  • DYHOYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMYHO ;
  • ***************************************
  • DRYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMDEN ;
  • DGXYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMMOX ;
  • DGYYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMMOY ;
  • DRETYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMRET ;
  • DYH2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMYH2 ;
  • DYO2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMYO2 ;
  • DYHOYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMYHO ;
  • ****************************************
  • DRYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMDEN ;
  • DGXYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMMOX ;
  • DGYYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMMOY ;
  • DRETYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMRET ;
  • DYH2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMYH2 ;
  • DYO2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMYO2 ;
  • DYHOYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMYHO ;
  • ****************************************************************
  • *****************************************************************
  • DRRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMDEN NOMDEN DELTA ;
  • DGXRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMDEN NOMMOX DELTA ;
  • DGYRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMDEN NOMMOY DELTA ;
  • DRETRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMDEN NOMRET DELTA ;
  • DYH2RN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMDEN NOMYH2 DELTA ;
  • DYO2RN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMDEN NOMYO2 DELTA ;
  • DYHORN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMDEN NOMYHO DELTA ;
  • **************************************
  • DRGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMMOX NOMDEN DELTA ;
  • DGXGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMMOX NOMMOX DELTA ;
  • DGYGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMMOX NOMMOY DELTA ;
  • DRETGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMMOX NOMRET DELTA ;
  • DYH2GXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMMOX NOMYH2 DELTA ;
  • DYO2GXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMMOX NOMYO2 DELTA ;
  • DYHOGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMMOX NOMYHO DELTA ;
  • **************************************
  • DRGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMMOY NOMDEN DELTA ;
  • DGXGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMMOY NOMMOX DELTA ;
  • DGYGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMMOY NOMMOY DELTA ;
  • DRETGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMMOY NOMRET DELTA ;
  • DYH2GYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMMOY NOMYH2 DELTA ;
  • DYO2GYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMMOY NOMYO2 DELTA ;
  • DYHOGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMMOY NOMYHO DELTA ;
  • **************************************
  • DRRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMRET NOMDEN DELTA ;
  • DGXRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMRET NOMMOX DELTA ;
  • DGYRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMRET NOMMOY DELTA ;
  • DRETRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMRET NOMRET DELTA ;
  • DYH2RETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMRET NOMYH2 DELTA ;
  • DYO2RETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMRET NOMYO2 DELTA ;
  • DYHORETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMRET NOMYHO DELTA ;
  • **************************************
  • DRYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYH2 NOMDEN DELTA ;
  • DGXYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYH2 NOMMOX DELTA ;
  • DGYYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYH2 NOMMOY DELTA ;
  • DRETYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYH2 NOMRET DELTA ;
  • DYH2YH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYH2 NOMYH2 DELTA ;
  • DYO2YH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYH2 NOMYO2 DELTA ;
  • DYHOYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYH2 NOMYHO DELTA ;
  • ***************************************
  • DRYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYO2 NOMDEN DELTA ;
  • DGXYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYO2 NOMMOX DELTA ;
  • DGYYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYO2 NOMMOY DELTA ;
  • DRETYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYO2 NOMRET DELTA ;
  • DYH2YO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYO2 NOMYH2 DELTA ;
  • DYO2YO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYO2 NOMYO2 DELTA ;
  • DYHOYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYO2 NOMYHO DELTA ;
  • ***************************************
  • DRYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYHO NOMDEN DELTA ;
  • DGXYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYHO NOMMOX DELTA ;
  • DGYYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYHO NOMMOY DELTA ;
  • DRETYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYHO NOMRET DELTA ;
  • DYH2YHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYHO NOMYH2 DELTA ;
  • DYO2YHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYHO NOMYO2 DELTA ;
  • DYHOYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN9 PCEN9 NOMYHO NOMYHO DELTA ;
  • **************************************************************
  • * Test des comparaisons jacobien exact-jacobien numerique
  • **************************************************************
  • 'SI' (('ABS'(DRR)) < DELTA) ;
  • ERR1 = ('ABS' (DRRN '-' DRR)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRRN '-' DRR)) '/' ('ABS' (DRR)) ;
  • 'FINSI' ;
  • ***********************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 1a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DGXR)) < DELTA) ;
  • ERR1 = ('ABS' (DGXRN '-' DGXR)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGXRN '-' DGXR)) '/' ('ABS' (DGXR)) ;
  • 'FINSI' ;
  • ***********************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 2a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DGYR)) < DELTA) ;
  • ERR1 = ('ABS' (DGYRN '-' DGYR)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGYRN '-' DGYR)) '/' ('ABS' (DGYR)) ;
  • 'FINSI' ;
  • ************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 3a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DRETR)) < DELTA) ;
  • ERR1 = ('ABS' (DRETRN '-' DRETR)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRETRN '-' DRETR)) '/' ('ABS' (DRETR)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 4a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYH2R)) < DELTA) ;
  • ERR1 = ('ABS' (DYH2RN '-' DYH2R)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYH2RN '-' DYH2R)) '/' ('ABS' (DYH2R)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 5a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYO2R)) < DELTA) ;
  • ERR1 = ('ABS' (DYO2RN '-' DYO2R)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYO2RN '-' DYO2R)) '/' ('ABS' (DYO2R)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 6a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYHOR)) < DELTA) ;
  • ERR1 = ('ABS' (DYHORN '-' DYHOR)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYHORN '-' DYHOR)) '/' ('ABS' (DYHOR)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 7a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • ***** Second line !!!!!!!!!!!!!!!!!!!!!!! **************
  • **************************************************************
  • **************************************************************
  • 'SI' (('ABS'(DRGX)) < DELTA) ;
  • ERR1 = ('ABS' (DRGXN '-' DRGX)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRGXN '-' DRGX)) '/' ('ABS' (DRGX)) ;
  • 'FINSI' ;
  • ***************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 8a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DGXGX)) < DELTA) ;
  • ERR1 = ('ABS' (DGXGXN '-' DGXGX)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGXGXN '-' DGXGX)) '/' ('ABS' (DGXGX)) ;
  • 'FINSI' ;
  • *****************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 9a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DGYGX)) < DELTA) ;
  • ERR1 = ('ABS' (DGYGXN '-' DGYGX)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGYGXN '-' DGYGX)) '/' ('ABS' (DGYGX)) ;
  • 'FINSI' ;
  • ******************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 10a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DRETGX)) < DELTA) ;
  • ERR1 = ('ABS' (DRETGXN '-' DRETGX)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRETGXN '-' DRETGX)) '/' ('ABS' (DRETGX)) ;
  • 'FINSI' ;
  • *******************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 11a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYH2GX)) < DELTA) ;
  • ERR1 = ('ABS' (DYH2GXN '-' DYH2GX)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYH2GXN '-' DYH2GX)) '/' ('ABS' (DYH2GX)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 12a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYO2GX)) < DELTA) ;
  • ERR1 = ('ABS' (DYO2GXN '-' DYO2GX)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYO2GXN '-' DYO2GX)) '/' ('ABS' (DYO2GX)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 13a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYHOGX)) < DELTA) ;
  • ERR1 = ('ABS' (DYHOGXN '-' DYHOGX)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYHOGXN '-' DYHOGX)) '/' ('ABS' (DYHOGX)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 14a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • ***** Third line !!!!!!!!!!!!!!!!!!!!!!! **************
  • **************************************************************
  • **************************************************************
  • 'SI' (('ABS'(DRGY)) < DELTA) ;
  • ERR1 = ('ABS' (DRGYN '-' DRGY)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRGYN '-' DRGY)) '/' ('ABS' (DRGY)) ;
  • 'FINSI' ;
  • *******************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 15a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DGXGY)) < DELTA) ;
  • ERR1 = ('ABS' (DGXGYN '-' DGXGY)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGXGYN '-' DGXGY)) '/' ('ABS' (DGXGY)) ;
  • 'FINSI' ;
  • ********************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 16a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DGYGY)) < DELTA) ;
  • ERR1 = ('ABS' (DGYGYN '-' DGYGY)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGYGYN '-' DGYGY)) '/' ('ABS' (DGYGY)) ;
  • 'FINSI' ;
  • ***************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 17a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DRETGY)) < DELTA) ;
  • ERR1 = ('ABS' (DRETGYN '-' DRETGY)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRETGYN '-' DRETGY)) '/' ('ABS' (DRETGY)) ;
  • 'FINSI' ;
  • ****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 18a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYH2GY)) < DELTA) ;
  • ERR1 = ('ABS' (DYH2GYN '-' DYH2GY)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYH2GYN '-' DYH2GY)) '/' ('ABS' (DYH2GY)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 19a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYO2GY)) < DELTA) ;
  • ERR1 = ('ABS' (DYO2GYN '-' DYO2GY)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYO2GYN '-' DYO2GY)) '/' ('ABS' (DYO2GY)) ;
  • 'FINSI' ;
  • * 'OPTION' DONN 5 ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 20a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYHOGY)) < DELTA) ;
  • ERR1 = ('ABS' (DYHOGYN '-' DYHOGY)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYHOGYN '-' DYHOGY)) '/' ('ABS' (DYHOGY)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 21a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • ******* Fouth line!!!!!!!!!!!!!!!!!!! ******************
  • **************************************************************
  • **************************************************************
  • 'SI' (('ABS'(DRRET)) < DELTA) ;
  • ERR1 = ('ABS' (DRRETN '-' DRRET)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRRETN '-' DRRET)) '/' ('ABS' (DRRET)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 22a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DGXRET)) < DELTA) ;
  • ERR1 = ('ABS' (DGXRETN '-' DGXRET)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGXRETN '-' DGXRET)) '/' ('ABS' (DGXRET)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 23a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DGYRET)) < DELTA) ;
  • ERR1 = ('ABS' (DGYRETN '-' DGYRET)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGYRETN '-' DGYRET)) '/' ('ABS' (DGYRET)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 24a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DRETRET)) < DELTA) ;
  • ERR1 = ('ABS' (DRETRETN '-' DRETRET)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRETRETN '-' DRETRET)) '/' ('ABS' (DRETRET)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 25a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYH2RET)) < DELTA) ;
  • ERR1 = ('ABS' (DYH2RETN '-' DYH2RET)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYH2RETN '-' DYH2RET)) '/' ('ABS' (DYH2RET)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 26a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYO2RET)) < DELTA) ;
  • ERR1 = ('ABS' (DYO2RETN '-' DYO2RET)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYO2RETN '-' DYO2RET)) '/' ('ABS' (DYO2RET)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 27a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYHORET)) < DELTA) ;
  • ERR1 = ('ABS' (DYHORETN '-' DYHORET)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYHORETN '-' DYHORET)) '/' ('ABS' (DYHORET)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 28a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • ******* Fith line!!!!!!!!!!!!!!!!!!! ******************
  • **************************************************************
  • 'SI' (('ABS'(DRYH2)) < DELTA) ;
  • ERR1 = ('ABS' (DRYH2N '-' DRYH2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRYH2N '-' DRYH2)) '/' ('ABS' (DRYH2)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 29a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DGXYH2)) < DELTA) ;
  • ERR1 = ('ABS' (DGXYH2N '-' DGXYH2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGXYH2N '-' DGXYH2)) '/' ('ABS' (DGXYH2)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 30a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DGYYH2)) < DELTA) ;
  • ERR1 = ('ABS' (DGYYH2N '-' DGYYH2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGYYH2N '-' DGYYH2)) '/' ('ABS' (DGYYH2)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 31a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DRETYH2)) < DELTA) ;
  • ERR1 = ('ABS' (DRETYH2N '-' DRETYH2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRETYH2N '-' DRETYH2)) '/' ('ABS' (DRETYH2)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 32a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYH2YH2)) < DELTA) ;
  • ERR1 = ('ABS' (DYH2YH2N '-' DYH2YH2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYH2YH2N '-' DYH2YH2)) '/' ('ABS' (DYH2YH2)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 33a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYO2YH2)) < DELTA) ;
  • ERR1 = ('ABS' (DYO2YH2N '-' DYO2YH2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYO2YH2N '-' DYO2YH2)) '/' ('ABS' (DYO2YH2)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 34a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYHOYH2)) < DELTA) ;
  • ERR1 = ('ABS' (DYHOYH2N '-' DYHOYH2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYHOYH2N '-' DYHOYH2)) '/' ('ABS' (DYHOYH2)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 35a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • ********* Sixth line !!!!!!!!!!!!!!!!!!!!!! ***************
  • **************************************************************
  • 'SI' (('ABS'(DRYO2)) < DELTA) ;
  • ERR1 = ('ABS' (DRYO2N '-' DRYO2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRYO2N '-' DRYO2)) '/' ('ABS' (DRYO2)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 36a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DGXYO2)) < DELTA) ;
  • ERR1 = ('ABS' (DGXYO2N '-' DGXYO2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGXYO2N '-' DGXYO2)) '/' ('ABS' (DGXYO2)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 37a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DGYYO2)) < DELTA) ;
  • ERR1 = ('ABS' (DGYYO2N '-' DGYYO2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGYYO2N '-' DGYYO2)) '/' ('ABS' (DGYYO2)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 38a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DRETYO2)) < DELTA) ;
  • ERR1 = ('ABS' (DRETYO2N '-' DRETYO2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRETYO2N '-' DRETYO2)) '/' ('ABS' (DRETYO2)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 39a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYH2YO2)) < DELTA) ;
  • ERR1 = ('ABS' (DYH2YO2N '-' DYH2YO2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYH2YO2N '-' DYH2YO2)) '/' ('ABS' (DYH2YO2)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 40a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYO2YO2)) < DELTA) ;
  • ERR1 = ('ABS' (DYO2YO2N '-' DYO2YO2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYO2YO2N '-' DYO2YO2)) '/' ('ABS' (DYO2YO2)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 41a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYHOYO2)) < DELTA) ;
  • ERR1 = ('ABS' (DYHOYO2N '-' DYHOYO2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYHOYO2N '-' DYHOYO2)) '/' ('ABS' (DYHOYO2)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 42a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • ********* Seventh line !!!!!!!!!!!!!!!!!!!!!! *************
  • **************************************************************
  • 'SI' (('ABS'(DRYHO)) < DELTA) ;
  • ERR1 = ('ABS' (DRYHON '-' DRYHO)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRYHON '-' DRYHO)) '/' ('ABS' (DRYHO)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 43a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DGXYHO)) < DELTA) ;
  • ERR1 = ('ABS' (DGXYHON '-' DGXYHO)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGXYHON '-' DGXYHO)) '/' ('ABS' (DGXYHO)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 44a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DGYYHO)) < DELTA) ;
  • ERR1 = ('ABS' (DGYYHON '-' DGYYHO)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGYYHON '-' DGYYHO)) '/' ('ABS' (DGYYHO)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 45a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DRETYHO)) < DELTA) ;
  • ERR1 = ('ABS' (DRETYHON '-' DRETYHO)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRETYHON '-' DRETYHO)) '/' ('ABS' (DRETYHO)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 46a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYH2YHO)) < DELTA) ;
  • ERR1 = ('ABS' (DYH2YHON '-' DYH2YHO)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYH2YHON '-' DYH2YHO)) '/' ('ABS' (DYH2YHO)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 47a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYO2YHO)) < DELTA) ;
  • ERR1 = ('ABS' (DYO2YHON '-' DYO2YHO)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYO2YHON '-' DYO2YHO)) '/' ('ABS' (DYO2YHO)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 48a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYHOYHO)) < DELTA) ;
  • ERR1 = ('ABS' (DYHOYHON '-' DYHOYHO)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYHOYHON '-' DYHOYHO)) '/' ('ABS' (DYHOYHO)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 49a';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *****************************************************
  • *****************************************************
  • *****************************************************
  • *****************************************************
  • ******* TEST2 ***************************************
  • *****************************************************
  • *****************************************************
  • *
  • * On compare le jacobien et la variation du residu
  • * en $DOM9 . 'CENTRE' par rapport à une variation
  • * infinitésimal en $DOM6 . 'CENTRE'
  • *
  •  
  • PCEN6 = 'POIN' 1 ('DOMA' $DOM6 CENTRE) ;
  •  
  • ****************************************************
  • * Le jacobien exact.
  • ****************************************************
  • DRR = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMDEN NOMDEN ;
  • DGXR = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMDEN NOMMOX ;
  • DGYR = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMDEN NOMMOY ;
  • DRETR = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMDEN NOMRET ;
  • DYH2R = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMDEN NOMYH2 ;
  • DYO2R = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMDEN NOMYO2 ;
  • DYHOR = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMDEN NOMYHO ;
  • ***************************************
  • DRGX = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOX NOMDEN ;
  • DGXGX = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOX NOMMOX ;
  • DGYGX = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOX NOMMOY ;
  • DRETGX = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOX NOMRET ;
  • DYH2GX = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOX NOMYH2 ;
  • DYO2GX = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOX NOMYO2 ;
  • DYHOGX = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOX NOMYHO ;
  • ***************************************
  • DRGY = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOY NOMDEN ;
  • DGXGY = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOY NOMMOX ;
  • DGYGY = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOY NOMMOY ;
  • DRETGY = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOY NOMRET ;
  • DYH2GY = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOY NOMYH2 ;
  • DYO2GY = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOY NOMYO2 ;
  • DYHOGY = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOY NOMYHO ;
  • ***************************************
  • DRRET = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMRET NOMDEN ;
  • DGXRET = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMRET NOMMOX ;
  • DGYRET = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMRET NOMMOY ;
  • DRETRET = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMRET NOMRET ;
  • DYH2RET = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMRET NOMYH2 ;
  • DYO2RET = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMRET NOMYO2 ;
  • DYHORET = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMRET NOMYHO ;
  • ***************************************
  • DRYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYH2 NOMDEN ;
  • DGXYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYH2 NOMMOX ;
  • DGYYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYH2 NOMMOY ;
  • DRETYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYH2 NOMRET ;
  • DYH2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYH2 NOMYH2 ;
  • DYO2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYH2 NOMYO2 ;
  • DYHOYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYH2 NOMYHO ;
  • ***************************************
  • DRYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYO2 NOMDEN ;
  • DGXYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYO2 NOMMOX ;
  • DGYYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYO2 NOMMOY ;
  • DRETYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYO2 NOMRET ;
  • DYH2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYO2 NOMYH2 ;
  • DYO2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYO2 NOMYO2 ;
  • DYHOYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYO2 NOMYHO ;
  • ****************************************
  • DRYHO = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYHO NOMDEN ;
  • DGXYHO = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYHO NOMMOX ;
  • DGYYHO = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYHO NOMMOY ;
  • DRETYHO = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYHO NOMRET ;
  • DYH2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYHO NOMYH2 ;
  • DYO2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYHO NOMYO2 ;
  • DYHOYHO = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYHO NOMYHO ;
  • *****************************************************************
  • DRRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMDEN NOMDEN DELTA ;
  • DGXRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMDEN NOMMOX DELTA ;
  • DGYRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMDEN NOMMOY DELTA ;
  • DRETRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMDEN NOMRET DELTA ;
  • DYH2RN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMDEN NOMYH2 DELTA ;
  • DYO2RN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMDEN NOMYO2 DELTA ;
  • DYHORN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMDEN NOMYHO DELTA ;
  • **************************************
  • DRGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMMOX NOMDEN DELTA ;
  • DGXGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMMOX NOMMOX DELTA ;
  • DGYGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMMOX NOMMOY DELTA ;
  • DRETGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMMOX NOMRET DELTA ;
  • DYH2GXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMMOX NOMYH2 DELTA ;
  • DYO2GXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMMOX NOMYO2 DELTA ;
  • DYHOGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMMOX NOMYHO DELTA ;
  • **************************************
  • DRGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMMOY NOMDEN DELTA ;
  • DGXGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMMOY NOMMOX DELTA ;
  • DGYGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMMOY NOMMOY DELTA ;
  • DRETGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMMOY NOMRET DELTA ;
  • DYH2GYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMMOY NOMYH2 DELTA ;
  • DYO2GYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMMOY NOMYO2 DELTA ;
  • DYHOGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMMOY NOMYHO DELTA ;
  • **************************************
  • DRRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMRET NOMDEN DELTA ;
  • DGXRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMRET NOMMOX DELTA ;
  • DGYRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMRET NOMMOY DELTA ;
  • DRETRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMRET NOMRET DELTA ;
  • DYH2RETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMRET NOMYH2 DELTA ;
  • DYO2RETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMRET NOMYO2 DELTA ;
  • DYHORETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMRET NOMYHO DELTA ;
  • **************************************
  • DRYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYH2 NOMDEN DELTA ;
  • DGXYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYH2 NOMMOX DELTA ;
  • DGYYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYH2 NOMMOY DELTA ;
  • DRETYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYH2 NOMRET DELTA ;
  • DYH2YH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYH2 NOMYH2 DELTA ;
  • DYO2YH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYH2 NOMYO2 DELTA ;
  • DYHOYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYH2 NOMYHO DELTA ;
  • ***************************************
  • DRYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYO2 NOMDEN DELTA ;
  • DGXYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYO2 NOMMOX DELTA ;
  • DGYYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYO2 NOMMOY DELTA ;
  • DRETYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYO2 NOMRET DELTA ;
  • DYH2YO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYO2 NOMYH2 DELTA ;
  • DYO2YO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYO2 NOMYO2 DELTA ;
  • DYHOYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYO2 NOMYHO DELTA ;
  • ***************************************
  • DRYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYHO NOMDEN DELTA ;
  • DGXYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYHO NOMMOX DELTA ;
  • DGYYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYHO NOMMOY DELTA ;
  • DRETYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYHO NOMRET DELTA ;
  • DYH2YHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYHO NOMYH2 DELTA ;
  • DYO2YHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYHO NOMYO2 DELTA ;
  • DYHOYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO MELLIM
  • UINF UPRI PCEN6 PCEN9 NOMYHO NOMYHO DELTA ;
  • **************************************************************
  • * Test des comparaisons jacobien exact-jacobien numerique
  • **************************************************************
  • 'SI' (('ABS'(DRR)) < DELTA) ;
  • ERR1 = ('ABS' (DRRN '-' DRR)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRRN '-' DRR)) '/' ('ABS' (DRR)) ;
  • 'FINSI' ;
  • ***********************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 1b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DGXR)) < DELTA) ;
  • ERR1 = ('ABS' (DGXRN '-' DGXR)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGXRN '-' DGXR)) '/' ('ABS' (DGXR)) ;
  • 'FINSI' ;
  • ***********************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 2b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DGYR)) < DELTA) ;
  • ERR1 = ('ABS' (DGYRN '-' DGYR)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGYRN '-' DGYR)) '/' ('ABS' (DGYR)) ;
  • 'FINSI' ;
  • ************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 3b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DRETR)) < DELTA) ;
  • ERR1 = ('ABS' (DRETRN '-' DRETR)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRETRN '-' DRETR)) '/' ('ABS' (DRETR)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 4b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYH2R)) < DELTA) ;
  • ERR1 = ('ABS' (DYH2RN '-' DYH2R)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYH2RN '-' DYH2R)) '/' ('ABS' (DYH2R)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 5b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYO2R)) < DELTA) ;
  • ERR1 = ('ABS' (DYO2RN '-' DYO2R)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYO2RN '-' DYO2R)) '/' ('ABS' (DYO2R)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 6b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYHOR)) < DELTA) ;
  • ERR1 = ('ABS' (DYHORN '-' DYHOR)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYHORN '-' DYHOR)) '/' ('ABS' (DYHOR)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 7b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • ***** Second line !!!!!!!!!!!!!!!!!!!!!!! **************
  • **************************************************************
  • **************************************************************
  • 'SI' (('ABS'(DRGX)) < DELTA) ;
  • ERR1 = ('ABS' (DRGXN '-' DRGX)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRGXN '-' DRGX)) '/' ('ABS' (DRGX)) ;
  • 'FINSI' ;
  • ***************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 8b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DGXGX)) < DELTA) ;
  • ERR1 = ('ABS' (DGXGXN '-' DGXGX)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGXGXN '-' DGXGX)) '/' ('ABS' (DGXGX)) ;
  • 'FINSI' ;
  • *****************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 9b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DGYGX)) < DELTA) ;
  • ERR1 = ('ABS' (DGYGXN '-' DGYGX)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGYGXN '-' DGYGX)) '/' ('ABS' (DGYGX)) ;
  • 'FINSI' ;
  • ******************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 10b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DRETGX)) < DELTA) ;
  • ERR1 = ('ABS' (DRETGXN '-' DRETGX)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRETGXN '-' DRETGX)) '/' ('ABS' (DRETGX)) ;
  • 'FINSI' ;
  • *******************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 11b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYH2GX)) < DELTA) ;
  • ERR1 = ('ABS' (DYH2GXN '-' DYH2GX)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYH2GXN '-' DYH2GX)) '/' ('ABS' (DYH2GX)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 12b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYO2GX)) < DELTA) ;
  • ERR1 = ('ABS' (DYO2GXN '-' DYO2GX)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYO2GXN '-' DYO2GX)) '/' ('ABS' (DYO2GX)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 13b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYHOGX)) < DELTA) ;
  • ERR1 = ('ABS' (DYHOGXN '-' DYHOGX)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYHOGXN '-' DYHOGX)) '/' ('ABS' (DYHOGX)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 14b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • ***** Third line !!!!!!!!!!!!!!!!!!!!!!! **************
  • **************************************************************
  • **************************************************************
  • 'SI' (('ABS'(DRGY)) < DELTA) ;
  • ERR1 = ('ABS' (DRGYN '-' DRGY)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRGYN '-' DRGY)) '/' ('ABS' (DRGY)) ;
  • 'FINSI' ;
  • *******************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 15b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DGXGY)) < DELTA) ;
  • ERR1 = ('ABS' (DGXGYN '-' DGXGY)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGXGYN '-' DGXGY)) '/' ('ABS' (DGXGY)) ;
  • 'FINSI' ;
  • ********************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 16b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DGYGY)) < DELTA) ;
  • ERR1 = ('ABS' (DGYGYN '-' DGYGY)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGYGYN '-' DGYGY)) '/' ('ABS' (DGYGY)) ;
  • 'FINSI' ;
  • ***************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 17b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DRETGY)) < DELTA) ;
  • ERR1 = ('ABS' (DRETGYN '-' DRETGY)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRETGYN '-' DRETGY)) '/' ('ABS' (DRETGY)) ;
  • 'FINSI' ;
  • ****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 18b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYH2GY)) < DELTA) ;
  • ERR1 = ('ABS' (DYH2GYN '-' DYH2GY)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYH2GYN '-' DYH2GY)) '/' ('ABS' (DYH2GY)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 19b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYO2GY)) < DELTA) ;
  • ERR1 = ('ABS' (DYO2GYN '-' DYO2GY)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYO2GYN '-' DYO2GY)) '/' ('ABS' (DYO2GY)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 20b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYHOGY)) < DELTA) ;
  • ERR1 = ('ABS' (DYHOGYN '-' DYHOGY)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYHOGYN '-' DYHOGY)) '/' ('ABS' (DYHOGY)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 21b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • ******* Fouth line!!!!!!!!!!!!!!!!!!! ******************
  • **************************************************************
  • **************************************************************
  • 'SI' (('ABS'(DRRET)) < DELTA) ;
  • ERR1 = ('ABS' (DRRETN '-' DRRET)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRRETN '-' DRRET)) '/' ('ABS' (DRRET)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 22b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DGXRET)) < DELTA) ;
  • ERR1 = ('ABS' (DGXRETN '-' DGXRET)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGXRETN '-' DGXRET)) '/' ('ABS' (DGXRET)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 23b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DGYRET)) < DELTA) ;
  • ERR1 = ('ABS' (DGYRETN '-' DGYRET)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGYRETN '-' DGYRET)) '/' ('ABS' (DGYRET)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 24b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DRETRET)) < DELTA) ;
  • ERR1 = ('ABS' (DRETRETN '-' DRETRET)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRETRETN '-' DRETRET)) '/' ('ABS' (DRETRET)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 25b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYH2RET)) < DELTA) ;
  • ERR1 = ('ABS' (DYH2RETN '-' DYH2RET)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYH2RETN '-' DYH2RET)) '/' ('ABS' (DYH2RET)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 26b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYO2RET)) < DELTA) ;
  • ERR1 = ('ABS' (DYO2RETN '-' DYO2RET)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYO2RETN '-' DYO2RET)) '/' ('ABS' (DYO2RET)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 27b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYHORET)) < DELTA) ;
  • ERR1 = ('ABS' (DYHORETN '-' DYHORET)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYHORETN '-' DYHORET)) '/' ('ABS' (DYHORET)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 28b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • ******* Fith line!!!!!!!!!!!!!!!!!!! ******************
  • **************************************************************
  • 'SI' (('ABS'(DRYH2)) < DELTA) ;
  • ERR1 = ('ABS' (DRYH2N '-' DRYH2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRYH2N '-' DRYH2)) '/' ('ABS' (DRYH2)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 29b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DGXYH2)) < DELTA) ;
  • ERR1 = ('ABS' (DGXYH2N '-' DGXYH2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGXYH2N '-' DGXYH2)) '/' ('ABS' (DGXYH2)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 30b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DGYYH2)) < DELTA) ;
  • ERR1 = ('ABS' (DGYYH2N '-' DGYYH2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGYYH2N '-' DGYYH2)) '/' ('ABS' (DGYYH2)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 31b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DRETYH2)) < DELTA) ;
  • ERR1 = ('ABS' (DRETYH2N '-' DRETYH2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRETYH2N '-' DRETYH2)) '/' ('ABS' (DRETYH2)) ;
  • 'FINSI' ;
  • *****************
  • * 'SI' (ERR1 > ERRTOL) ;
  • * 'MESSAGE' 'Problem 32b';
  • * 'ERREUR' 5 ;
  • * 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYH2YH2)) < DELTA) ;
  • ERR1 = ('ABS' (DYH2YH2N '-' DYH2YH2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYH2YH2N '-' DYH2YH2)) '/' ('ABS' (DYH2YH2)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 33b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYO2YH2)) < DELTA) ;
  • ERR1 = ('ABS' (DYO2YH2N '-' DYO2YH2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYO2YH2N '-' DYO2YH2)) '/' ('ABS' (DYO2YH2)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 34b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYHOYH2)) < DELTA) ;
  • ERR1 = ('ABS' (DYHOYH2N '-' DYHOYH2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYHOYH2N '-' DYHOYH2)) '/' ('ABS' (DYHOYH2)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 35b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • ********* Sixth line !!!!!!!!!!!!!!!!!!!!!! ***************
  • **************************************************************
  • 'SI' (('ABS'(DRYO2)) < DELTA) ;
  • ERR1 = ('ABS' (DRYO2N '-' DRYO2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRYO2N '-' DRYO2)) '/' ('ABS' (DRYO2)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 36b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DGXYO2)) < DELTA) ;
  • ERR1 = ('ABS' (DGXYO2N '-' DGXYO2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGXYO2N '-' DGXYO2)) '/' ('ABS' (DGXYO2)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 37b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DGYYO2)) < DELTA) ;
  • ERR1 = ('ABS' (DGYYO2N '-' DGYYO2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGYYO2N '-' DGYYO2)) '/' ('ABS' (DGYYO2)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 38b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DRETYO2)) < DELTA) ;
  • ERR1 = ('ABS' (DRETYO2N '-' DRETYO2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRETYO2N '-' DRETYO2)) '/' ('ABS' (DRETYO2)) ;
  • 'FINSI' ;
  • *****************
  • * 'SI' (ERR1 > ERRTOL) ;
  • * 'MESSAGE' 'Problem 39b';
  • * 'ERREUR' 5 ;
  • * 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYH2YO2)) < DELTA) ;
  • ERR1 = ('ABS' (DYH2YO2N '-' DYH2YO2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYH2YO2N '-' DYH2YO2)) '/' ('ABS' (DYH2YO2)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 40b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYO2YO2)) < DELTA) ;
  • ERR1 = ('ABS' (DYO2YO2N '-' DYO2YO2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYO2YO2N '-' DYO2YO2)) '/' ('ABS' (DYO2YO2)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 41b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYHOYO2)) < DELTA) ;
  • ERR1 = ('ABS' (DYHOYO2N '-' DYHOYO2)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYHOYO2N '-' DYHOYO2)) '/' ('ABS' (DYHOYO2)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 42b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • ********* Seventh line !!!!!!!!!!!!!!!!!!!!!! *************
  • **************************************************************
  • 'SI' (('ABS'(DRYHO)) < DELTA) ;
  • ERR1 = ('ABS' (DRYHON '-' DRYHO)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRYHON '-' DRYHO)) '/' ('ABS' (DRYHO)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 43b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DGXYHO)) < DELTA) ;
  • ERR1 = ('ABS' (DGXYHON '-' DGXYHO)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGXYHON '-' DGXYHO)) '/' ('ABS' (DGXYHO)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 44b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *************************************************************
  • 'SI' (('ABS'(DGYYHO)) < DELTA) ;
  • ERR1 = ('ABS' (DGYYHON '-' DGYYHO)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DGYYHON '-' DGYYHO)) '/' ('ABS' (DGYYHO)) ;
  • 'FINSI' ;
  • *****************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 45b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DRETYHO)) < DELTA) ;
  • ERR1 = ('ABS' (DRETYHON '-' DRETYHO)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DRETYHON '-' DRETYHO)) '/' ('ABS' (DRETYHO)) ;
  • 'FINSI' ;
  • *****************
  • * 'SI' (ERR1 > ERRTOL) ;
  • * 'MESSAGE' 'Problem 46b';
  • * 'ERREUR' 5 ;
  • * 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYH2YHO)) < DELTA) ;
  • ERR1 = ('ABS' (DYH2YHON '-' DYH2YHO)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYH2YHON '-' DYH2YHO)) '/' ('ABS' (DYH2YHO)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 47b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYO2YHO)) < DELTA) ;
  • ERR1 = ('ABS' (DYO2YHON '-' DYO2YHO)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYO2YHON '-' DYO2YHO)) '/' ('ABS' (DYO2YHO)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 48b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **************************************************************
  • 'SI' (('ABS'(DYHOYHO)) < DELTA) ;
  • ERR1 = ('ABS' (DYHOYHON '-' DYHOYHO)) ;
  • 'SINON' ;
  • ERR1 = ('ABS' (DYHOYHON '-' DYHOYHO)) '/' ('ABS' (DYHOYHO)) ;
  • 'FINSI' ;
  • **************************
  • 'SI' (ERR1 > ERRTOL) ;
  • 'MESSAGE' 'Problem 49b';
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **********************************************************
  • *****************************************************
  • *****************************************************
  • ******* TEST3 ***************************************
  • *****************************************************
  • *****************************************************
  • *
  • * On observe que la variation du residu
  • * en $DOM9 . 'CENTRE' par rapport à une variation
  • * infinitésimal en $DOM7 . 'CENTRE' doit etre nulle.
  • *
  •  
  • PCEN7 = 'POIN' 1 ('DOMA' $DOM7 CENTRE) ;
  •  
  • *
  • * Le jacobien exact.
  • *
  •  
  • DRR = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMDEN NOMDEN ;
  • DGXR = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMDEN NOMMOX ;
  • DGYR = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMDEN NOMMOY ;
  • DRETR = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMDEN NOMRET ;
  • DYH2R = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMDEN NOMYH2 ;
  • DYO2R = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMDEN NOMYO2 ;
  • DYHOR = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMDEN NOMYHO ;
  • ***************************************
  • DRGX = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOX NOMDEN ;
  • DGXGX = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOX NOMMOX ;
  • DGYGX = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOX NOMMOY ;
  • DRETGX = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOX NOMRET ;
  • DYH2GX = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOX NOMYH2 ;
  • DYO2GX = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOX NOMYO2 ;
  • DYHOGX = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOX NOMYHO ;
  • ***************************************
  • DRGY = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOY NOMDEN ;
  • DGXGY = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOY NOMMOX ;
  • DGYGY = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOY NOMMOY ;
  • DRETGY = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOY NOMRET ;
  • DYH2GY = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOY NOMYH2 ;
  • DYO2GY = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOY NOMYO2 ;
  • DYHOGY = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOY NOMYHO ;
  • ***************************************
  • DRRET = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMRET NOMDEN ;
  • DGXRET = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMRET NOMMOX ;
  • DGYRET = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMRET NOMMOY ;
  • DRETRET = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMRET NOMRET ;
  • DYH2RET = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMRET NOMYH2 ;
  • DYO2RET = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMRET NOMYO2 ;
  • DYHORET = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMRET NOMYHO ;
  • ***************************************
  • DRYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYH2 NOMDEN ;
  • DGXYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYH2 NOMMOX ;
  • DGYYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYH2 NOMMOY ;
  • DRETYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYH2 NOMRET ;
  • DYH2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYH2 NOMYH2 ;
  • DYO2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYH2 NOMYO2 ;
  • DYHOYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYH2 NOMYHO ;
  • ***************************************
  • DRYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYO2 NOMDEN ;
  • DGXYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYO2 NOMMOX ;
  • DGYYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYO2 NOMMOY ;
  • DRETYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYO2 NOMRET ;
  • DYH2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYO2 NOMYH2 ;
  • DYO2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYO2 NOMYO2 ;
  • DYHOYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYO2 NOMYHO ;
  • ****************************************
  • DRYHO = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYHO NOMDEN ;
  • DGXYHO = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYHO NOMMOX ;
  • DGYYHO = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYHO NOMMOY ;
  • DRETYHO = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYHO NOMRET ;
  • DYH2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYHO NOMYH2 ;
  • DYO2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYHO NOMYO2 ;
  • DYHOYHO = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYHO NOMYHO ;
  •  
  • **********************************************
  • 'SI' (('MAXIMUM'
  • ('PROG' DRR DGXR DGYR DRETR
  • DYH2R DYO2R DYHOR
  • DRGX DGXGX DGYGX DRETGX
  • DYH2GX DYO2GX DYHOGX ) 'ABS' )
  • > ERRTOL) ;
  • 'MESSAGE' 'Problem final 1' ;
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *********************************************
  • 'SI' (('MAXIMUM'
  • ('PROG' DRGY DGXGY DGYGY DRETGY
  • DYH2GY DYO2GY DYHOGY
  • DRRET DGXRET DGYRET DRETRET
  • DYH2RET DYO2RET DYHORET ) 'ABS' )
  • > ERRTOL) ;
  • 'MESSAGE' 'Problem final 2' ;
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • *********************************************
  • 'SI' (('MAXIMUM'
  • ('PROG' DRYH2 DGXYH2 DGYYH2 DRETYH2
  • DYH2YH2 DYO2YH2 DYHOYH2
  • DRYO2 DGXYO2 DGYYO2 DRETYO2
  • DYH2YO2 DYO2YO2 DYHOYO2 ) 'ABS' )
  • > ERRTOL) ;
  • 'MESSAGE' 'Problem final 3' ;
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **********************************************
  • 'SI' (('MAXIMUM'
  • ('PROG' DRYHO DGXYHO DGYYHO DRETYHO
  • DYH2YHO DYO2YHO DYHOYHO ) 'ABS' )
  • > ERRTOL) ;
  • 'MESSAGE' 'Problem final 4' ;
  • 'ERREUR' 5 ;
  • 'FINSI' ;
  • **********************************************
  •  
  • 'FIN' ;
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  • © Cast3M 2003 - Tous droits réservés.
    Mentions légales