Télécharger la plaquette Cast3M


fr en

Télécharger konmsp_impl2D.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  konmsp_impl2D.dgibi
  2. ***********************************************************
  3. ***********************************************************
  4. **** APPROCHE VF "Cell-Centred Formulation" pour la    ****
  5. **** solution des                                      ****
  6. **** Equations d'Euler pour un gaz parfait multiespes  ****
  7. **** Implicit: calcul du jacobien du residu            ****
  8. ****                                                   ****
  9. **** Cas gaz multiespes, "calorically perfect"         ****
  10. ****                                                   ****
  11. **** Methodes: AUSM+                                   ****
  12. ****                                                   ****
  13. **** S. Kudriakov  DRN/DMT/SEMT/LTMF    DEC  2002      ****
  14. ***********************************************************
  15. ***********************************************************
  16.  
  17.  'OPTION'  'DIME' 2 ;
  18.  'OPTION'  'ELEM' QUA4 ;
  19.  'OPTION'  'ECHO'  1 ;
  20.  'OPTION'  'TRAC' 'X' ;
  21.  
  22. *
  23. *** GRAPH
  24. *
  25.  
  26.  GRAPH = FAUX ;
  27. * GRAPH = VRAI ;
  28.  
  29.  ERRTOL = 5.0D-4 ;
  30.  DELTA = 0.000001 ;
  31.  METO = 'AUSMPLUS' ;
  32. ***************************
  33. ***** DOMAINE SPATIAL  ****
  34. ***************************
  35.  
  36.  
  37.  A0 = 0.0D0 0.0D0;
  38.  A1 = 1.0D0 0.0D0;
  39.  A2 = 2.0D0 0.0D0;
  40.  A3 = 3.0D0 0.0D0;
  41. ************************ 
  42.  A0A1 = A0 'DROIT' 1 A1;
  43.  A1A2 = A1 'DROIT' 1 A2;
  44.  A2A3 = A2 'DROIT'  1 A3;
  45. ************************
  46.  DOM1 = 'TRANSLATION' A0A1 1 (0.0 1.0) ;
  47.  DOM2 = DOM1 'PLUS' (0.0 1.0) ;
  48.  DOM3 = DOM2 'PLUS' (0.0 1.0) ;
  49.  DOM4 = 'TRANSLATION' A1A2 1 (0.0 1.0) ;
  50.  DOM5 = DOM4 'PLUS' (0.0 1.0) ;
  51.  DOM6 = DOM5 'PLUS' (0.0 1.0) ;
  52.  DOM7 = 'TRANSLATION' A2A3 1 (0.0 1.0) ;
  53.  DOM8 = DOM7 'PLUS' (0.0 1.0) ;
  54.  DOM9 = DOM8 'PLUS' (0.0 1.0) ;
  55.  
  56.  
  57.  ANGLE = 30.0D0 ;
  58.  ORIG1 = 0.0D0 0.0D0 ;
  59.  
  60. 'MESSAGE';
  61. 'MESSAGE' (CHAIN 'ANGLE = ' ANGLE);
  62. 'MESSAGE';
  63.  
  64.  DOM1   = DOM1   'TOURNER' ANGLE ORIG1 ;
  65.  DOM2   = DOM2   'TOURNER' ANGLE ORIG1 ;
  66.  DOM3   = DOM3   'TOURNER' ANGLE ORIG1 ;
  67.  DOM4   = DOM4   'TOURNER' ANGLE ORIG1 ;
  68.  DOM5   = DOM5   'TOURNER' ANGLE ORIG1 ;
  69.  DOM6   = DOM6   'TOURNER' ANGLE ORIG1 ;
  70.  DOM7   = DOM7   'TOURNER' ANGLE ORIG1 ;
  71.  DOM8   = DOM8   'TOURNER' ANGLE ORIG1 ;
  72.  DOM9   = DOM9   'TOURNER' ANGLE ORIG1 ;
  73.  
  74.  
  75.  DOMTOT = DOM1 'ET' DOM2 'ET' DOM3 'ET' DOM4 'ET' DOM5 'ET'
  76.           DOM6 'ET' DOM7 'ET' DOM8 'ET' DOM9 'ELIMINATION' 0.0001 ;
  77.  
  78. ***********************************************************
  79. *** These things are necessary in order to work with    ***
  80. ***             each square separately                  ***
  81. ***********************************************************
  82.  
  83.  $DOMTOT = 'MODELISER' DOMTOT 'EULER';
  84.  
  85.  $DOM1 = 'MODELISER'  DOM1 'EULER';
  86.  $DOM2 = 'MODELISER'  DOM2 'EULER';
  87.  $DOM3 = 'MODELISER'  DOM3 'EULER';
  88.  $DOM4 = 'MODELISER'  DOM4 'EULER';
  89.  $DOM5 = 'MODELISER'  DOM5 'EULER';
  90.  $DOM6 = 'MODELISER'  DOM6 'EULER';
  91.  $DOM7 = 'MODELISER'  DOM7 'EULER';
  92.  $DOM8 = 'MODELISER'  DOM8 'EULER';
  93.  $DOM9 = 'MODELISER'  DOM9 'EULER';
  94.  
  95.  TDOMTOT = 'DOMA'  $DOMTOT 'VF';
  96.  
  97.  TDOM1 = 'DOMA'  $DOM1 'VF';
  98.  TDOM2 = 'DOMA'  $DOM2 'VF';
  99.  TDOM3 = 'DOMA'  $DOM3 'VF';
  100.  TDOM4 = 'DOMA'  $DOM4 'VF';
  101.  TDOM5 = 'DOMA'  $DOM5 'VF';
  102.  TDOM6 = 'DOMA'  $DOM6 'VF';
  103.  TDOM7 = 'DOMA'  $DOM7 'VF';
  104.  TDOM8 = 'DOMA'  $DOM8 'VF';
  105.  TDOM9 = 'DOMA'  $DOM9 'VF';
  106.  
  107.  MDOM1 = TDOM1 . 'QUAF' ;
  108.  MDOM2 = TDOM2 . 'QUAF' ;
  109.  MDOM3 = TDOM3 . 'QUAF' ;
  110.  MDOM4 = TDOM4 . 'QUAF' ;
  111.  MDOM5 = TDOM5 . 'QUAF' ;
  112.  MDOM6 = TDOM6 . 'QUAF' ;
  113.  MDOM7 = TDOM7 . 'QUAF' ;
  114.  MDOM8 = TDOM8 . 'QUAF' ;
  115.  MDOM9 = TDOM9 . 'QUAF' ;          
  116.  
  117. **** old stuff  $DOMTOT = 'DOMA' DOMTOT ;
  118.  
  119.  MDOMTOT = TDOMTOT . 'QUAF' ;
  120.  
  121.  'ELIMINATION'  (MDOMTOT ET MDOM1) 0.0001 ;
  122.  'ELIMINATION'  (MDOMTOT ET MDOM2) 0.0001 ;
  123.  'ELIMINATION'  (MDOMTOT ET MDOM3) 0.0001 ;
  124.  'ELIMINATION'  (MDOMTOT ET MDOM4) 0.0001 ;
  125.  'ELIMINATION'  (MDOMTOT ET MDOM5) 0.0001 ;
  126.  'ELIMINATION'  (MDOMTOT ET MDOM6) 0.0001 ;
  127.  'ELIMINATION'  (MDOMTOT ET MDOM7) 0.0001 ;
  128.  'ELIMINATION'  (MDOMTOT ET MDOM8) 0.0001 ;
  129.  'ELIMINATION'  (MDOMTOT ET MDOM9) 0.0001 ;
  130.  
  131.  
  132.  
  133.  
  134. ********************************************************
  135. *******    MODELE DU GAS   ***************************** 
  136. ********************************************************
  137.  
  138.  NESP = 4;
  139.  
  140. *
  141. *** GAS: H_2, O_2, H_2O, N_2
  142. *
  143. *   CP, CV en J/Kg/K @ T = 3000
  144. *
  145.  
  146.  PGAS = 'TABLE' ;
  147.  
  148.  PGAS . 'CP' = 'TABLE' ;
  149.  PGAS . 'CP' . 'H2  '  = .18729066D+05 ;
  150.  PGAS . 'CP' . 'O2  '  = .11886820D+04 ;
  151.  PGAS . 'CP' . 'H2O '  = .31209047D+04 ;
  152.  PGAS . 'CP' . 'N2  '  = .12993995D+04 ;
  153.  
  154.  
  155.  
  156.  PGAS . 'CV' = 'TABLE' ;
  157.  PGAS . 'CV' . 'H2  '  = .14571861D+05 ;
  158.  PGAS . 'CV' . 'O2  '  = .92885670D+03 ;
  159.  PGAS . 'CV' . 'H2O '  = .26589930D+04 ;
  160.  PGAS . 'CV' . 'N2  '  =  .10024563D+04;
  161.  
  162.  
  163.  
  164. *
  165. **** Especes qui sont dans les equations d'Euler
  166. *
  167.  
  168.  PGAS . 'ESPEULE' = 'MOTS' 'H2  ' 'O2  ' 'H2O ' ;
  169.  
  170. *
  171. **** Espece qui n'y est pas
  172. *
  173.  
  174.  
  175.  PGAS . 'ESPNEULE' = 'N2  ';
  176.  
  177. ******************************************************
  178. ************  Les CHPOINTs   *************************
  179. ******************************************************
  180.  
  181.  MCEN = 'DOMA' $DOMTOT CENTRE ;
  182.  
  183. *** density ** 
  184.  RN = 'BRUI' 'BLAN' 'UNIF' 1.11 0.5 MCEN ;
  185. *** pressure **
  186.  PN = 'BRUI' 'BLAN' 'UNIF' 1234.1 800 MCEN ;
  187. ***  CP  ****
  188.  CP = ('MANUEL' 'CHPO' MCEN 1 'H2'
  189.        (PGAS . 'CP' . 'H2  ') 'NATU' 'DISCRET') 'ET' 
  190.       ('MANUEL' 'CHPO' MCEN 1 'O2'
  191.       (PGAS . 'CP' . 'O2  ') 'NATU' 'DISCRET') 'ET'
  192.        ('MANUEL' 'CHPO' MCEN  1 'H2O'
  193.        (PGAS . 'CP' . 'H2O ') 'NATU' 'DISCRET') 'ET'
  194.        ('MANUEL' 'CHPO' MCEN  1 'N2'
  195.        (PGAS . 'CP' . 'N2  ') 'NATU' 'DISCRET') ;
  196. ***  CV ****
  197.  CV = ('MANUEL' 'CHPO' MCEN 1 'H2'
  198.        (PGAS . 'CV' . 'H2  ') 'NATU' 'DISCRET') 'ET' 
  199.       ('MANUEL' 'CHPO' MCEN 1 'O2'
  200.       (PGAS . 'CV' . 'O2  ') 'NATU' 'DISCRET') 'ET'
  201.        ('MANUEL' 'CHPO' MCEN  1 'H2O'
  202.        (PGAS . 'CV' . 'H2O ') 'NATU' 'DISCRET') 'ET'
  203.        ('MANUEL' 'CHPO' MCEN  1 'N2'
  204.        (PGAS . 'CV' . 'N2  ') 'NATU' 'DISCRET') ;       
  205. *** mass fractions ****
  206.  FF = 'BRUI' 'BLAN' 'UNIF' 0.2 0.05 MCEN  ;
  207.  CHF = 'NOMC' 'H2' FF 'NATU' 'DISCRET' ; 
  208.  
  209.  FS = 'BRUI' 'BLAN' 'UNIF' 0.3 0.04 MCEN ;
  210.  CHS = 'NOMC' 'O2' FS 'NATU' 'DISCRET' ; 
  211.  
  212.  FT = 'BRUI' 'BLAN' 'UNIF' 0.15 0.045 MCEN ;
  213.  CHT = 'NOMC' 'H2O' FT 'NATU' 'DISCRET' ; 
  214.  
  215.  YNB = CHF 'ET' CHS ;
  216.  YN = YNB 'ET' CHT ;
  217.  
  218.  UNCH =  'MANUEL' 'CHPO' MCEN 1 'SCAL' 1.0 ;
  219.  
  220.  YNESP = UNCH '-' ('PSCAL' UNCH YN ('MOTS' 'SCAL' 'SCAL' 'SCAL')
  221.          (PGAS . 'ESPEULE')) ;
  222.  YNESP = 'NOMC' 'N2  ' YNESP 'NATU' 'DISCRET' ;
  223.  
  224.  YNTOT = YN 'ET' YNESP ;
  225.  
  226.  RYN = YN '*' RN  (PGAS . 'ESPEULE') ('MOTS' 'SCAL' 'SCAL' 'SCAL')
  227.        (PGAS . 'ESPEULE') ;
  228. *** gamma  ***************
  229.  LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
  230.  CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
  231.  CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
  232.  GAMN = CPTOT '/' CVTOT ;
  233. **** velocities ********** 
  234.  CSONN = (GAMN '*' PN) '/' RN ;
  235.  CSONN = 'KOPS' CSONN '**' 0.5 ;
  236.  UXN  = 0.5 *  CSONN ;
  237.  UYN  = -0.2 * CSONN ;
  238.  
  239. ****************************************************************
  240. ****************************************************************
  241. ****************************************************************
  242.  
  243.  GN = ('NOMC' (RN '*' UXN) 'UX' 'NATU' 'DISCRET') 'ET'
  244.        ('NOMC' (RN '*' UYN) 'UY' 'NATU' 'DISCRET') ;
  245.  
  246.  
  247.  
  248.  ECIN = 0.5D0 '*' RN '*' ((UXN  '*' UXN) '+' (UYN  '*' UYN)) ;
  249.  RETN  = 'NOMC' 'SCAL' ((PN  '/' (GAMN  '-' 1.0)) '+' ECIN)
  250.    'NATU' 'DISCRET' ;
  251.  
  252.  VIT PRES TTEMP FRAM GAMM  = 'PRIM' 'PERFMULT' PGAS RN GN RETN RYN  ;
  253. ****************************************************************
  254. *********  Application of the Operator 'KONV' ******************
  255. ****************************************************************
  256.  NOMDEN = 'RN  ' ;
  257.  NOMMOX = 'RUXN' ;
  258.  NOMMOY = 'RUYN' ;
  259.  NOMRET = 'RETN' ;
  260.  NOMYH2 = 'RYH2' ;
  261.  NOMYO2 = 'RYO2' ;
  262.  NOMYHO = 'RYHO' ;
  263.  
  264.  LINCO = 'MOTS'  NOMDEN  NOMMOX  NOMMOY  NOMRET
  265.              NOMYH2 NOMYO2  NOMYHO ;
  266. ******************************************************************
  267.  MRN MUN MPN MYN MGAMN = 'PRET' 'PERFMULT' 1 1 $DOMTOT
  268.                     RN VIT PN YN GAMN ;
  269.  
  270.  RRESID RFLOT = 'KONV' 'VF' 'PERFMULT' 'RESI' METO LINCO
  271.                 $DOMTOT MRN MUN MPN MGAMN MYN PGAS ;                   
  272.  
  273. ******************************************************************
  274. ******************************************************************             
  275.  
  276.   RMAT = 'KONV' 'VF' 'PERFMULT' 'JACOCONS' METO  $DOMTOT PGAS LINCO RN
  277.             VIT PN YN ;
  278. ******************************************************************
  279.  
  280.  PCEN9 = 'POIN' 1 ('DOMA' $DOM9 CENTRE) ;
  281.  
  282. *****************************************************************
  283. ****             Beggining of the procedure            **********
  284. ***************************************************************** 
  285.  
  286.  'DEBPROC' JAEXAC RMAT*'MATRIK' $DOMA*'MMODEL' LISTINCO*'LISTMOTS'
  287.     PPRIM*'POINT' PDUAL*'POINT' MOTPRI*'MOT' MOTDUA*'MOT' ;
  288.  
  289.    MCEN = 'DOMA' $DOMA CENTRE ; 
  290.    ELT1 = 'MANUEL' 'POI1' PPRIM ;
  291.    NDIM = 'DIME'  LISTINCO ;
  292.    CHPUN = 'MANUEL' 'CHPO' ELT1 1 MOTPRI 1.0 'NATURE' 'DISCRET' ;
  293. ********************** 
  294.   'REPETER'  BL1 NDIM ;
  295.     MOTCEL = 'EXTRAIRE' LISTINCO &BL1 ;
  296.     CHPUN = CHPUN 'ET' ('MANUEL' 'CHPO' MCEN 1 MOTCEL
  297.      0.0 'NATURE' 'DISCRET') ; 
  298.   'FIN' BL1 ;
  299. *********************
  300.    CHCH = 'KOPS'  RMAT 'MULT' CHPUN ;
  301.    SCAL = 'EXTRAIRE' CHCH PDUAL MOTDUA ;
  302.  'FINPROC' SCAL ;
  303.  
  304. ****************************************************************
  305. *****    End of the procedure  JAEXAC                   ********
  306. ****************************************************************
  307.  
  308. ****************************************************************
  309. ****************   Procedure JACNUM     ************************ 
  310. ****************************************************************
  311.  'DEBPROC' JACNUM $DOMA*'MMODEL' NOMMET*'MOT' RN*'CHPOINT' 
  312.     GN*'CHPOINT' RETN*'CHPOINT' RYN*'CHPOINT' GAMN*'CHPOINT'
  313.     PGAS*'TABLE ' LISTINCO*'LISTMOTS' 
  314.     PPRIM*'POINT' PDUAL*'POINT' 
  315.     MOTPRI*'MOT' MOTDUA*'MOT' EPSILON*'FLOTTANT';
  316.  
  317.  
  318.     RNCEL   = 'COPIER'  RN ;
  319.     GNCEL   = 'COPIER'  GN ;
  320.     RETNCEL = 'COPIER'  RETN ;
  321.     RYNCEL  = 'COPIER'  RYN ;   
  322.  
  323.     VIT PRES TTEMP FRAM GAMM  = 'PRIM' 'PERFMULT'
  324.                         PGAS RNCEL GNCEL RETNCEL RYNCEL  ;
  325.  
  326.     MRN MUN MPN MYN MGAMN = 'PRET' 'PERFMULT' 1 1 $DOMTOT
  327.                     RNCEL VIT PRES FRAM GAMN ;
  328.  
  329.     CHPRESO RFLOT = 'KONV' 'VF' 'PERFMULT' 'RESI' NOMMET
  330.             LISTINCO $DOMA MRN MUN MPN MGAMN MYN PGAS ;
  331.  
  332.     VAL0 = 'EXTRAIRE' CHPRESO PDUAL MOTDUA ;
  333. ******************************************
  334.     dens0 = 'EXTRAIRE' RN PPRIM 'SCAL' ;
  335.  
  336.     CN2 = GAMN '*' (PRES '/' RN) ;
  337.  
  338.     cson0 = ('EXTRAIRE' CN2 PPRIM 'SCAL') '**' 0.5 ;
  339.  
  340.     ret0 = ('EXTRAIRE' RETN PPRIM 'SCAL') '**' 0.5 ;
  341. ******************************************    
  342.     NDIM = 'DIME' LISTINCO ;
  343.     'REPETER' BL1 NDIM ;
  344.        MOTCEL = 'EXTRAIRE' LISTINCO &BL1 ;
  345.        'SI' ('EGA' MOTCEL MOTPRI) ;
  346.           ICEL = &BL1 ;
  347.           'QUITTER' BL1 ;
  348.        'FINSI' ;
  349.     'FIN' BL1 ;
  350.  
  351.     'SI' (ICEL > NDIM) ;
  352.        'MESSAGE' 'Procedure JACNUM' ;
  353.        'MESSAGE' 'MOTPRI = ??? ';
  354.        'ERREUR' 21 ;
  355.     'FINSI' ;
  356.  
  357.     ELT1 = 'MANUEL' 'POI1' PPRIM ;
  358. ****************************************************    
  359. *   ICEL = 1 -> On perturbe la densité
  360. ****************************************************
  361.     'SI' ('EGA' ICEL 1) ;
  362.        DELTATOT = (EPSILON * dens0) ;
  363.        RNCEL = ('MANUEL' 'CHPO' ELT1 1 'SCAL' DELTATOT
  364.           'NATURE' 'DISCRET') 'ET' RN ;
  365.        GNCEL = 'COPIER' GN ;
  366.        RETNCEL = 'COPIER' RETN ;
  367.        RYNCEL = 'COPIER' RYN ;
  368.     'FINSI' ;
  369. ****************************************************
  370. *   ICEL = 2 -> On perturbe la q.d.m. long  l'ax x
  371. ****************************************************
  372.     'SI' ('EGA' ICEL 2) ;
  373.        DELTATOT = (EPSILON * dens0 * cson0) ;
  374.        GNCEL = ('MANUEL' 'CHPO' ELT1 1 'UX' DELTATOT
  375.           'NATURE' 'DISCRET') 'ET' GN ;
  376.        RNCEL = 'COPIER' RN ;
  377.        RETNCEL = 'COPIER' RETN ;
  378.        RYNCEL = 'COPIER' RYN ;
  379.     'FINSI' ;
  380. ****************************************************
  381. *   ICEL = 3 -> On perturbe la q.d.m. long  l'ax y
  382. ****************************************************
  383.     'SI' ('EGA' ICEL 3) ;
  384.        DELTATOT = (EPSILON * dens0 * cson0) ;
  385.        GNCEL = ('MANUEL' 'CHPO' ELT1 1 'UY' DELTATOT
  386.           'NATURE' 'DISCRET') 'ET' GN ;
  387.        RNCEL = 'COPIER' RN ;
  388.        RETNCEL = 'COPIER' RETN ;
  389.        RYNCEL = 'COPIER' RYN ; 
  390.     'FINSI' ;
  391. ****************************************************
  392. *   ICEL = 4 -> On perturbe l'énergie totale
  393. ****************************************************
  394.     'SI' ('EGA' ICEL 4) ;
  395.        DELTATOT = (EPSILON * ret0) ;
  396.        RETNCEL = ('MANUEL' 'CHPO' ELT1 1 'SCAL' DELTATOT
  397.           'NATURE' 'DISCRET') 'ET' RETN ;
  398.        RNCEL = 'COPIER' RN ;
  399.        GNCEL = 'COPIER' GN ;
  400.        RYNCEL = 'COPIER' RYN ;
  401.     'FINSI' ;
  402. ****************************************************
  403. *   ICEL = 5 -> On perturbe (\rho Y_{H2})
  404. ****************************************************
  405.     'SI' ('EGA' ICEL 5) ;
  406.        DELTATOT = (EPSILON * ret0) ;
  407.        RYNCEL = ('MANUEL' 'CHPO' ELT1 1 'H2' DELTATOT
  408.           'NATURE' 'DISCRET') 'ET' RYN ;
  409.        RNCEL = 'COPIER' RN ;
  410.        GNCEL = 'COPIER' GN ;
  411.        RETNCEL = 'COPIER' RETN ;
  412.     'FINSI' ;
  413. ****************************************************
  414. *   ICEL = 6 -> On perturbe (\rho Y_{O2})
  415. ****************************************************
  416.     'SI' ('EGA' ICEL 6) ;
  417.        DELTATOT = (EPSILON * ret0) ;
  418.        RYNCEL = ('MANUEL' 'CHPO' ELT1 1 'O2' DELTATOT
  419.           'NATURE' 'DISCRET') 'ET' RYN ;
  420.        RNCEL = 'COPIER' RN ;
  421.        GNCEL = 'COPIER' GN ;
  422.        RETNCEL = 'COPIER' RETN ;
  423.     'FINSI' ;
  424. ****************************************************
  425. *   ICEL = 7 -> On perturbe (\rho Y_{HO})
  426. ****************************************************
  427.     'SI' ('EGA' ICEL 7) ;
  428.        DELTATOT = (EPSILON * ret0) ;
  429.        RYNCEL = ('MANUEL' 'CHPO' ELT1 1 'H2O' DELTATOT
  430.           'NATURE' 'DISCRET') 'ET' RYN ;
  431.        RNCEL = 'COPIER' RN ;
  432.        GNCEL = 'COPIER' GN ;
  433.        RETNCEL = 'COPIER' RETN ;
  434.     'FINSI' ;            
  435. ****************************************************
  436.     VIT PRES TTEMP FRAM GAMM  = 'PRIM' 'PERFMULT'
  437.                         PGAS RNCEL GNCEL RETNCEL RYNCEL  ;
  438.  
  439.     MRNF MUNF MPNF MYNF MGAMNF = 'PRET' 'PERFMULT' 1 1 $DOMTOT
  440.                     RNCEL VIT PRES FRAM GAMM ;
  441.  
  442.     CHPRES1 RFLOT = 'KONV' 'VF' 'PERFMULT' 'RESI' NOMMET
  443.             LISTINCO $DOMA MRNF MUNF MPNF MGAMNF MYNF PGAS ;
  444.  
  445.     VAL1 = 'EXTRAIRE' CHPRES1 PDUAL MOTDUA ;                                  
  446.  
  447.  'FINPROC' ((VAL1 '-' VAL0) '/' DELTATOT) ;
  448. ****************************************************************
  449. ****************************************************************
  450.   DRR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMDEN ;
  451.   DGXR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMMOX ;
  452.   DGYR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMMOY ;
  453.   DRETR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMRET ;
  454.   DYH2R = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMYH2 ;
  455.   DYO2R = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMYO2 ;
  456.   DYHOR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMYHO ;
  457. ***************************************
  458.   DRGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMDEN ;
  459.   DGXGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMMOX ;
  460.   DGYGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMMOY ;
  461.   DRETGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMRET ;
  462.   DYH2GX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMYH2 ;
  463.   DYO2GX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMYO2 ;
  464.   DYHOGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMYHO ;
  465. ***************************************
  466.   DRGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMDEN ;
  467.   DGXGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMMOX ;
  468.   DGYGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMMOY ;
  469.   DRETGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMRET ;
  470.   DYH2GY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMYH2 ;
  471.   DYO2GY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMYO2 ;
  472.   DYHOGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMYHO ;  
  473. ***************************************
  474.   DRRET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMDEN ;
  475.   DGXRET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMMOX ;
  476.   DGYRET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMMOY ;
  477.   DRETRET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMRET ;
  478.   DYH2RET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMYH2 ;
  479.   DYO2RET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMYO2 ;
  480.   DYHORET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMYHO ;
  481. ***************************************
  482.   DRYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMDEN ;
  483.   DGXYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMMOX ;
  484.   DGYYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMMOY ;
  485.   DRETYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMRET ;
  486.   DYH2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMYH2 ;
  487.   DYO2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMYO2 ;
  488.   DYHOYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMYHO ;
  489. ***************************************
  490.   DRYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMDEN ;
  491.   DGXYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMMOX ;
  492.   DGYYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMMOY ;
  493.   DRETYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMRET ;
  494.   DYH2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMYH2 ;
  495.   DYO2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMYO2 ;
  496.   DYHOYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMYHO ;
  497. ****************************************
  498.   DRYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMDEN ;
  499.   DGXYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMMOX ;
  500.   DGYYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMMOY ;
  501.   DRETYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMRET ;
  502.   DYH2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMYH2 ;
  503.   DYO2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMYO2 ;
  504.   DYHOYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMYHO ;      
  505. ****************************************************************
  506. *****************************************************************  
  507.   DRRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  508.       PCEN9 NOMDEN NOMDEN DELTA ;  
  509.   DGXRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  510.       PCEN9 NOMDEN NOMMOX DELTA ;      
  511.   DGYRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  512.       PCEN9 NOMDEN NOMMOY DELTA ;
  513.   DRETRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  514.       PCEN9 NOMDEN NOMRET DELTA ;
  515.   DYH2RN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  516.       PCEN9 NOMDEN NOMYH2 DELTA ;
  517.   DYO2RN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  518.       PCEN9 NOMDEN NOMYO2 DELTA ;
  519.   DYHORN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  520.       PCEN9 NOMDEN NOMYHO DELTA ;
  521. **************************************
  522.   DRGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  523.       PCEN9 NOMMOX NOMDEN DELTA ;  
  524.   DGXGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  525.       PCEN9 NOMMOX NOMMOX DELTA ;      
  526.   DGYGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  527.       PCEN9 NOMMOX NOMMOY DELTA ;
  528.   DRETGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  529.       PCEN9 NOMMOX NOMRET DELTA ;
  530.   DYH2GXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  531.       PCEN9 NOMMOX NOMYH2 DELTA ;
  532.   DYO2GXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  533.       PCEN9 NOMMOX NOMYO2 DELTA ;
  534.   DYHOGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  535.       PCEN9 NOMMOX NOMYHO DELTA ;
  536. **************************************
  537.   DRGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  538.       PCEN9 NOMMOY NOMDEN DELTA ;  
  539.   DGXGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  540.       PCEN9 NOMMOY NOMMOX DELTA ;      
  541.   DGYGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  542.       PCEN9 NOMMOY NOMMOY DELTA ;
  543.   DRETGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  544.       PCEN9 NOMMOY NOMRET DELTA ;
  545.   DYH2GYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  546.       PCEN9 NOMMOY NOMYH2 DELTA ;
  547.   DYO2GYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  548.       PCEN9 NOMMOY NOMYO2 DELTA ;
  549.   DYHOGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  550.       PCEN9 NOMMOY NOMYHO DELTA ;
  551. **************************************
  552.   DRRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  553.       PCEN9 NOMRET NOMDEN DELTA ;  
  554.   DGXRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  555.       PCEN9 NOMRET NOMMOX DELTA ;      
  556.   DGYRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  557.       PCEN9 NOMRET NOMMOY DELTA ;
  558.   DRETRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  559.       PCEN9 NOMRET NOMRET DELTA ;
  560.   DYH2RETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  561.       PCEN9 NOMRET NOMYH2 DELTA ;
  562.   DYO2RETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  563.       PCEN9 NOMRET NOMYO2 DELTA ;
  564.   DYHORETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  565.       PCEN9 NOMRET NOMYHO DELTA ;
  566. **************************************
  567.   DRYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  568.       PCEN9 NOMYH2 NOMDEN DELTA ;  
  569.   DGXYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  570.       PCEN9 NOMYH2 NOMMOX DELTA ;      
  571.   DGYYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  572.       PCEN9 NOMYH2 NOMMOY DELTA ;
  573.   DRETYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  574.       PCEN9 NOMYH2 NOMRET DELTA ;
  575.   DYH2YH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  576.       PCEN9 NOMYH2 NOMYH2 DELTA ;
  577.   DYO2YH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  578.       PCEN9 NOMYH2 NOMYO2 DELTA ;
  579.   DYHOYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  580.       PCEN9 NOMYH2 NOMYHO DELTA ;
  581. ***************************************
  582.   DRYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  583.       PCEN9 NOMYO2 NOMDEN DELTA ;  
  584.   DGXYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  585.       PCEN9 NOMYO2 NOMMOX DELTA ;      
  586.   DGYYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  587.       PCEN9 NOMYO2 NOMMOY DELTA ;
  588.   DRETYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  589.       PCEN9 NOMYO2 NOMRET DELTA ;
  590.   DYH2YO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  591.       PCEN9 NOMYO2 NOMYH2 DELTA ;
  592.   DYO2YO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  593.       PCEN9 NOMYO2 NOMYO2 DELTA ;
  594.   DYHOYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  595.       PCEN9 NOMYO2 NOMYHO DELTA ;
  596. ***************************************
  597.   DRYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  598.       PCEN9 NOMYHO NOMDEN DELTA ;  
  599.   DGXYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  600.       PCEN9 NOMYHO NOMMOX DELTA ;      
  601.   DGYYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  602.       PCEN9 NOMYHO NOMMOY DELTA ;
  603.   DRETYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  604.       PCEN9 NOMYHO NOMRET DELTA ;
  605.   DYH2YHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  606.       PCEN9 NOMYHO NOMYH2 DELTA ;
  607.   DYO2YHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  608.       PCEN9 NOMYHO NOMYO2 DELTA ;
  609.   DYHOYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  610.       PCEN9 NOMYHO NOMYHO DELTA ;      
  611. **************************************************************
  612. * Test des comparaisons jacobien exact-jacobien numerique
  613. **************************************************************
  614.  'SI' (('ABS'(DRR)) < DELTA) ;
  615.     ERR1 = ('ABS' (DRRN '-' DRR)) ;
  616.  'SINON' ;  
  617.     ERR1 = ('ABS' (DRRN '-' DRR)) '/' ('ABS' (DRR)) ;
  618.  'FINSI' ; 
  619. ***********************
  620.  'SI' (ERR1 > ERRTOL) ;
  621.       'MESSAGE' 'Problem 1a';
  622.       'ERREUR' 5 ;
  623.  'FINSI' ;
  624. *************************************************************
  625.  'SI' (('ABS'(DGXR)) < DELTA) ;
  626.     ERR1 = ('ABS' (DGXRN '-' DGXR)) ;
  627.  'SINON' ;  
  628.     ERR1 = ('ABS' (DGXRN '-' DGXR)) '/' ('ABS' (DGXR)) ;
  629.  'FINSI' ; 
  630. ***********************
  631.  'SI' (ERR1 > ERRTOL) ;
  632.       'MESSAGE' 'Problem 2a';
  633.       'ERREUR' 5 ;
  634.  'FINSI' ;      
  635. **************************************************************
  636.  'SI' (('ABS'(DGYR)) < DELTA) ;
  637.     ERR1 = ('ABS' (DGYRN '-' DGYR)) ;
  638.  'SINON' ;  
  639.     ERR1 = ('ABS' (DGYRN '-' DGYR)) '/' ('ABS' (DGYR)) ;
  640.  'FINSI' ; 
  641. ************************
  642.  'SI' (ERR1 > ERRTOL) ;
  643.       'MESSAGE' 'Problem 3a';
  644.        'ERREUR' 5 ;
  645.  'FINSI' ;
  646. **************************************************************
  647.  'SI' (('ABS'(DRETR)) < DELTA) ;
  648.     ERR1 = ('ABS' (DRETRN '-' DRETR)) ;
  649.  'SINON' ;  
  650.     ERR1 = ('ABS' (DRETRN '-' DRETR)) '/' ('ABS' (DRETR)) ;
  651.  'FINSI' ;
  652. **************************
  653.  'SI' (ERR1 > ERRTOL) ;
  654.       'MESSAGE' 'Problem 4a';
  655.       'ERREUR' 5 ;
  656.  'FINSI' ;
  657. **************************************************************
  658.  'SI' (('ABS'(DYH2R)) < DELTA) ;
  659.     ERR1 = ('ABS' (DYH2RN '-' DYH2R)) ;
  660.  'SINON' ;  
  661.     ERR1 = ('ABS' (DYH2RN '-' DYH2R)) '/' ('ABS' (DYH2R)) ;
  662.  'FINSI' ;
  663. **************************
  664.  'SI' (ERR1 > ERRTOL) ;
  665.       'MESSAGE' 'Problem 5a';
  666.       'ERREUR' 5 ;
  667.  'FINSI' ;
  668. **************************************************************
  669.  'SI' (('ABS'(DYO2R)) < DELTA) ;
  670.     ERR1 = ('ABS' (DYO2RN '-' DYO2R)) ;
  671.  'SINON' ;  
  672.     ERR1 = ('ABS' (DYO2RN '-' DYO2R)) '/' ('ABS' (DYO2R)) ;
  673.  'FINSI' ;
  674. **************************
  675.  'SI' (ERR1 > ERRTOL) ;
  676.       'MESSAGE' 'Problem 6a';
  677.       'ERREUR' 5 ;
  678.  'FINSI' ;
  679. **************************************************************
  680.  'SI' (('ABS'(DYHOR)) < DELTA) ;
  681.     ERR1 = ('ABS' (DYHORN '-' DYHOR)) ;
  682.  'SINON' ;  
  683.     ERR1 = ('ABS' (DYHORN '-' DYHOR)) '/' ('ABS' (DYHOR)) ;
  684.  'FINSI' ;
  685. **************************
  686.  'SI' (ERR1 > ERRTOL) ;
  687.       'MESSAGE' 'Problem 7a';
  688.       'ERREUR' 5 ;
  689.  'FINSI' ;
  690. **************************************************************
  691. *****    Second line !!!!!!!!!!!!!!!!!!!!!!!    **************
  692. **************************************************************
  693. **************************************************************
  694.  'SI' (('ABS'(DRGX)) < DELTA) ;
  695.     ERR1 = ('ABS' (DRGXN '-' DRGX)) ;
  696.  'SINON' ;   
  697.     ERR1 = ('ABS' (DRGXN '-' DRGX)) '/' ('ABS' (DRGX))  ;
  698.  'FINSI' ;
  699. ***************************
  700.  'SI' (ERR1 > ERRTOL) ;
  701.       'MESSAGE' 'Problem 8a';
  702.       'ERREUR' 5 ;
  703.  'FINSI' ;
  704. *************************************************************
  705.  'SI' (('ABS'(DGXGX)) < DELTA) ;
  706.    ERR1 = ('ABS' (DGXGXN '-' DGXGX)) ;
  707.  'SINON' ;
  708.    ERR1 = ('ABS' (DGXGXN '-' DGXGX)) '/' ('ABS' (DGXGX)) ;
  709.  'FINSI' ;  
  710. *****************************   
  711.  'SI' (ERR1 > ERRTOL) ;
  712.         'MESSAGE' 'Problem 9a';
  713.         'ERREUR' 5 ;
  714.  'FINSI' ;
  715. *************************************************************
  716.  'SI' (('ABS'(DGYGX)) < DELTA) ;
  717.     ERR1 = ('ABS' (DGYGXN '-' DGYGX)) ;
  718.  'SINON' ;
  719.     ERR1 = ('ABS' (DGYGXN '-' DGYGX)) '/' ('ABS' (DGYGX)) ;
  720.  'FINSI' ;
  721. ******************************
  722.  'SI' (ERR1 > ERRTOL) ;
  723.         'MESSAGE' 'Problem 10a';
  724.         'ERREUR' 5 ;
  725.  'FINSI' ;
  726. *************************************************************
  727.  'SI' (('ABS'(DRETGX)) < DELTA) ;
  728.     ERR1 = ('ABS' (DRETGXN '-' DRETGX)) ;
  729.  'SINON' ;
  730.     ERR1 = ('ABS' (DRETGXN '-' DRETGX)) '/' ('ABS' (DRETGX)) ;
  731.  'FINSI' ;
  732. *******************************    
  733.  'SI' (ERR1 > ERRTOL) ;
  734.         'MESSAGE' 'Problem 11a';
  735.         'ERREUR' 5 ;
  736.  'FINSI' ;
  737. **************************************************************
  738.  'SI' (('ABS'(DYH2GX)) < DELTA) ;
  739.     ERR1 = ('ABS' (DYH2GXN '-' DYH2GX)) ;
  740.  'SINON' ;  
  741.     ERR1 = ('ABS' (DYH2GXN '-' DYH2GX)) '/' ('ABS' (DYH2GX)) ;
  742.  'FINSI' ;
  743. **************************
  744.  'SI' (ERR1 > ERRTOL) ;
  745.       'MESSAGE' 'Problem 12a';
  746.       'ERREUR' 5 ;
  747.  'FINSI' ;
  748. **************************************************************
  749.  'SI' (('ABS'(DYO2GX)) < DELTA) ;
  750.     ERR1 = ('ABS' (DYO2GXN '-' DYO2GX)) ;
  751.  'SINON' ;  
  752.     ERR1 = ('ABS' (DYO2GXN '-' DYO2GX)) '/' ('ABS' (DYO2GX)) ;
  753.  'FINSI' ;
  754. **************************
  755.  'SI' (ERR1 > ERRTOL) ;
  756.       'MESSAGE' 'Problem 13a';
  757.       'ERREUR' 5 ;
  758.  'FINSI' ;
  759. **************************************************************
  760.  'SI' (('ABS'(DYHOGX)) < DELTA) ;
  761.     ERR1 = ('ABS' (DYHOGXN '-' DYHOGX)) ;
  762.  'SINON' ;  
  763.     ERR1 = ('ABS' (DYHOGXN '-' DYHOGX)) '/' ('ABS' (DYHOGX)) ;
  764.  'FINSI' ;
  765. **************************
  766.  'SI' (ERR1 > ERRTOL) ;
  767.       'MESSAGE' 'Problem 14a';
  768.       'ERREUR' 5 ;
  769.  'FINSI' ;
  770. **************************************************************
  771. *****    Third line !!!!!!!!!!!!!!!!!!!!!!!    **************
  772. **************************************************************
  773. **************************************************************
  774.  'SI' (('ABS'(DRGY)) < DELTA) ;
  775.     ERR1 = ('ABS' (DRGYN '-' DRGY)) ;
  776.  'SINON' ;
  777.     ERR1 = ('ABS' (DRGYN '-' DRGY)) '/' ('ABS' (DRGY)) ;
  778.  'FINSI' ;
  779. *******************************
  780.  'SI' (ERR1 > ERRTOL) ;
  781.         'MESSAGE' 'Problem 15a';
  782.         'ERREUR' 5 ;
  783.  'FINSI' ;
  784. **************************************************************
  785.  'SI' (('ABS'(DGXGY)) < DELTA) ;
  786.     ERR1 = ('ABS' (DGXGYN '-' DGXGY)) ;
  787.  'SINON' ;
  788.     ERR1 = ('ABS' (DGXGYN '-' DGXGY))  '/' ('ABS' (DGXGY)) ;
  789.  'FINSI' ;
  790. ********************************
  791.  'SI' (ERR1 > ERRTOL) ;
  792.         'MESSAGE' 'Problem 16a';
  793.         'ERREUR' 5 ;
  794.  'FINSI' ;
  795. **************************************************************
  796.  'SI' (('ABS'(DGYGY)) < DELTA) ;
  797.     ERR1 = ('ABS' (DGYGYN '-' DGYGY)) ;
  798.  'SINON' ;
  799.     ERR1 = ('ABS' (DGYGYN '-' DGYGY)) '/' ('ABS' (DGYGY)) ;
  800.  'FINSI' ;
  801. ***************
  802.  'SI' (ERR1 > ERRTOL) ;
  803.         'MESSAGE' 'Problem 17a';
  804.         'ERREUR' 5 ;
  805.  'FINSI' ;
  806. *************************************************************
  807.  'SI' (('ABS'(DRETGY)) < DELTA) ;
  808.     ERR1 = ('ABS' (DRETGYN '-' DRETGY)) ;
  809.  'SINON' ;
  810.     ERR1 = ('ABS' (DRETGYN '-' DRETGY)) '/' ('ABS' (DRETGY)) ;
  811.  'FINSI' ;
  812. ****************     
  813.  'SI' (ERR1 > ERRTOL) ;
  814.         'MESSAGE' 'Problem 18a';
  815.         'ERREUR' 5 ;
  816.  'FINSI' ;
  817. **************************************************************
  818.  'SI' (('ABS'(DYH2GY)) < DELTA) ;
  819.     ERR1 = ('ABS' (DYH2GYN '-' DYH2GY)) ;
  820.  'SINON' ;  
  821.     ERR1 = ('ABS' (DYH2GYN '-' DYH2GY)) '/' ('ABS' (DYH2GY)) ;
  822.  'FINSI' ;
  823. **************************
  824.  'SI' (ERR1 > ERRTOL) ;
  825.       'MESSAGE' 'Problem 19a';
  826.       'ERREUR' 5 ;
  827.  'FINSI' ;
  828. **************************************************************
  829.  'SI' (('ABS'(DYO2GY)) < DELTA) ;
  830.     ERR1 = ('ABS' (DYO2GYN '-' DYO2GY)) ;
  831.  'SINON' ;  
  832.     ERR1 = ('ABS' (DYO2GYN '-' DYO2GY)) '/' ('ABS' (DYO2GY)) ;
  833.  'FINSI' ;
  834. **************************
  835.  'SI' (ERR1 > ERRTOL) ;
  836.       'MESSAGE' 'Problem 20a';
  837.       'ERREUR' 5 ;
  838.  'FINSI' ;
  839. **************************************************************
  840.  'SI' (('ABS'(DYHOGY)) < DELTA) ;
  841.     ERR1 = ('ABS' (DYHOGYN '-' DYHOGY)) ;
  842.  'SINON' ;  
  843.     ERR1 = ('ABS' (DYHOGYN '-' DYHOGY)) '/' ('ABS' (DYHOGY)) ;
  844.  'FINSI' ;
  845. **************************
  846.  'SI' (ERR1 > ERRTOL) ;
  847.       'MESSAGE' 'Problem 21a';
  848.       'ERREUR' 5 ;
  849.  'FINSI' ;
  850. **************************************************************
  851. *******     Fouth line!!!!!!!!!!!!!!!!!!!   ******************
  852. **************************************************************
  853. **************************************************************
  854.  'SI' (('ABS'(DRRET)) < DELTA) ;
  855.     ERR1 = ('ABS' (DRRETN '-' DRRET)) ;
  856.  'SINON' ;
  857.     ERR1 = ('ABS' (DRRETN '-' DRRET)) '/' ('ABS' (DRRET)) ;
  858.  'FINSI' ;
  859. *****************
  860.  'SI' (ERR1 > ERRTOL) ;
  861.         'MESSAGE' 'Problem 22a';
  862.         'ERREUR' 5 ;
  863.  'FINSI' ;
  864. **************************************************************
  865.  'SI' (('ABS'(DGXRET)) < DELTA) ;
  866.     ERR1 = ('ABS' (DGXRETN '-' DGXRET)) ;
  867.  'SINON' ;
  868.     ERR1 = ('ABS' (DGXRETN '-' DGXRET)) '/' ('ABS' (DGXRET)) ;
  869.  'FINSI' ;
  870. *****************
  871.  'SI' (ERR1 > ERRTOL) ;
  872.         'MESSAGE' 'Problem 23a';
  873.         'ERREUR' 5 ;
  874.  'FINSI' ;
  875. *************************************************************
  876.  'SI' (('ABS'(DGYRET)) < DELTA) ;
  877.     ERR1 = ('ABS' (DGYRETN '-' DGYRET)) ; 
  878.  'SINON' ;
  879.     ERR1 = ('ABS' (DGYRETN '-' DGYRET)) '/' ('ABS' (DGYRET)) ;
  880.  'FINSI' ;
  881. *****************    
  882.  'SI' (ERR1 > ERRTOL) ;
  883.         'MESSAGE' 'Problem 24a';
  884.         'ERREUR' 5 ;
  885.  'FINSI' ;
  886. **************************************************************
  887.  'SI' (('ABS'(DRETRET)) < DELTA) ;
  888.     ERR1 = ('ABS' (DRETRETN '-' DRETRET)) ;
  889.  'SINON' ;
  890.     ERR1 = ('ABS' (DRETRETN '-' DRETRET)) '/' ('ABS' (DRETRET)) ;
  891.  'FINSI' ;
  892. *****************
  893.  'SI' (ERR1 > ERRTOL) ;
  894.         'MESSAGE' 'Problem 25a';
  895.         'ERREUR' 5 ;
  896.  'FINSI' ;
  897. **************************************************************
  898.  'SI' (('ABS'(DYH2RET)) < DELTA) ;
  899.     ERR1 = ('ABS' (DYH2RETN '-' DYH2RET)) ;
  900.  'SINON' ;  
  901.     ERR1 = ('ABS' (DYH2RETN '-' DYH2RET)) '/' ('ABS' (DYH2RET)) ;
  902.  'FINSI' ;
  903. **************************
  904.  'SI' (ERR1 > ERRTOL) ;
  905.       'MESSAGE' 'Problem 26a';
  906.       'ERREUR' 5 ;
  907.  'FINSI' ;
  908. **************************************************************
  909.  'SI' (('ABS'(DYO2RET)) < DELTA) ;
  910.     ERR1 = ('ABS' (DYO2RETN '-' DYO2RET)) ;
  911.  'SINON' ;  
  912.     ERR1 = ('ABS' (DYO2RETN '-' DYO2RET)) '/' ('ABS' (DYO2RET)) ;
  913.  'FINSI' ;
  914. **************************
  915.  'SI' (ERR1 > ERRTOL) ;
  916.       'MESSAGE' 'Problem 27a';
  917.       'ERREUR' 5 ;
  918.  'FINSI' ;
  919. **************************************************************
  920.  'SI' (('ABS'(DYHORET)) < DELTA) ;
  921.     ERR1 = ('ABS' (DYHORETN '-' DYHORET)) ;
  922.  'SINON' ;  
  923.     ERR1 = ('ABS' (DYHORETN '-' DYHORET)) '/' ('ABS' (DYHORET)) ;
  924.  'FINSI' ;
  925. **************************
  926.  'SI' (ERR1 > ERRTOL) ;
  927.       'MESSAGE' 'Problem 28a';
  928.       'ERREUR' 5 ;
  929.  'FINSI' ;
  930. **************************************************************
  931. *******     Fith line!!!!!!!!!!!!!!!!!!!   ******************
  932. **************************************************************
  933.  'SI' (('ABS'(DRYH2)) < DELTA) ;
  934.     ERR1 = ('ABS' (DRYH2N '-' DRYH2)) ;
  935.  'SINON' ;
  936.     ERR1 = ('ABS' (DRYH2N '-' DRYH2)) '/' ('ABS' (DRYH2)) ;
  937.  'FINSI' ;
  938. *****************
  939.  'SI' (ERR1 > ERRTOL) ;
  940.         'MESSAGE' 'Problem 29a';
  941.         'ERREUR' 5 ;
  942.  'FINSI' ;
  943. **************************************************************
  944.  'SI' (('ABS'(DGXYH2)) < DELTA) ;
  945.     ERR1 = ('ABS' (DGXYH2N '-' DGXYH2)) ;
  946.  'SINON' ;
  947.     ERR1 = ('ABS' (DGXYH2N '-' DGXYH2)) '/' ('ABS' (DGXYH2)) ;
  948.  'FINSI' ;
  949. *****************
  950.  'SI' (ERR1 > ERRTOL) ;
  951.         'MESSAGE' 'Problem 30a';
  952.         'ERREUR' 5 ;
  953.  'FINSI' ;
  954. *************************************************************
  955.  'SI' (('ABS'(DGYYH2)) < DELTA) ;
  956.     ERR1 = ('ABS' (DGYYH2N '-' DGYYH2)) ; 
  957.  'SINON' ;
  958.     ERR1 = ('ABS' (DGYYH2N '-' DGYYH2)) '/' ('ABS' (DGYYH2)) ;
  959.  'FINSI' ;
  960. *****************    
  961.  'SI' (ERR1 > ERRTOL) ;
  962.         'MESSAGE' 'Problem 31a';
  963.         'ERREUR' 5 ;
  964.  'FINSI' ;
  965. **************************************************************
  966.  'SI' (('ABS'(DRETYH2)) < DELTA) ;
  967.     ERR1 = ('ABS' (DRETYH2N '-' DRETYH2)) ;
  968.  'SINON' ;
  969.     ERR1 = ('ABS' (DRETYH2N '-' DRETYH2)) '/' ('ABS' (DRETYH2)) ;
  970.  'FINSI' ;
  971. *****************
  972.  'SI' (ERR1 > ERRTOL) ;
  973.         'MESSAGE' 'Problem 32a';
  974.         'ERREUR' 5 ;
  975.  'FINSI' ;
  976. **************************************************************
  977.  'SI' (('ABS'(DYH2YH2)) < DELTA) ;
  978.     ERR1 = ('ABS' (DYH2YH2N '-' DYH2YH2)) ;
  979.  'SINON' ;  
  980.     ERR1 = ('ABS' (DYH2YH2N '-' DYH2YH2)) '/' ('ABS' (DYH2YH2)) ;
  981.  'FINSI' ;
  982. **************************
  983.  'SI' (ERR1 > ERRTOL) ;
  984.       'MESSAGE' 'Problem 33a';
  985.       'ERREUR' 5 ;
  986.  'FINSI' ;
  987. **************************************************************
  988.  'SI' (('ABS'(DYO2YH2)) < DELTA) ;
  989.     ERR1 = ('ABS' (DYO2YH2N '-' DYO2YH2)) ;
  990.  'SINON' ;  
  991.     ERR1 = ('ABS' (DYO2YH2N '-' DYO2YH2)) '/' ('ABS' (DYO2YH2)) ;
  992.  'FINSI' ;
  993. **************************
  994.  'SI' (ERR1 > ERRTOL) ;
  995.       'MESSAGE' 'Problem 34a';
  996.       'ERREUR' 5 ;
  997.  'FINSI' ;
  998. **************************************************************
  999.  'SI' (('ABS'(DYHOYH2)) < DELTA) ;
  1000.     ERR1 = ('ABS' (DYHOYH2N '-' DYHOYH2)) ;
  1001.  'SINON' ;  
  1002.     ERR1 = ('ABS' (DYHOYH2N '-' DYHOYH2)) '/' ('ABS' (DYHOYH2)) ;
  1003.  'FINSI' ;
  1004. **************************
  1005.  'SI' (ERR1 > ERRTOL) ;
  1006.       'MESSAGE' 'Problem 35a';
  1007.       'ERREUR' 5 ;
  1008.  'FINSI' ;
  1009. **************************************************************
  1010. *********  Sixth line !!!!!!!!!!!!!!!!!!!!!!   ***************
  1011. **************************************************************
  1012.  'SI' (('ABS'(DRYO2)) < DELTA) ;
  1013.     ERR1 = ('ABS' (DRYO2N '-' DRYO2)) ;
  1014.  'SINON' ;
  1015.     ERR1 = ('ABS' (DRYO2N '-' DRYO2)) '/' ('ABS' (DRYO2)) ;
  1016.  'FINSI' ;
  1017. *****************
  1018.  'SI' (ERR1 > ERRTOL) ;
  1019.         'MESSAGE' 'Problem 36a';
  1020.         'ERREUR' 5 ;
  1021.  'FINSI' ;
  1022. **************************************************************
  1023.  'SI' (('ABS'(DGXYO2)) < DELTA) ;
  1024.     ERR1 = ('ABS' (DGXYO2N '-' DGXYO2)) ;
  1025.  'SINON' ;
  1026.     ERR1 = ('ABS' (DGXYO2N '-' DGXYO2)) '/' ('ABS' (DGXYO2)) ;
  1027.  'FINSI' ;
  1028. *****************
  1029.  'SI' (ERR1 > ERRTOL) ;
  1030.         'MESSAGE' 'Problem 37a';
  1031.         'ERREUR' 5 ;
  1032.  'FINSI' ;
  1033. *************************************************************
  1034.  'SI' (('ABS'(DGYYO2)) < DELTA) ;
  1035.     ERR1 = ('ABS' (DGYYO2N '-' DGYYO2)) ; 
  1036.  'SINON' ;
  1037.     ERR1 = ('ABS' (DGYYO2N '-' DGYYO2)) '/' ('ABS' (DGYYO2)) ;
  1038.  'FINSI' ;
  1039. *****************    
  1040.  'SI' (ERR1 > ERRTOL) ;
  1041.         'MESSAGE' 'Problem 38a';
  1042.         'ERREUR' 5 ;
  1043.  'FINSI' ;
  1044. **************************************************************
  1045.  'SI' (('ABS'(DRETYO2)) < DELTA) ;
  1046.     ERR1 = ('ABS' (DRETYO2N '-' DRETYO2)) ;
  1047.  'SINON' ;
  1048.     ERR1 = ('ABS' (DRETYO2N '-' DRETYO2)) '/' ('ABS' (DRETYO2)) ;
  1049.  'FINSI' ;
  1050. *****************
  1051.  'SI' (ERR1 > ERRTOL) ;
  1052.         'MESSAGE' 'Problem 39a';
  1053.         'ERREUR' 5 ;
  1054.  'FINSI' ;
  1055. **************************************************************
  1056.  'SI' (('ABS'(DYH2YO2)) < DELTA) ;
  1057.     ERR1 = ('ABS' (DYH2YO2N '-' DYH2YO2)) ;
  1058.  'SINON' ;  
  1059.     ERR1 = ('ABS' (DYH2YO2N '-' DYH2YO2)) '/' ('ABS' (DYH2YO2)) ;
  1060.  'FINSI' ;
  1061. **************************
  1062.  'SI' (ERR1 > ERRTOL) ;
  1063.       'MESSAGE' 'Problem 40a';
  1064.       'ERREUR' 5 ;
  1065.  'FINSI' ;
  1066. **************************************************************
  1067.  'SI' (('ABS'(DYO2YO2)) < DELTA) ;
  1068.     ERR1 = ('ABS' (DYO2YO2N '-' DYO2YO2)) ;
  1069.  'SINON' ;  
  1070.     ERR1 = ('ABS' (DYO2YO2N '-' DYO2YO2)) '/' ('ABS' (DYO2YO2)) ;
  1071.  'FINSI' ;
  1072. **************************
  1073.  'SI' (ERR1 > ERRTOL) ;
  1074.       'MESSAGE' 'Problem 41a';
  1075.       'ERREUR' 5 ;
  1076.  'FINSI' ;
  1077. **************************************************************
  1078.  'SI' (('ABS'(DYHOYO2)) < DELTA) ;
  1079.     ERR1 = ('ABS' (DYHOYO2N '-' DYHOYO2)) ;
  1080.  'SINON' ;  
  1081.     ERR1 = ('ABS' (DYHOYO2N '-' DYHOYO2)) '/' ('ABS' (DYHOYO2)) ;
  1082.  'FINSI' ;
  1083. **************************
  1084.  'SI' (ERR1 > ERRTOL) ;
  1085.       'MESSAGE' 'Problem 42a';
  1086.       'ERREUR' 5 ;
  1087.  'FINSI' ;
  1088. **************************************************************
  1089. *********  Seventh line !!!!!!!!!!!!!!!!!!!!!!   *************
  1090. **************************************************************
  1091.  'SI' (('ABS'(DRYHO)) < DELTA) ;
  1092.     ERR1 = ('ABS' (DRYHON '-' DRYHO)) ;
  1093.  'SINON' ;
  1094.     ERR1 = ('ABS' (DRYHON '-' DRYHO)) '/' ('ABS' (DRYHO)) ;
  1095.  'FINSI' ;
  1096. *****************
  1097.  'SI' (ERR1 > ERRTOL) ;
  1098.         'MESSAGE' 'Problem 43a';
  1099.         'ERREUR' 5 ;
  1100.  'FINSI' ;
  1101. **************************************************************
  1102.  'SI' (('ABS'(DGXYHO)) < DELTA) ;
  1103.     ERR1 = ('ABS' (DGXYHON '-' DGXYHO)) ;
  1104.  'SINON' ;
  1105.     ERR1 = ('ABS' (DGXYHON '-' DGXYHO)) '/' ('ABS' (DGXYHO)) ;
  1106.  'FINSI' ;
  1107. *****************
  1108.  'SI' (ERR1 > ERRTOL) ;
  1109.         'MESSAGE' 'Problem 44a';
  1110.         'ERREUR' 5 ;
  1111.  'FINSI' ;
  1112. *************************************************************
  1113.  'SI' (('ABS'(DGYYHO)) < DELTA) ;
  1114.     ERR1 = ('ABS' (DGYYHON '-' DGYYHO)) ; 
  1115.  'SINON' ;
  1116.     ERR1 = ('ABS' (DGYYHON '-' DGYYHO)) '/' ('ABS' (DGYYHO)) ;
  1117.  'FINSI' ;
  1118. *****************    
  1119.  'SI' (ERR1 > ERRTOL) ;
  1120.         'MESSAGE' 'Problem 45a';
  1121.         'ERREUR' 5 ;
  1122.  'FINSI' ;
  1123. **************************************************************
  1124.  'SI' (('ABS'(DRETYHO)) < DELTA) ;
  1125.     ERR1 = ('ABS' (DRETYHON '-' DRETYHO)) ;
  1126.  'SINON' ;
  1127.     ERR1 = ('ABS' (DRETYHON '-' DRETYHO)) '/' ('ABS' (DRETYHO)) ;
  1128.  'FINSI' ;
  1129. *****************
  1130.  'SI' (ERR1 > ERRTOL) ;
  1131.         'MESSAGE' 'Problem 46a';
  1132.         'ERREUR' 5 ;
  1133.  'FINSI' ;
  1134. **************************************************************
  1135.  'SI' (('ABS'(DYH2YHO)) < DELTA) ;
  1136.     ERR1 = ('ABS' (DYH2YHON '-' DYH2YHO)) ;
  1137.  'SINON' ;  
  1138.     ERR1 = ('ABS' (DYH2YHON '-' DYH2YHO)) '/' ('ABS' (DYH2YHO)) ;
  1139.  'FINSI' ;
  1140. **************************
  1141.  'SI' (ERR1 > ERRTOL) ;
  1142.       'MESSAGE' 'Problem 47a';
  1143.       'ERREUR' 5 ;
  1144.  'FINSI' ;
  1145. **************************************************************
  1146.  'SI' (('ABS'(DYO2YHO)) < DELTA) ;
  1147.     ERR1 = ('ABS' (DYO2YHON '-' DYO2YHO)) ;
  1148.  'SINON' ;  
  1149.     ERR1 = ('ABS' (DYO2YHON '-' DYO2YHO)) '/' ('ABS' (DYO2YHO)) ;
  1150.  'FINSI' ;
  1151. **************************
  1152.  'SI' (ERR1 > ERRTOL) ;
  1153.       'MESSAGE' 'Problem 48a';
  1154.       'ERREUR' 5 ;
  1155.  'FINSI' ;
  1156. **************************************************************
  1157.  'SI' (('ABS'(DYHOYHO)) < DELTA) ;
  1158.     ERR1 = ('ABS' (DYHOYHON '-' DYHOYHO)) ;
  1159.  'SINON' ;  
  1160.     ERR1 = ('ABS' (DYHOYHON '-' DYHOYHO)) '/' ('ABS' (DYHOYHO)) ;
  1161.  'FINSI' ;
  1162. **************************
  1163.  'SI' (ERR1 > ERRTOL) ;
  1164.       'MESSAGE' 'Problem 49a';
  1165.       'ERREUR' 5 ;
  1166.  'FINSI' ;
  1167. *****************************************************
  1168. *****************************************************
  1169. *****************************************************
  1170. *****************************************************
  1171. ******* TEST2 ***************************************
  1172. *****************************************************
  1173. *****************************************************
  1174. *
  1175. * On compare le jacobien et la variation du residu
  1176. * en $DOM9 . 'CENTRE' par rapport à une variation
  1177. * infinitésimal en $DOM6 . 'CENTRE'
  1178. *
  1179.  
  1180.  PCEN6 = 'POIN' 1 ('DOMA' $DOM6 CENTRE) ;
  1181.  
  1182. ****************************************************
  1183. * Le jacobien exact. 
  1184. ****************************************************
  1185.   DRR = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMDEN NOMDEN ;
  1186.   DGXR = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMDEN NOMMOX ;
  1187.   DGYR = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMDEN NOMMOY ;
  1188.   DRETR = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMDEN NOMRET ;
  1189.   DYH2R = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMDEN NOMYH2 ;
  1190.   DYO2R = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMDEN NOMYO2 ;
  1191.   DYHOR = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMDEN NOMYHO ;
  1192. ***************************************
  1193.   DRGX = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOX NOMDEN ;
  1194.   DGXGX = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOX NOMMOX ;
  1195.   DGYGX = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOX NOMMOY ;
  1196.   DRETGX = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOX NOMRET ;
  1197.   DYH2GX = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOX NOMYH2 ;
  1198.   DYO2GX = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOX NOMYO2 ;
  1199.   DYHOGX = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOX NOMYHO ;
  1200. ***************************************
  1201.   DRGY = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOY NOMDEN ;
  1202.   DGXGY = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOY NOMMOX ;
  1203.   DGYGY = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOY NOMMOY ;
  1204.   DRETGY = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOY NOMRET ;
  1205.   DYH2GY = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOY NOMYH2 ;
  1206.   DYO2GY = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOY NOMYO2 ;
  1207.   DYHOGY = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMMOY NOMYHO ;  
  1208. ***************************************
  1209.   DRRET = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMRET NOMDEN ;
  1210.   DGXRET = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMRET NOMMOX ;
  1211.   DGYRET = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMRET NOMMOY ;
  1212.   DRETRET = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMRET NOMRET ;
  1213.   DYH2RET = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMRET NOMYH2 ;
  1214.   DYO2RET = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMRET NOMYO2 ;
  1215.   DYHORET = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMRET NOMYHO ;
  1216. ***************************************
  1217.   DRYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYH2 NOMDEN ;
  1218.   DGXYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYH2 NOMMOX ;
  1219.   DGYYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYH2 NOMMOY ;
  1220.   DRETYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYH2 NOMRET ;
  1221.   DYH2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYH2 NOMYH2 ;
  1222.   DYO2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYH2 NOMYO2 ;
  1223.   DYHOYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYH2 NOMYHO ;
  1224. ***************************************
  1225.   DRYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYO2 NOMDEN ;
  1226.   DGXYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYO2 NOMMOX ;
  1227.   DGYYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYO2 NOMMOY ;
  1228.   DRETYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYO2 NOMRET ;
  1229.   DYH2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYO2 NOMYH2 ;
  1230.   DYO2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYO2 NOMYO2 ;
  1231.   DYHOYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYO2 NOMYHO ;
  1232. ****************************************
  1233.   DRYHO = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYHO NOMDEN ;
  1234.   DGXYHO = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYHO NOMMOX ;
  1235.   DGYYHO = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYHO NOMMOY ;
  1236.   DRETYHO = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYHO NOMRET ;
  1237.   DYH2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYHO NOMYH2 ;
  1238.   DYO2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYHO NOMYO2 ;
  1239.   DYHOYHO = JAEXAC RMAT $DOMTOT LINCO PCEN6 PCEN9 NOMYHO NOMYHO ;
  1240. *****************************************************************  
  1241.   DRRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1242.       PCEN9 NOMDEN NOMDEN DELTA ;  
  1243.   DGXRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1244.       PCEN9 NOMDEN NOMMOX DELTA ;      
  1245.   DGYRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1246.       PCEN9 NOMDEN NOMMOY DELTA ;
  1247.   DRETRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1248.       PCEN9 NOMDEN NOMRET DELTA ;
  1249.   DYH2RN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1250.       PCEN9 NOMDEN NOMYH2 DELTA ;
  1251.   DYO2RN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1252.       PCEN9 NOMDEN NOMYO2 DELTA ;
  1253.   DYHORN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1254.       PCEN9 NOMDEN NOMYHO DELTA ;
  1255. **************************************
  1256.   DRGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1257.       PCEN9 NOMMOX NOMDEN DELTA ;  
  1258.   DGXGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1259.       PCEN9 NOMMOX NOMMOX DELTA ;      
  1260.   DGYGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1261.       PCEN9 NOMMOX NOMMOY DELTA ;
  1262.   DRETGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1263.       PCEN9 NOMMOX NOMRET DELTA ;
  1264.   DYH2GXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1265.       PCEN9 NOMMOX NOMYH2 DELTA ;
  1266.   DYO2GXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1267.       PCEN9 NOMMOX NOMYO2 DELTA ;
  1268.   DYHOGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1269.       PCEN9 NOMMOX NOMYHO DELTA ;
  1270. **************************************
  1271.   DRGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1272.       PCEN9 NOMMOY NOMDEN DELTA ;  
  1273.   DGXGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1274.       PCEN9 NOMMOY NOMMOX DELTA ;      
  1275.   DGYGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1276.       PCEN9 NOMMOY NOMMOY DELTA ;
  1277.   DRETGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1278.       PCEN9 NOMMOY NOMRET DELTA ;
  1279.   DYH2GYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1280.       PCEN9 NOMMOY NOMYH2 DELTA ;
  1281.   DYO2GYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1282.       PCEN9 NOMMOY NOMYO2 DELTA ;
  1283.   DYHOGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1284.       PCEN9 NOMMOY NOMYHO DELTA ;
  1285. **************************************
  1286.   DRRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1287.       PCEN9 NOMRET NOMDEN DELTA ;  
  1288.   DGXRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1289.       PCEN9 NOMRET NOMMOX DELTA ;      
  1290.   DGYRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1291.       PCEN9 NOMRET NOMMOY DELTA ;
  1292.   DRETRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1293.       PCEN9 NOMRET NOMRET DELTA ;
  1294.   DYH2RETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1295.       PCEN9 NOMRET NOMYH2 DELTA ;
  1296.   DYO2RETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1297.       PCEN9 NOMRET NOMYO2 DELTA ;
  1298.   DYHORETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1299.       PCEN9 NOMRET NOMYHO DELTA ;
  1300. **************************************
  1301.   DRYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1302.       PCEN9 NOMYH2 NOMDEN DELTA ;  
  1303.   DGXYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1304.       PCEN9 NOMYH2 NOMMOX DELTA ;      
  1305.   DGYYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1306.       PCEN9 NOMYH2 NOMMOY DELTA ;
  1307.   DRETYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1308.       PCEN9 NOMYH2 NOMRET DELTA ;
  1309.   DYH2YH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1310.       PCEN9 NOMYH2 NOMYH2 DELTA ;
  1311.   DYO2YH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1312.       PCEN9 NOMYH2 NOMYO2 DELTA ;
  1313.   DYHOYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1314.       PCEN9 NOMYH2 NOMYHO DELTA ;
  1315. ***************************************
  1316.   DRYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1317.       PCEN9 NOMYO2 NOMDEN DELTA ;  
  1318.   DGXYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1319.       PCEN9 NOMYO2 NOMMOX DELTA ;      
  1320.   DGYYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1321.       PCEN9 NOMYO2 NOMMOY DELTA ;
  1322.   DRETYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1323.       PCEN9 NOMYO2 NOMRET DELTA ;
  1324.   DYH2YO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1325.       PCEN9 NOMYO2 NOMYH2 DELTA ;
  1326.   DYO2YO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1327.       PCEN9 NOMYO2 NOMYO2 DELTA ;
  1328.   DYHOYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1329.       PCEN9 NOMYO2 NOMYHO DELTA ;
  1330. ***************************************
  1331.   DRYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1332.       PCEN9 NOMYHO NOMDEN DELTA ;  
  1333.   DGXYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1334.       PCEN9 NOMYHO NOMMOX DELTA ;      
  1335.   DGYYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1336.       PCEN9 NOMYHO NOMMOY DELTA ;
  1337.   DRETYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1338.       PCEN9 NOMYHO NOMRET DELTA ;
  1339.   DYH2YHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1340.       PCEN9 NOMYHO NOMYH2 DELTA ;
  1341.   DYO2YHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1342.       PCEN9 NOMYHO NOMYO2 DELTA ;
  1343.   DYHOYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN6
  1344.       PCEN9 NOMYHO NOMYHO DELTA ;      
  1345. **************************************************************
  1346. * Test des comparaisons jacobien exact-jacobien numerique
  1347. ************************************************************** 
  1348.   'SI' (('ABS'(DRR)) < DELTA) ;
  1349.     ERR1 = ('ABS' (DRRN '-' DRR)) ;
  1350.  'SINON' ;  
  1351.     ERR1 = ('ABS' (DRRN '-' DRR)) '/' ('ABS' (DRR)) ;
  1352.  'FINSI' ; 
  1353. ***********************
  1354.  'SI' (ERR1 > ERRTOL) ;
  1355.       'MESSAGE' 'Problem 1b';
  1356.       'ERREUR' 5 ;
  1357.  'FINSI' ;
  1358. *************************************************************
  1359.  'SI' (('ABS'(DGXR)) < DELTA) ;
  1360.     ERR1 = ('ABS' (DGXRN '-' DGXR)) ;
  1361.  'SINON' ;  
  1362.     ERR1 = ('ABS' (DGXRN '-' DGXR)) '/' ('ABS' (DGXR)) ;
  1363.  'FINSI' ; 
  1364. ***********************
  1365.  'SI' (ERR1 > ERRTOL) ;
  1366.       'MESSAGE' 'Problem 2b';
  1367.       'ERREUR' 5 ;
  1368.  'FINSI' ;      
  1369. **************************************************************
  1370.  'SI' (('ABS'(DGYR)) < DELTA) ;
  1371.     ERR1 = ('ABS' (DGYRN '-' DGYR)) ;
  1372.  'SINON' ;  
  1373.     ERR1 = ('ABS' (DGYRN '-' DGYR)) '/' ('ABS' (DGYR)) ;
  1374.  'FINSI' ; 
  1375. ************************
  1376.  'SI' (ERR1 > ERRTOL) ;
  1377.       'MESSAGE' 'Problem 3b';
  1378.        'ERREUR' 5 ;
  1379.  'FINSI' ;
  1380. **************************************************************
  1381.  'SI' (('ABS'(DRETR)) < DELTA) ;
  1382.     ERR1 = ('ABS' (DRETRN '-' DRETR)) ;
  1383.  'SINON' ;  
  1384.     ERR1 = ('ABS' (DRETRN '-' DRETR)) '/' ('ABS' (DRETR)) ;
  1385.  'FINSI' ;
  1386. **************************
  1387.  'SI' (ERR1 > ERRTOL) ;
  1388.       'MESSAGE' 'Problem 4b';
  1389.       'ERREUR' 5 ;
  1390.  'FINSI' ;
  1391. **************************************************************
  1392.  'SI' (('ABS'(DYH2R)) < DELTA) ;
  1393.     ERR1 = ('ABS' (DYH2RN '-' DYH2R)) ;
  1394.  'SINON' ;  
  1395.     ERR1 = ('ABS' (DYH2RN '-' DYH2R)) '/' ('ABS' (DYH2R)) ;
  1396.  'FINSI' ;
  1397. **************************
  1398.  'SI' (ERR1 > ERRTOL) ;
  1399.       'MESSAGE' 'Problem 5b';
  1400.       'ERREUR' 5 ;
  1401.  'FINSI' ;
  1402. **************************************************************
  1403.  'SI' (('ABS'(DYO2R)) < DELTA) ;
  1404.     ERR1 = ('ABS' (DYO2RN '-' DYO2R)) ;
  1405.  'SINON' ;  
  1406.     ERR1 = ('ABS' (DYO2RN '-' DYO2R)) '/' ('ABS' (DYO2R)) ;
  1407.  'FINSI' ;
  1408. **************************
  1409.  'SI' (ERR1 > ERRTOL) ;
  1410.       'MESSAGE' 'Problem 6b';
  1411.       'ERREUR' 5 ;
  1412.  'FINSI' ;
  1413. **************************************************************
  1414.  'SI' (('ABS'(DYHOR)) < DELTA) ;
  1415.     ERR1 = ('ABS' (DYHORN '-' DYHOR)) ;
  1416.  'SINON' ;  
  1417.     ERR1 = ('ABS' (DYHORN '-' DYHOR)) '/' ('ABS' (DYHOR)) ;
  1418.  'FINSI' ;
  1419. **************************
  1420.  'SI' (ERR1 > ERRTOL) ;
  1421.       'MESSAGE' 'Problem 7b';
  1422.       'ERREUR' 5 ;
  1423.  'FINSI' ;
  1424. **************************************************************
  1425. *****    Second line !!!!!!!!!!!!!!!!!!!!!!!    **************
  1426. **************************************************************
  1427. **************************************************************
  1428.  'SI' (('ABS'(DRGX)) < DELTA) ;
  1429.     ERR1 = ('ABS' (DRGXN '-' DRGX)) ;
  1430.  'SINON' ;   
  1431.     ERR1 = ('ABS' (DRGXN '-' DRGX)) '/' ('ABS' (DRGX))  ;
  1432.  'FINSI' ;
  1433. ***************************
  1434.  'SI' (ERR1 > ERRTOL) ;
  1435.       'MESSAGE' 'Problem 8b';
  1436.       'ERREUR' 5 ;
  1437.  'FINSI' ;
  1438. *************************************************************
  1439.  'SI' (('ABS'(DGXGX)) < DELTA) ;
  1440.    ERR1 = ('ABS' (DGXGXN '-' DGXGX)) ;
  1441.  'SINON' ;
  1442.    ERR1 = ('ABS' (DGXGXN '-' DGXGX)) '/' ('ABS' (DGXGX)) ;
  1443.  'FINSI' ;  
  1444. *****************************   
  1445.  'SI' (ERR1 > ERRTOL) ;
  1446.         'MESSAGE' 'Problem 9b';
  1447.         'ERREUR' 5 ;
  1448.  'FINSI' ;
  1449. *************************************************************
  1450.  'SI' (('ABS'(DGYGX)) < DELTA) ;
  1451.     ERR1 = ('ABS' (DGYGXN '-' DGYGX)) ;
  1452.  'SINON' ;
  1453.     ERR1 = ('ABS' (DGYGXN '-' DGYGX)) '/' ('ABS' (DGYGX)) ;
  1454.  'FINSI' ;
  1455. ******************************
  1456.  'SI' (ERR1 > ERRTOL) ;
  1457.         'MESSAGE' 'Problem 10b';
  1458.         'ERREUR' 5 ;
  1459.  'FINSI' ;
  1460. *************************************************************
  1461.  'SI' (('ABS'(DRETGX)) < DELTA) ;
  1462.     ERR1 = ('ABS' (DRETGXN '-' DRETGX)) ;
  1463.  'SINON' ;
  1464.     ERR1 = ('ABS' (DRETGXN '-' DRETGX)) '/' ('ABS' (DRETGX)) ;
  1465.  'FINSI' ;
  1466. *******************************    
  1467.  'SI' (ERR1 > ERRTOL) ;
  1468.         'MESSAGE' 'Problem 11b';
  1469.         'ERREUR' 5 ;
  1470.  'FINSI' ;
  1471. **************************************************************
  1472.  'SI' (('ABS'(DYH2GX)) < DELTA) ;
  1473.     ERR1 = ('ABS' (DYH2GXN '-' DYH2GX)) ;
  1474.  'SINON' ;  
  1475.     ERR1 = ('ABS' (DYH2GXN '-' DYH2GX)) '/' ('ABS' (DYH2GX)) ;
  1476.  'FINSI' ;
  1477. **************************
  1478.  'SI' (ERR1 > ERRTOL) ;
  1479.       'MESSAGE' 'Problem 12b';
  1480.       'ERREUR' 5 ;
  1481.  'FINSI' ;
  1482. **************************************************************
  1483.  'SI' (('ABS'(DYO2GX)) < DELTA) ;
  1484.     ERR1 = ('ABS' (DYO2GXN '-' DYO2GX)) ;
  1485.  'SINON' ;  
  1486.     ERR1 = ('ABS' (DYO2GXN '-' DYO2GX)) '/' ('ABS' (DYO2GX)) ;
  1487.  'FINSI' ;
  1488. **************************
  1489.  'SI' (ERR1 > ERRTOL) ;
  1490.       'MESSAGE' 'Problem 13b';
  1491.       'ERREUR' 5 ;
  1492.  'FINSI' ;
  1493. **************************************************************
  1494.  'SI' (('ABS'(DYHOGX)) < DELTA) ;
  1495.     ERR1 = ('ABS' (DYHOGXN '-' DYHOGX)) ;
  1496.  'SINON' ;  
  1497.     ERR1 = ('ABS' (DYHOGXN '-' DYHOGX)) '/' ('ABS' (DYHOGX)) ;
  1498.  'FINSI' ;
  1499. **************************
  1500.  'SI' (ERR1 > ERRTOL) ;
  1501.       'MESSAGE' 'Problem 14b';
  1502.       'ERREUR' 5 ;
  1503.  'FINSI' ;
  1504. **************************************************************
  1505. *****    Third line !!!!!!!!!!!!!!!!!!!!!!!    **************
  1506. **************************************************************
  1507. **************************************************************
  1508.  'SI' (('ABS'(DRGY)) < DELTA) ;
  1509.     ERR1 = ('ABS' (DRGYN '-' DRGY)) ;
  1510.  'SINON' ;
  1511.     ERR1 = ('ABS' (DRGYN '-' DRGY)) '/' ('ABS' (DRGY)) ;
  1512.  'FINSI' ;
  1513. *******************************
  1514.  'SI' (ERR1 > ERRTOL) ;
  1515.         'MESSAGE' 'Problem 15b';
  1516.         'ERREUR' 5 ;
  1517.  'FINSI' ;
  1518. **************************************************************
  1519.  'SI' (('ABS'(DGXGY)) < DELTA) ;
  1520.     ERR1 = ('ABS' (DGXGYN '-' DGXGY)) ;
  1521.  'SINON' ;
  1522.     ERR1 = ('ABS' (DGXGYN '-' DGXGY))  '/' ('ABS' (DGXGY)) ;
  1523.  'FINSI' ;
  1524. ********************************
  1525.  'SI' (ERR1 > ERRTOL) ;
  1526.         'MESSAGE' 'Problem 16b';
  1527.         'ERREUR' 5 ;
  1528.  'FINSI' ;
  1529. **************************************************************
  1530.  'SI' (('ABS'(DGYGY)) < DELTA) ;
  1531.     ERR1 = ('ABS' (DGYGYN '-' DGYGY)) ;
  1532.  'SINON' ;
  1533.     ERR1 = ('ABS' (DGYGYN '-' DGYGY)) '/' ('ABS' (DGYGY)) ;
  1534.  'FINSI' ;
  1535. ***************
  1536.  'SI' (ERR1 > ERRTOL) ;
  1537.         'MESSAGE' 'Problem 17b';
  1538.         'ERREUR' 5 ;
  1539.  'FINSI' ;
  1540. *************************************************************
  1541.  'SI' (('ABS'(DRETGY)) < DELTA) ;
  1542.     ERR1 = ('ABS' (DRETGYN '-' DRETGY)) ;
  1543.  'SINON' ;
  1544.     ERR1 = ('ABS' (DRETGYN '-' DRETGY)) '/' ('ABS' (DRETGY)) ;
  1545.  'FINSI' ;
  1546. ****************     
  1547.  'SI' (ERR1 > ERRTOL) ;
  1548.         'MESSAGE' 'Problem 18b';
  1549.         'ERREUR' 5 ;
  1550.  'FINSI' ;
  1551. **************************************************************
  1552.  'SI' (('ABS'(DYH2GY)) < DELTA) ;
  1553.     ERR1 = ('ABS' (DYH2GYN '-' DYH2GY)) ;
  1554.  'SINON' ;  
  1555.     ERR1 = ('ABS' (DYH2GYN '-' DYH2GY)) '/' ('ABS' (DYH2GY)) ;
  1556.  'FINSI' ;
  1557. **************************
  1558.  'SI' (ERR1 > ERRTOL) ;
  1559.       'MESSAGE' 'Problem 19b';
  1560.       'ERREUR' 5 ;
  1561.  'FINSI' ;
  1562. **************************************************************
  1563.  'SI' (('ABS'(DYO2GY)) < DELTA) ;
  1564.     ERR1 = ('ABS' (DYO2GYN '-' DYO2GY)) ;
  1565.  'SINON' ;  
  1566.     ERR1 = ('ABS' (DYO2GYN '-' DYO2GY)) '/' ('ABS' (DYO2GY)) ;
  1567.  'FINSI' ;
  1568. **************************
  1569.  'SI' (ERR1 > ERRTOL) ;
  1570.       'MESSAGE' 'Problem 20b';
  1571.       'ERREUR' 5 ;
  1572.  'FINSI' ;
  1573. **************************************************************
  1574.  'SI' (('ABS'(DYHOGY)) < DELTA) ;
  1575.     ERR1 = ('ABS' (DYHOGYN '-' DYHOGY)) ;
  1576.  'SINON' ;  
  1577.     ERR1 = ('ABS' (DYHOGYN '-' DYHOGY)) '/' ('ABS' (DYHOGY)) ;
  1578.  'FINSI' ;
  1579. **************************
  1580.  'SI' (ERR1 > ERRTOL) ;
  1581.       'MESSAGE' 'Problem 21b';
  1582.       'ERREUR' 5 ;
  1583.  'FINSI' ;
  1584. **************************************************************
  1585. *******     Fouth line!!!!!!!!!!!!!!!!!!!   ******************
  1586. **************************************************************
  1587. **************************************************************
  1588.  'SI' (('ABS'(DRRET)) < DELTA) ;
  1589.     ERR1 = ('ABS' (DRRETN '-' DRRET)) ;
  1590.  'SINON' ;
  1591.     ERR1 = ('ABS' (DRRETN '-' DRRET)) '/' ('ABS' (DRRET)) ;
  1592.  'FINSI' ;
  1593. *****************
  1594.  'SI' (ERR1 > ERRTOL) ;
  1595.         'MESSAGE' 'Problem 22b';
  1596.         'ERREUR' 5 ;
  1597.  'FINSI' ;
  1598. **************************************************************
  1599.  'SI' (('ABS'(DGXRET)) < DELTA) ;
  1600.     ERR1 = ('ABS' (DGXRETN '-' DGXRET)) ;
  1601.  'SINON' ;
  1602.     ERR1 = ('ABS' (DGXRETN '-' DGXRET)) '/' ('ABS' (DGXRET)) ;
  1603.  'FINSI' ;
  1604. *****************
  1605.  'SI' (ERR1 > ERRTOL) ;
  1606.         'MESSAGE' 'Problem 23b';
  1607.         'ERREUR' 5 ;
  1608.  'FINSI' ;
  1609. *************************************************************
  1610.  'SI' (('ABS'(DGYRET)) < DELTA) ;
  1611.     ERR1 = ('ABS' (DGYRETN '-' DGYRET)) ; 
  1612.  'SINON' ;
  1613.     ERR1 = ('ABS' (DGYRETN '-' DGYRET)) '/' ('ABS' (DGYRET)) ;
  1614.  'FINSI' ;
  1615. *****************    
  1616.  'SI' (ERR1 > ERRTOL) ;
  1617.         'MESSAGE' 'Problem 24b';
  1618.         'ERREUR' 5 ;
  1619.  'FINSI' ;
  1620. **************************************************************
  1621.  'SI' (('ABS'(DRETRET)) < DELTA) ;
  1622.     ERR1 = ('ABS' (DRETRETN '-' DRETRET)) ;
  1623.  'SINON' ;
  1624.     ERR1 = ('ABS' (DRETRETN '-' DRETRET)) '/' ('ABS' (DRETRET)) ;
  1625.  'FINSI' ;
  1626. *****************
  1627.  'SI' (ERR1 > ERRTOL) ;
  1628.         'MESSAGE' 'Problem 25b';
  1629.         'ERREUR' 5 ;
  1630.  'FINSI' ;
  1631. **************************************************************
  1632.  'SI' (('ABS'(DYH2RET)) < DELTA) ;
  1633.     ERR1 = ('ABS' (DYH2RETN '-' DYH2RET)) ;
  1634.  'SINON' ;  
  1635.     ERR1 = ('ABS' (DYH2RETN '-' DYH2RET)) '/' ('ABS' (DYH2RET)) ;
  1636.  'FINSI' ;
  1637. **************************
  1638.  'SI' (ERR1 > ERRTOL) ;
  1639.       'MESSAGE' 'Problem 26b';
  1640.       'ERREUR' 5 ;
  1641.  'FINSI' ;
  1642. **************************************************************
  1643.  'SI' (('ABS'(DYO2RET)) < DELTA) ;
  1644.     ERR1 = ('ABS' (DYO2RETN '-' DYO2RET)) ;
  1645.  'SINON' ;  
  1646.     ERR1 = ('ABS' (DYO2RETN '-' DYO2RET)) '/' ('ABS' (DYO2RET)) ;
  1647.  'FINSI' ;
  1648. **************************
  1649.  'SI' (ERR1 > ERRTOL) ;
  1650.       'MESSAGE' 'Problem 27b';
  1651.       'ERREUR' 5 ;
  1652.  'FINSI' ;
  1653. **************************************************************
  1654.  'SI' (('ABS'(DYHORET)) < DELTA) ;
  1655.     ERR1 = ('ABS' (DYHORETN '-' DYHORET)) ;
  1656.  'SINON' ;  
  1657.     ERR1 = ('ABS' (DYHORETN '-' DYHORET)) '/' ('ABS' (DYHORET)) ;
  1658.  'FINSI' ;
  1659. **************************
  1660.  'SI' (ERR1 > ERRTOL) ;
  1661.       'MESSAGE' 'Problem 28b';
  1662.       'ERREUR' 5 ;
  1663.  'FINSI' ;
  1664. **************************************************************
  1665. *******     Fith line!!!!!!!!!!!!!!!!!!!   ******************
  1666. **************************************************************
  1667.  'SI' (('ABS'(DRYH2)) < DELTA) ;
  1668.     ERR1 = ('ABS' (DRYH2N '-' DRYH2)) ;
  1669.  'SINON' ;
  1670.     ERR1 = ('ABS' (DRYH2N '-' DRYH2)) '/' ('ABS' (DRYH2)) ;
  1671.  'FINSI' ;
  1672. *****************
  1673.  'SI' (ERR1 > ERRTOL) ;
  1674.         'MESSAGE' 'Problem 29b';
  1675.         'ERREUR' 5 ;
  1676.  'FINSI' ;
  1677. **************************************************************
  1678.  'SI' (('ABS'(DGXYH2)) < DELTA) ;
  1679.     ERR1 = ('ABS' (DGXYH2N '-' DGXYH2)) ;
  1680.  'SINON' ;
  1681.     ERR1 = ('ABS' (DGXYH2N '-' DGXYH2)) '/' ('ABS' (DGXYH2)) ;
  1682.  'FINSI' ;
  1683. *****************
  1684.  'SI' (ERR1 > ERRTOL) ;
  1685.         'MESSAGE' 'Problem 30b';
  1686.         'ERREUR' 5 ;
  1687.  'FINSI' ;
  1688. *************************************************************
  1689.  'SI' (('ABS'(DGYYH2)) < DELTA) ;
  1690.     ERR1 = ('ABS' (DGYYH2N '-' DGYYH2)) ; 
  1691.  'SINON' ;
  1692.     ERR1 = ('ABS' (DGYYH2N '-' DGYYH2)) '/' ('ABS' (DGYYH2)) ;
  1693.  'FINSI' ;
  1694. *****************    
  1695.  'SI' (ERR1 > ERRTOL) ;
  1696.         'MESSAGE' 'Problem 31b';
  1697.         'ERREUR' 5 ;
  1698.  'FINSI' ;
  1699. **************************************************************
  1700.  'SI' (('ABS'(DRETYH2)) < DELTA) ;
  1701.     ERR1 = ('ABS' (DRETYH2N '-' DRETYH2)) ;
  1702.  'SINON' ;
  1703.     ERR1 = ('ABS' (DRETYH2N '-' DRETYH2)) '/' ('ABS' (DRETYH2)) ;
  1704.  'FINSI' ;
  1705. *****************
  1706.  'SI' (ERR1 > ERRTOL) ;
  1707.         'MESSAGE' 'Problem 32b';
  1708.         'ERREUR' 5 ;
  1709.  'FINSI' ;
  1710. **************************************************************
  1711.  'SI' (('ABS'(DYH2YH2)) < DELTA) ;
  1712.     ERR1 = ('ABS' (DYH2YH2N '-' DYH2YH2)) ;
  1713.  'SINON' ;  
  1714.     ERR1 = ('ABS' (DYH2YH2N '-' DYH2YH2)) '/' ('ABS' (DYH2YH2)) ;
  1715.  'FINSI' ;
  1716. **************************
  1717.  'SI' (ERR1 > ERRTOL) ;
  1718.       'MESSAGE' 'Problem 33b';
  1719.       'ERREUR' 5 ;
  1720.  'FINSI' ;
  1721. **************************************************************
  1722.  'SI' (('ABS'(DYO2YH2)) < DELTA) ;
  1723.     ERR1 = ('ABS' (DYO2YH2N '-' DYO2YH2)) ;
  1724.  'SINON' ;  
  1725.     ERR1 = ('ABS' (DYO2YH2N '-' DYO2YH2)) '/' ('ABS' (DYO2YH2)) ;
  1726.  'FINSI' ;
  1727. **************************
  1728.  'SI' (ERR1 > ERRTOL) ;
  1729.       'MESSAGE' 'Problem 34b';
  1730.       'ERREUR' 5 ;
  1731.  'FINSI' ;
  1732. **************************************************************
  1733.  'SI' (('ABS'(DYHOYH2)) < DELTA) ;
  1734.     ERR1 = ('ABS' (DYHOYH2N '-' DYHOYH2)) ;
  1735.  'SINON' ;  
  1736.     ERR1 = ('ABS' (DYHOYH2N '-' DYHOYH2)) '/' ('ABS' (DYHOYH2)) ;
  1737.  'FINSI' ;
  1738. **************************
  1739.  'SI' (ERR1 > ERRTOL) ;
  1740.       'MESSAGE' 'Problem 35b';
  1741.       'ERREUR' 5 ;
  1742.  'FINSI' ;
  1743. **************************************************************
  1744. *********  Sixth line !!!!!!!!!!!!!!!!!!!!!!   ***************
  1745. **************************************************************
  1746.  'SI' (('ABS'(DRYO2)) < DELTA) ;
  1747.     ERR1 = ('ABS' (DRYO2N '-' DRYO2)) ;
  1748.  'SINON' ;
  1749.     ERR1 = ('ABS' (DRYO2N '-' DRYO2)) '/' ('ABS' (DRYO2)) ;
  1750.  'FINSI' ;
  1751. *****************
  1752.  'SI' (ERR1 > ERRTOL) ;
  1753.         'MESSAGE' 'Problem 36b';
  1754.         'ERREUR' 5 ;
  1755.  'FINSI' ;
  1756. **************************************************************
  1757.  'SI' (('ABS'(DGXYO2)) < DELTA) ;
  1758.     ERR1 = ('ABS' (DGXYO2N '-' DGXYO2)) ;
  1759.  'SINON' ;
  1760.     ERR1 = ('ABS' (DGXYO2N '-' DGXYO2)) '/' ('ABS' (DGXYO2)) ;
  1761.  'FINSI' ;
  1762. *****************
  1763.  'SI' (ERR1 > ERRTOL) ;
  1764.         'MESSAGE' 'Problem 37b';
  1765.         'ERREUR' 5 ;
  1766.  'FINSI' ;
  1767. *************************************************************
  1768.  'SI' (('ABS'(DGYYO2)) < DELTA) ;
  1769.     ERR1 = ('ABS' (DGYYO2N '-' DGYYO2)) ; 
  1770.  'SINON' ;
  1771.     ERR1 = ('ABS' (DGYYO2N '-' DGYYO2)) '/' ('ABS' (DGYYO2)) ;
  1772.  'FINSI' ;
  1773. *****************    
  1774.  'SI' (ERR1 > ERRTOL) ;
  1775.         'MESSAGE' 'Problem 38b';
  1776.         'ERREUR' 5 ;
  1777.  'FINSI' ;
  1778. **************************************************************
  1779.  'SI' (('ABS'(DRETYO2)) < DELTA) ;
  1780.     ERR1 = ('ABS' (DRETYO2N '-' DRETYO2)) ;
  1781.  'SINON' ;
  1782.     ERR1 = ('ABS' (DRETYO2N '-' DRETYO2)) '/' ('ABS' (DRETYO2)) ;
  1783.  'FINSI' ;
  1784. *****************
  1785.  'SI' (ERR1 > ERRTOL) ;
  1786.         'MESSAGE' 'Problem 39b';
  1787.         'ERREUR' 5 ;
  1788.  'FINSI' ;
  1789. **************************************************************
  1790.  'SI' (('ABS'(DYH2YO2)) < DELTA) ;
  1791.     ERR1 = ('ABS' (DYH2YO2N '-' DYH2YO2)) ;
  1792.  'SINON' ;  
  1793.     ERR1 = ('ABS' (DYH2YO2N '-' DYH2YO2)) '/' ('ABS' (DYH2YO2)) ;
  1794.  'FINSI' ;
  1795. **************************
  1796.  'SI' (ERR1 > ERRTOL) ;
  1797.       'MESSAGE' 'Problem 40b';
  1798.       'ERREUR' 5 ;
  1799.  'FINSI' ;
  1800. **************************************************************
  1801.  'SI' (('ABS'(DYO2YO2)) < DELTA) ;
  1802.     ERR1 = ('ABS' (DYO2YO2N '-' DYO2YO2)) ;
  1803.  'SINON' ;  
  1804.     ERR1 = ('ABS' (DYO2YO2N '-' DYO2YO2)) '/' ('ABS' (DYO2YO2)) ;
  1805.  'FINSI' ;
  1806. **************************
  1807.  'SI' (ERR1 > ERRTOL) ;
  1808.       'MESSAGE' 'Problem 41b';
  1809.       'ERREUR' 5 ;
  1810.  'FINSI' ;
  1811. **************************************************************
  1812.  'SI' (('ABS'(DYHOYO2)) < DELTA) ;
  1813.     ERR1 = ('ABS' (DYHOYO2N '-' DYHOYO2)) ;
  1814.  'SINON' ;  
  1815.     ERR1 = ('ABS' (DYHOYO2N '-' DYHOYO2)) '/' ('ABS' (DYHOYO2)) ;
  1816.  'FINSI' ;
  1817. **************************
  1818.  'SI' (ERR1 > ERRTOL) ;
  1819.       'MESSAGE' 'Problem 42b';
  1820.       'ERREUR' 5 ;
  1821.  'FINSI' ;
  1822. **************************************************************
  1823. *********  Seventh line !!!!!!!!!!!!!!!!!!!!!!   *************
  1824. **************************************************************
  1825.  'SI' (('ABS'(DRYHO)) < DELTA) ;
  1826.     ERR1 = ('ABS' (DRYHON '-' DRYHO)) ;
  1827.  'SINON' ;
  1828.     ERR1 = ('ABS' (DRYHON '-' DRYHO)) '/' ('ABS' (DRYHO)) ;
  1829.  'FINSI' ;
  1830. *****************
  1831.  'SI' (ERR1 > ERRTOL) ;
  1832.         'MESSAGE' 'Problem 43b';
  1833.         'ERREUR' 5 ;
  1834.  'FINSI' ;
  1835. **************************************************************
  1836.  'SI' (('ABS'(DGXYHO)) < DELTA) ;
  1837.     ERR1 = ('ABS' (DGXYHON '-' DGXYHO)) ;
  1838.  'SINON' ;
  1839.     ERR1 = ('ABS' (DGXYHON '-' DGXYHO)) '/' ('ABS' (DGXYHO)) ;
  1840.  'FINSI' ;
  1841. *****************
  1842.  'SI' (ERR1 > ERRTOL) ;
  1843.         'MESSAGE' 'Problem 44b';
  1844.         'ERREUR' 5 ;
  1845.  'FINSI' ;
  1846. *************************************************************
  1847.  'SI' (('ABS'(DGYYHO)) < DELTA) ;
  1848.     ERR1 = ('ABS' (DGYYHON '-' DGYYHO)) ; 
  1849.  'SINON' ;
  1850.     ERR1 = ('ABS' (DGYYHON '-' DGYYHO)) '/' ('ABS' (DGYYHO)) ;
  1851.  'FINSI' ;
  1852. *****************    
  1853.  'SI' (ERR1 > ERRTOL) ;
  1854.         'MESSAGE' 'Problem 45b';
  1855.         'ERREUR' 5 ;
  1856.  'FINSI' ;
  1857. **************************************************************
  1858.  'SI' (('ABS'(DRETYHO)) < DELTA) ;
  1859.     ERR1 = ('ABS' (DRETYHON '-' DRETYHO)) ;
  1860.  'SINON' ;
  1861.     ERR1 = ('ABS' (DRETYHON '-' DRETYHO)) '/' ('ABS' (DRETYHO)) ;
  1862.  'FINSI' ;
  1863. *****************
  1864.  'SI' (ERR1 > ERRTOL) ;
  1865.         'MESSAGE' 'Problem 46b';
  1866.         'ERREUR' 5 ;
  1867.  'FINSI' ;
  1868. **************************************************************
  1869.  'SI' (('ABS'(DYH2YHO)) < DELTA) ;
  1870.     ERR1 = ('ABS' (DYH2YHON '-' DYH2YHO)) ;
  1871.  'SINON' ;  
  1872.     ERR1 = ('ABS' (DYH2YHON '-' DYH2YHO)) '/' ('ABS' (DYH2YHO)) ;
  1873.  'FINSI' ;
  1874. **************************
  1875.  'SI' (ERR1 > ERRTOL) ;
  1876.       'MESSAGE' 'Problem 47b';
  1877.       'ERREUR' 5 ;
  1878.  'FINSI' ;
  1879. **************************************************************
  1880.  'SI' (('ABS'(DYO2YHO)) < DELTA) ;
  1881.     ERR1 = ('ABS' (DYO2YHON '-' DYO2YHO)) ;
  1882.  'SINON' ;  
  1883.     ERR1 = ('ABS' (DYO2YHON '-' DYO2YHO)) '/' ('ABS' (DYO2YHO)) ;
  1884.  'FINSI' ;
  1885. **************************
  1886.  'SI' (ERR1 > ERRTOL) ;
  1887.       'MESSAGE' 'Problem 48b';
  1888.       'ERREUR' 5 ;
  1889.  'FINSI' ;
  1890. **************************************************************
  1891.  'SI' (('ABS'(DYHOYHO)) < DELTA) ;
  1892.     ERR1 = ('ABS' (DYHOYHON '-' DYHOYHO)) ;
  1893.  'SINON' ;  
  1894.     ERR1 = ('ABS' (DYHOYHON '-' DYHOYHO)) '/' ('ABS' (DYHOYHO)) ;
  1895.  'FINSI' ;
  1896. **************************
  1897.  'SI' (ERR1 > ERRTOL) ;
  1898.       'MESSAGE' 'Problem 49b';
  1899.       'ERREUR' 5 ;
  1900.  'FINSI' ;
  1901. **********************************************************
  1902. *****************************************************
  1903. *****************************************************
  1904. ******* TEST3 ***************************************
  1905. *****************************************************
  1906. *****************************************************
  1907. *
  1908. * On observe que la variation du residu
  1909. * en $DOM9 . 'CENTRE' par rapport à une variation
  1910. * infinitésimal en $DOM7 . 'CENTRE' doit etre nulle.
  1911. * 
  1912.  
  1913.  PCEN7 = 'POIN' 1 ('DOMA' $DOM7 CENTRE) ;
  1914.  
  1915. *
  1916. * Le jacobien exact. 
  1917. *
  1918.  
  1919.   DRR = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMDEN NOMDEN ;
  1920.   DGXR = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMDEN NOMMOX ;
  1921.   DGYR = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMDEN NOMMOY ;
  1922.   DRETR = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMDEN NOMRET ;
  1923.   DYH2R = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMDEN NOMYH2 ;
  1924.   DYO2R = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMDEN NOMYO2 ;
  1925.   DYHOR = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMDEN NOMYHO ;
  1926. ***************************************
  1927.   DRGX = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOX NOMDEN ;
  1928.   DGXGX = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOX NOMMOX ;
  1929.   DGYGX = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOX NOMMOY ;
  1930.   DRETGX = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOX NOMRET ;
  1931.   DYH2GX = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOX NOMYH2 ;
  1932.   DYO2GX = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOX NOMYO2 ;
  1933.   DYHOGX = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOX NOMYHO ;
  1934. ***************************************
  1935.   DRGY = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOY NOMDEN ;
  1936.   DGXGY = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOY NOMMOX ;
  1937.   DGYGY = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOY NOMMOY ;
  1938.   DRETGY = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOY NOMRET ;
  1939.   DYH2GY = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOY NOMYH2 ;
  1940.   DYO2GY = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOY NOMYO2 ;
  1941.   DYHOGY = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMMOY NOMYHO ;  
  1942. ***************************************
  1943.   DRRET = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMRET NOMDEN ;
  1944.   DGXRET = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMRET NOMMOX ;
  1945.   DGYRET = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMRET NOMMOY ;
  1946.   DRETRET = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMRET NOMRET ;
  1947.   DYH2RET = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMRET NOMYH2 ;
  1948.   DYO2RET = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMRET NOMYO2 ;
  1949.   DYHORET = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMRET NOMYHO ;
  1950. ***************************************
  1951.   DRYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYH2 NOMDEN ;
  1952.   DGXYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYH2 NOMMOX ;
  1953.   DGYYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYH2 NOMMOY ;
  1954.   DRETYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYH2 NOMRET ;
  1955.   DYH2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYH2 NOMYH2 ;
  1956.   DYO2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYH2 NOMYO2 ;
  1957.   DYHOYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYH2 NOMYHO ;
  1958. ***************************************
  1959.   DRYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYO2 NOMDEN ;
  1960.   DGXYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYO2 NOMMOX ;
  1961.   DGYYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYO2 NOMMOY ;
  1962.   DRETYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYO2 NOMRET ;
  1963.   DYH2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYO2 NOMYH2 ;
  1964.   DYO2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYO2 NOMYO2 ;
  1965.   DYHOYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYO2 NOMYHO ;
  1966. ****************************************
  1967.   DRYHO = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYHO NOMDEN ;
  1968.   DGXYHO = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYHO NOMMOX ;
  1969.   DGYYHO = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYHO NOMMOY ;
  1970.   DRETYHO = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYHO NOMRET ;
  1971.   DYH2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYHO NOMYH2 ;
  1972.   DYO2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYHO NOMYO2 ;
  1973.   DYHOYHO = JAEXAC RMAT $DOMTOT LINCO PCEN7 PCEN9 NOMYHO NOMYHO ;
  1974.  
  1975. **********************************************
  1976.  'SI' (('MAXIMUM'
  1977.     ('PROG'  DRR  DGXR DGYR DRETR
  1978.              DYH2R DYO2R DYHOR
  1979.              DRGX DGXGX DGYGX DRETGX
  1980.              DYH2GX DYO2GX DYHOGX ) 'ABS' )
  1981.       > ERRTOL) ;
  1982.       'MESSAGE' 'Problem final 1' ;
  1983.       'ERREUR' 5 ;
  1984.  'FINSI' ;    
  1985. *********************************************
  1986.  'SI' (('MAXIMUM'
  1987.     ('PROG'  DRGY DGXGY DGYGY DRETGY
  1988.              DYH2GY DYO2GY DYHOGY
  1989.              DRRET DGXRET DGYRET DRETRET
  1990.              DYH2RET DYO2RET DYHORET ) 'ABS' )
  1991.       > ERRTOL) ;
  1992.       'MESSAGE' 'Problem final 2' ;
  1993.       'ERREUR' 5 ;
  1994.  'FINSI' ;
  1995. *********************************************                  
  1996.   'SI' (('MAXIMUM'
  1997.     ('PROG'  DRYH2 DGXYH2 DGYYH2 DRETYH2
  1998.              DYH2YH2 DYO2YH2 DYHOYH2
  1999.              DRYO2 DGXYO2 DGYYO2 DRETYO2
  2000.              DYH2YO2 DYO2YO2 DYHOYO2 ) 'ABS' )
  2001.       > ERRTOL) ;
  2002.       'MESSAGE' 'Problem final 3' ;
  2003.       'ERREUR' 5 ;
  2004.  'FINSI' ;
  2005. **********************************************
  2006.   'SI' (('MAXIMUM'
  2007.     ('PROG'  DRYHO DGXYHO DGYYHO DRETYHO
  2008.              DYH2YHO DYO2YHO DYHOYHO ) 'ABS' )
  2009.       > ERRTOL) ;
  2010.       'MESSAGE' 'Problem final 4' ;
  2011.       'ERREUR' 5 ;
  2012.  'FINSI' ;
  2013. **********************************************
  2014.  
  2015.  'FIN' ;
  2016.  
  2017.  
  2018.  
  2019.  
  2020.  
  2021.  
  2022.  
  2023.  
  2024.  
  2025.  
  2026.  

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