Télécharger la plaquette Cast3M


fr en

Télécharger konmsp_impl3D.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  konmsp_impl3D.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ***********************************************************
  5. ***********************************************************
  6. **** APPROCHE VF "Cell-Centred Formulation" pour la    ****
  7. **** solution des                                      **** 
  8. **** Equations d'Euler pour un gaz                     ****
  9. ****          parfait multiespecies                    ****
  10. **** OPERATEURS PRIM, PRET, KONV                       ****
  11. **** Implicit: calcul du jacobien du residu            ****
  12. ****                                                   ****
  13. **** Cas gaz multiespece, "calorically perfect"        ****
  14. **** 3D                                                ****
  15. ****                                                   ****
  16. **** Methodes: AUSM+                                   ****
  17. ****                                                   **** 
  18. **** This procedure compares the jacobian exact        ****
  19. **** with jacobian coputed numerically                 ****
  20. **** (using numerical flux function)                   ****
  21. **** S.KUDRIAKOV SFME/LTMF               DEC  2002     ****
  22. ***********************************************************
  23. ***********************************************************
  24.  
  25.  'OPTION'  'DIME' 3 
  26.      'ELEM' 'CUB8' 
  27.      'ECHO'  0
  28.      'TRAC' 'X' ;
  29.  
  30. *
  31. *** GRAPH
  32. *
  33.  
  34.  GRAPH = VRAI ;
  35.  GRAPH = FAUX ;
  36.  
  37.  ERRTOL = 3.0D-4 ;
  38.  DELTA = 0.000001 ;
  39.  METO = 'AUSMPLUS' ; 
  40. ***************************
  41. ***** DOMAINE SPATIAL  ****
  42. ***************************
  43.  
  44.  
  45.  A0 = 0.0D0 0.0D0 0.0D0;
  46.  A1 = 1.0D0 0.0D0 0.0D0;
  47.  A2 = 1.0D0 1.0D0 0.0D0;
  48.  A3 = 0.0D0 1.0D0 0.0D0;
  49.  
  50.  SUR1 = 'MANUEL' 'QUA4' A0 A1 A2 A3 ;
  51.  DOM1 = SUR1 'VOLUME' 'TRAN' 1 (0.0 0.0 1.0) ;
  52.  
  53.  DOM2 = DOM1 'PLUS' (1.0 0.0 0.0) ;
  54.  DOM3 = DOM1 'PLUS' (-1.0 0.0 0.0) ;
  55.  DOM4 = DOM1 'PLUS' (0.0 1.0 0.0) ;
  56.  DOM5 = DOM1 'PLUS' (0.0 -1.0 0.0) ;
  57.  DOM6 = DOM1 'PLUS' (1.0 1.0 0.0) ;
  58.  DOM7 = DOM1 'PLUS' (-1.0 -1.0 0.0) ;
  59.  DOM8 = DOM1 'PLUS' (1.0 -1.0 0.0) ;
  60.  DOM9 = DOM1 'PLUS' (-1.0 1.0 0.0) ;
  61.  
  62.  DOM10 = DOM1 'PLUS' (0.0 0.0 1.0) ;
  63.  DOM11 = DOM1 'PLUS' (1.0 0.0 1.0) ;
  64.  DOM12 = DOM1 'PLUS' (-1.0 0.0 1.0) ;
  65.  DOM13 = DOM1 'PLUS' (0.0 1.0 1.0) ;
  66.  DOM14 = DOM1 'PLUS' (0.0 -1.0 1.0) ;
  67.  DOM15 = DOM1 'PLUS' (1.0 1.0 1.0) ;
  68.  DOM16 = DOM1 'PLUS' (-1.0 -1.0 1.0) ;
  69.  DOM17 = DOM1 'PLUS' (1.0 -1.0 1.0) ;
  70.  DOM18 = DOM1 'PLUS' (-1.0 1.0 1.0) ;
  71.  
  72.  DOM19 = DOM1 'PLUS' (0.0 0.0 -1.0) ;
  73.  DOM20 = DOM1 'PLUS' (1.0 0.0 -1.0) ;
  74.  DOM21 = DOM1 'PLUS' (-1.0 0.0 -1.0) ;
  75.  DOM22 = DOM1 'PLUS' (0.0 1.0 -1.0) ;
  76.  DOM23 = DOM1 'PLUS' (0.0 -1.0 -1.0) ;
  77.  DOM24 = DOM1 'PLUS' (1.0 1.0 -1.0) ;
  78.  DOM25 = DOM1 'PLUS' (-1.0 -1.0 -1.0) ;
  79.  DOM26 = DOM1 'PLUS' (1.0 -1.0 -1.0) ;
  80.  DOM27 = DOM1 'PLUS' (-1.0 1.0 -1.0) ;
  81.  
  82.  DOMTOT = DOM1 'ET' DOM2 'ET' DOM3 'ET' DOM4 'ET' DOM5
  83.    'ET' DOM6 'ET' DOM7 'ET' DOM8 'ET' DOM9 'ET' DOM10
  84.    'ET' DOM11 'ET' DOM12 'ET' DOM13 'ET' DOM14 'ET' DOM15
  85.    'ET' DOM16 'ET' DOM17 'ET' DOM18 'ET' DOM19 'ET' DOM20
  86.    'ET' DOM21 'ET' DOM22 'ET' DOM23 'ET' DOM24 'ET' DOM25
  87.    'ET' DOM26 'ET' DOM27 ;
  88.  
  89.  'ELIMINATION' DOMTOT 0.0001 ;
  90.  
  91. *
  92. **** Perturbation du domaine
  93. *
  94.  
  95.  CHPBRU = ('NOMC' 'UX' ('BRUI' 'BLAN' 'UNIF' 0.0  0.1 DOMTOT)
  96.             'NATU' 'DISCRET') 'ET'
  97.             ('NOMC' 'UY' ('BRUI' 'BLAN' 'UNIF' 0.0  0.15 DOMTOT)
  98.             'NATU' 'DISCRET') 'ET'
  99.             ('NOMC' 'UZ' ('BRUI' 'BLAN' 'UNIF' 0.0  0.13 DOMTOT)
  100.             'NATU' 'DISCRET');
  101.  
  102.  'FORME' CHPBRU ;
  103.  
  104.  
  105.  $DOMTOT = 'MODELISER' DOMTOT 'EULER';
  106.  
  107.  $DOM1  = 'MODELISER'  DOM1  'EULER';
  108.  $DOM27 = 'MODELISER'  DOM27 'EULER';
  109.  
  110.  TDOMTOT = 'DOMA'  $DOMTOT 'VF';
  111.  TDOM1  = 'DOMA'  $DOM1  'VF';
  112.  TDOM27 = 'DOMA'  $DOM27 'VF';
  113.  
  114.  MDOMTOT = TDOMTOT . 'QUAF' ;
  115.  MDOM1 = TDOM1 . 'QUAF' ;
  116.  MDOM27 = TDOM27 . 'QUAF' ;
  117.  
  118. **** old stuff  $DOMTOT = 'DOMA' DOMTOT ;
  119.  
  120.  'ELIMINATION'  (MDOMTOT ET MDOM1) 0.0001 ;
  121.  'ELIMINATION'  (MDOMTOT ET MDOM27) 0.0001 ;
  122.  
  123.  MCEN = 'DOMA' $DOMTOT CENTRE ;
  124.  
  125. ********************************************************
  126. *******    MODELE DU GAS   ***************************** 
  127. ********************************************************
  128.  
  129.  NESP = 4;
  130.  
  131. *
  132. *** GAS: H_2, O_2, H_2O, N_2
  133. *
  134. *   CP, CV en J/Kg/K @ T = 3000
  135. *
  136.  
  137.  PGAS = 'TABLE' ;
  138.  
  139.  PGAS . 'CP' = 'TABLE' ;
  140.  PGAS . 'CP' . 'H2  '  = .18729066D+05 ;
  141.  PGAS . 'CP' . 'O2  '  = .11886820D+04 ;
  142.  PGAS . 'CP' . 'H2O '  = .31209047D+04 ;
  143.  PGAS . 'CP' . 'N2  '  = .12993995D+04 ;
  144.  
  145.  
  146.  
  147.  PGAS . 'CV' = 'TABLE' ;
  148.  PGAS . 'CV' . 'H2  '  = .14571861D+05 ;
  149.  PGAS . 'CV' . 'O2  '  = .92885670D+03 ;
  150.  PGAS . 'CV' . 'H2O '  = .26589930D+04 ;
  151.  PGAS . 'CV' . 'N2  '  =  .10024563D+04;
  152.  
  153.  
  154.  
  155. *
  156. **** Especes qui sont dans les equations d'Euler
  157. *
  158.  
  159.  PGAS . 'ESPEULE' = 'MOTS' 'H2  ' 'O2  ' 'H2O ' ;
  160.  
  161. *
  162. **** Espece qui n'y est pas
  163. *
  164.  
  165.  
  166.  PGAS . 'ESPNEULE' = 'N2  ';
  167.  
  168.  
  169. ***************************************************
  170. ***** Densité, pression, vitesse, gamma ***********
  171. ***************************************************
  172.  
  173.  RN = 'BRUI' 'BLAN' 'UNIF' 1.11 0.5 ('DOMA' $DOMTOT 'CENTRE') ;
  174.  PN = 'BRUI' 'BLAN' 'UNIF' 123.41 80 ('DOMA' $DOMTOT 'CENTRE') ;
  175.  
  176. ***  CP  ****
  177.  CP = ('MANUEL' 'CHPO' MCEN 1 'H2'
  178.        (PGAS . 'CP' . 'H2  ') 'NATU' 'DISCRET') 'ET' 
  179.       ('MANUEL' 'CHPO' MCEN 1 'O2'
  180.       (PGAS . 'CP' . 'O2  ') 'NATU' 'DISCRET') 'ET'
  181.        ('MANUEL' 'CHPO' MCEN  1 'H2O'
  182.        (PGAS . 'CP' . 'H2O ') 'NATU' 'DISCRET') 'ET'
  183.        ('MANUEL' 'CHPO' MCEN  1 'N2'
  184.        (PGAS . 'CP' . 'N2  ') 'NATU' 'DISCRET') ;
  185. ***  CV ****
  186.  CV = ('MANUEL' 'CHPO' MCEN 1 'H2'
  187.        (PGAS . 'CV' . 'H2  ') 'NATU' 'DISCRET') 'ET' 
  188.       ('MANUEL' 'CHPO' MCEN 1 'O2'
  189.       (PGAS . 'CV' . 'O2  ') 'NATU' 'DISCRET') 'ET'
  190.        ('MANUEL' 'CHPO' MCEN  1 'H2O'
  191.        (PGAS . 'CV' . 'H2O ') 'NATU' 'DISCRET') 'ET'
  192.        ('MANUEL' 'CHPO' MCEN  1 'N2'
  193.        (PGAS . 'CV' . 'N2  ') 'NATU' 'DISCRET') ;       
  194. *** mass fractions ****
  195.  FF = 'BRUI' 'BLAN' 'UNIF' 0.2 0.05 MCEN  ;
  196.  CHF = 'NOMC' 'H2' FF 'NATU' 'DISCRET' ; 
  197.  
  198.  FS = 'BRUI' 'BLAN' 'UNIF' 0.3 0.04 MCEN ;
  199.  CHS = 'NOMC' 'O2' FS 'NATU' 'DISCRET' ; 
  200.  
  201.  FT = 'BRUI' 'BLAN' 'UNIF' 0.15 0.045 MCEN ;
  202.  CHT = 'NOMC' 'H2O' FT 'NATU' 'DISCRET' ; 
  203.  
  204.  YNB = CHF 'ET' CHS ;
  205.  YN = YNB 'ET' CHT ;
  206.  
  207.  UNCH =  'MANUEL' 'CHPO' MCEN 1 'SCAL' 1.0 ;
  208.  
  209.  YNESP = UNCH '-' ('PSCAL' UNCH YN ('MOTS' 'SCAL' 'SCAL' 'SCAL')
  210.          (PGAS . 'ESPEULE')) ;
  211.  YNESP = 'NOMC' 'N2  ' YNESP 'NATU' 'DISCRET' ;
  212.  
  213.  YNTOT = YN 'ET' YNESP ;
  214.  
  215.  RYN = YN '*' RN  (PGAS . 'ESPEULE') ('MOTS' 'SCAL' 'SCAL' 'SCAL')
  216.        (PGAS . 'ESPEULE') ;
  217. *** gamma  ***************
  218.  LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
  219.  CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
  220.  CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
  221.  GAMN = CPTOT '/' CVTOT ;
  222. **** velocities ********** 
  223.  CSONN = (GAMN '*' PN) '/' RN ;
  224.  CSONN = 'KOPS' CSONN '**' 0.5 ; 
  225.  UXN  = 1.8 *  CSONN ;
  226.  UYN  = 0.1  * CSONN ;
  227.  UZN  = -0.2  * CSONN ;
  228.  
  229. ***************************************************
  230. *************************************************** 
  231.  GN = ('NOMC' (RN '*' UXN) 'UX' 'NATU' 'DISCRET') 'ET'
  232.        ('NOMC' (RN '*' UYN) 'UY' 'NATU' 'DISCRET') 'ET'
  233.        ('NOMC' (RN '*' UZN) 'UZ' 'NATU' 'DISCRET') ;
  234.  
  235.  ECIN = 0.5D0 '*' RN '*' ((UXN  '*' UXN) '+' (UYN  '*' UYN) '+'
  236.       (UZN * UZN)) ;
  237.  RETN  = 'NOMC' 'SCAL' ((PN  '/' (GAMN  '-' 1.0)) '+' ECIN)
  238.     'NATU' 'DISCRET' ;
  239.  
  240.  
  241.  VIT PRES TTEMP FRAM GAMM  = 'PRIM' 'PERFMULT' PGAS RN GN RETN RYN ;
  242.  
  243.  
  244.  ERR1 = 'MAXIMUM' (PRES '-' PN) 'ABS' ;
  245.  ERR2 = 'MAXIMUM' (GN '-' (RN '*' VIT ('MOTS' 'SCAL' 'SCAL' 'SCAL')
  246.          ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) 'ABS' ;
  247.  
  248.  
  249.  'SI' (('MAXIMUM' ('PROG' ERR1 ERR2 ) 'ABS' ) > 1.0D-6) ;
  250.       'MESSAGE' 'Problem 0' ;
  251.       'ERREUR' 5 ;
  252.  'FINSI' ;
  253.  
  254.  
  255.  'SI' GRAPH;
  256.       'TRACER' (('DOMA' $DOMTOT 'MAILLAGE') 'ET' 
  257.           ('DOMA' $DOMTOT 'CENTRE')) 'TITRE' 'Domaine et centre' ;
  258.  'FINSI' ;
  259.  
  260. ****************************************************************
  261. *********  Application of the Operator 'KONV' ******************
  262. ****************************************************************
  263.  NOMDEN = 'RN  ' ;
  264.  NOMMOX = 'RUXN' ;
  265.  NOMMOY = 'RUYN' ;
  266.  NOMMOZ = 'RUZN' ;
  267.  NOMRET = 'RETN' ;
  268.  NOMYH2 = 'RYH2' ;
  269.  NOMYO2 = 'RYO2' ;
  270.  NOMYHO = 'RYHO' ;
  271.  
  272.  LINCO = 'MOTS'  NOMDEN  NOMMOX  NOMMOY NOMMOZ  NOMRET
  273.              NOMYH2 NOMYO2  NOMYHO ;
  274. ******************************************************************
  275.  MRN MUN MPN MYN MGAMN = 'PRET' 'PERFMULT' 1 1 $DOMTOT
  276.                     RN VIT PN YN GAMN ;
  277.  
  278.  RRESID RFLOT = 'KONV' 'VF' 'PERFMULT' 'RESI' METO LINCO
  279.                 $DOMTOT MRN MUN MPN MGAMN MYN PGAS ;                   
  280.  
  281. ******************************************************************
  282. ******************************************************************             
  283.  RMAT = 'KONV' 'VF' 'PERFMULT' 'JACOCONS'
  284.         METO  $DOMTOT PGAS LINCO RN VIT PN YN ;
  285.  
  286. ****************************************************
  287. * Creating a mesh WITHOUT mesh
  288. ****************************************************
  289.   P1 = 0.5 0.5 0.5;
  290.   ELP1 = 'MANU' 'POI1' P1 ;
  291.   MELLIM = 'DIFF' ELP1 ELP1 ;
  292.  
  293.  
  294. *
  295. **** Les variables conservative
  296. *
  297. *    RN0 (densité)
  298. *    GN0 (quantité de mouvement)
  299. *    RETN0 (énergie totale par unité de volume)
  300. *
  301. *    sont definiés
  302. *
  303.  
  304. *****************************************************
  305. *****************************************************
  306. *****************************************************
  307. ******** PROCEDURES *********************************
  308. *****************************************************
  309. *****************************************************
  310. *****************************************************************
  311. ****             Beggining of the procedure            **********
  312. ***************************************************************** 
  313.  
  314.  'DEBPROC' JAEXAC RMAT*'MATRIK' $DOMA*'MMODEL' LISTINCO*'LISTMOTS'
  315.     PPRIM*'POINT' PDUAL*'POINT' MOTPRI*'MOT' MOTDUA*'MOT' ;
  316.  
  317.    MCEN = 'DOMA' $DOMA CENTRE ; 
  318.    ELT1 = 'MANUEL' 'POI1' PPRIM ;
  319.    NDIM = 'DIME'  LISTINCO ;
  320.    CHPUN = 'MANUEL' 'CHPO' ELT1 1 MOTPRI 1.0 'NATURE' 'DISCRET' ;
  321. ********************** 
  322.   'REPETER'  BL1 NDIM ;
  323.     MOTCEL = 'EXTRAIRE' LISTINCO &BL1 ;
  324.     CHPUN = CHPUN 'ET' ('MANUEL' 'CHPO' MCEN 1 MOTCEL
  325.      0.0 'NATURE' 'DISCRET') ; 
  326.   'FIN' BL1 ;
  327. *********************
  328.    CHCH = 'KOPS'  RMAT 'MULT' CHPUN ;
  329.    SCAL = 'EXTRAIRE' CHCH PDUAL MOTDUA ;
  330.  'FINPROC' SCAL ;
  331.  
  332. ****************************************************************
  333. *****    End of the procedure  JAEXAC                   ********
  334. ****************************************************************
  335. ****************************************************************
  336. ****************   Procedure JACNUM     ************************ 
  337. ****************************************************************
  338.  'DEBPROC' JACNUM $DOMA*'MMODEL' NOMMET*'MOT' RN*'CHPOINT' 
  339.     GN*'CHPOINT' RETN*'CHPOINT' RYN*'CHPOINT' GAMN*'CHPOINT'
  340.     PGAS*'TABLE ' LISTINCO*'LISTMOTS' 
  341.     PPRIM*'POINT' PDUAL*'POINT' 
  342.     MOTPRI*'MOT' MOTDUA*'MOT' EPSILON*'FLOTTANT';
  343.  
  344.  
  345.     RNCEL   = 'COPIER'  RN ;
  346.     GNCEL   = 'COPIER'  GN ;
  347.     RETNCEL = 'COPIER'  RETN ;
  348.     RYNCEL  = 'COPIER'  RYN ;   
  349.  
  350.     VIT PRES TTEMP FRAM GAMM  = 'PRIM' 'PERFMULT'
  351.                         PGAS RNCEL GNCEL RETNCEL RYNCEL  ;
  352.  
  353.     MRN MUN MPN MYN MGAMN = 'PRET' 'PERFMULT' 1 1 $DOMTOT
  354.                     RNCEL VIT PRES FRAM GAMN ;
  355.  
  356.     CHPRESO RFLOT = 'KONV' 'VF' 'PERFMULT' 'RESI' NOMMET
  357.             LISTINCO $DOMA MRN MUN MPN MGAMN MYN PGAS ;
  358.  
  359.     VAL0 = 'EXTRAIRE' CHPRESO PDUAL MOTDUA ;
  360. ******************************************
  361.     dens0 = 'EXTRAIRE' RN PPRIM 'SCAL' ;
  362.  
  363.     CN2 = GAMN '*' (PRES '/' RN) ;
  364.  
  365.     cson0 = ('EXTRAIRE' CN2 PPRIM 'SCAL') '**' 0.5 ;
  366.  
  367.     ret0 = ('EXTRAIRE' RETN PPRIM 'SCAL') '**' 0.5 ;
  368. ******************************************    
  369.     NDIM = 'DIME' LISTINCO ;
  370.     'REPETER' BL1 NDIM ;
  371.        MOTCEL = 'EXTRAIRE' LISTINCO &BL1 ;
  372.        'SI' ('EGA' MOTCEL MOTPRI) ;
  373.           ICEL = &BL1 ;
  374.           'QUITTER' BL1 ;
  375.        'FINSI' ;
  376.     'FIN' BL1 ;
  377.  
  378.     'SI' (ICEL > NDIM) ;
  379.        'MESSAGE' 'Procedure JACNUM' ;
  380.        'MESSAGE' 'MOTPRI = ??? ';
  381.        'ERREUR' 21 ;
  382.     'FINSI' ;
  383.  
  384.     ELT1 = 'MANUEL' 'POI1' PPRIM ;
  385. ****************************************************    
  386. *   ICEL = 1 -> On perturbe la densité
  387. ****************************************************
  388.     'SI' ('EGA' ICEL 1) ;
  389.        DELTATOT = (EPSILON * dens0) ;
  390.        RNCEL = ('MANUEL' 'CHPO' ELT1 1 'SCAL' DELTATOT
  391.           'NATURE' 'DISCRET') 'ET' RN ;
  392.        GNCEL = 'COPIER' GN ;
  393.        RETNCEL = 'COPIER' RETN ;
  394.        RYNCEL = 'COPIER' RYN ;
  395.     'FINSI' ;
  396. ****************************************************
  397. *   ICEL = 2 -> On perturbe la q.d.m. long  l'ax x
  398. ****************************************************
  399.     'SI' ('EGA' ICEL 2) ;
  400.        DELTATOT = (EPSILON * dens0 * cson0) ;
  401.        GNCEL = ('MANUEL' 'CHPO' ELT1 1 'UX' DELTATOT
  402.           'NATURE' 'DISCRET') 'ET' GN ;
  403.        RNCEL = 'COPIER' RN ;
  404.        RETNCEL = 'COPIER' RETN ;
  405.        RYNCEL = 'COPIER' RYN ;
  406.     'FINSI' ;
  407. ****************************************************
  408. *   ICEL = 3 -> On perturbe la q.d.m. long  l'ax y
  409. ****************************************************
  410.     'SI' ('EGA' ICEL 3) ;
  411.        DELTATOT = (EPSILON * dens0 * cson0) ;
  412.        GNCEL = ('MANUEL' 'CHPO' ELT1 1 'UY' DELTATOT
  413.           'NATURE' 'DISCRET') 'ET' GN ;
  414.        RNCEL = 'COPIER' RN ;
  415.        RETNCEL = 'COPIER' RETN ;
  416.        RYNCEL = 'COPIER' RYN ; 
  417.     'FINSI' ;
  418. ****************************************************
  419. *   ICEL = 4 -> On perturbe la q.d.m. long  l'ax y
  420. ****************************************************
  421.     'SI' ('EGA' ICEL 4) ;
  422.        DELTATOT = (EPSILON * dens0 * cson0) ;
  423.        GNCEL = ('MANUEL' 'CHPO' ELT1 1 'UZ' DELTATOT
  424.           'NATURE' 'DISCRET') 'ET' GN ;
  425.        RNCEL = 'COPIER' RN ;
  426.        RETNCEL = 'COPIER' RETN ;
  427.        RYNCEL = 'COPIER' RYN ; 
  428.     'FINSI' ;       
  429. ****************************************************
  430. *   ICEL = 5 -> On perturbe l'énergie totale
  431. ****************************************************
  432.     'SI' ('EGA' ICEL 5) ;
  433.        DELTATOT = (EPSILON * ret0) ;
  434.        RETNCEL = ('MANUEL' 'CHPO' ELT1 1 'SCAL' DELTATOT
  435.           'NATURE' 'DISCRET') 'ET' RETN ;
  436.        RNCEL = 'COPIER' RN ;
  437.        GNCEL = 'COPIER' GN ;
  438.        RYNCEL = 'COPIER' RYN ;
  439.     'FINSI' ;
  440. ****************************************************
  441. *   ICEL = 6 -> On perturbe (\rho Y_{H2})
  442. ****************************************************
  443.     'SI' ('EGA' ICEL 6) ;
  444.        DELTATOT = (EPSILON * ret0) ;
  445.        RYNCEL = ('MANUEL' 'CHPO' ELT1 1 'H2' DELTATOT
  446.           'NATURE' 'DISCRET') 'ET' RYN ;
  447.        RNCEL = 'COPIER' RN ;
  448.        GNCEL = 'COPIER' GN ;
  449.        RETNCEL = 'COPIER' RETN ;
  450.     'FINSI' ;
  451. ****************************************************
  452. *   ICEL = 7 -> On perturbe (\rho Y_{O2})
  453. ****************************************************
  454.     'SI' ('EGA' ICEL 7) ;
  455.        DELTATOT = (EPSILON * ret0) ;
  456.        RYNCEL = ('MANUEL' 'CHPO' ELT1 1 'O2' DELTATOT
  457.           'NATURE' 'DISCRET') 'ET' RYN ;
  458.        RNCEL = 'COPIER' RN ;
  459.        GNCEL = 'COPIER' GN ;
  460.        RETNCEL = 'COPIER' RETN ;
  461.     'FINSI' ;
  462. ****************************************************
  463. *   ICEL = 8 -> On perturbe (\rho Y_{HO})
  464. ****************************************************
  465.     'SI' ('EGA' ICEL 8) ;
  466.        DELTATOT = (EPSILON * ret0) ;
  467.        RYNCEL = ('MANUEL' 'CHPO' ELT1 1 'H2O' DELTATOT
  468.           'NATURE' 'DISCRET') 'ET' RYN ;
  469.        RNCEL = 'COPIER' RN ;
  470.        GNCEL = 'COPIER' GN ;
  471.        RETNCEL = 'COPIER' RETN ;
  472.     'FINSI' ;            
  473. ****************************************************
  474.     VIT PRES TTEMP FRAM GAMM  = 'PRIM' 'PERFMULT'
  475.                         PGAS RNCEL GNCEL RETNCEL RYNCEL  ;
  476.  
  477.     MRNF MUNF MPNF MYNF MGAMNF = 'PRET' 'PERFMULT' 1 1 $DOMTOT
  478.                     RNCEL VIT PRES FRAM GAMM ;
  479.  
  480.     CHPRES1 RFLOT = 'KONV' 'VF' 'PERFMULT' 'RESI' NOMMET
  481.             LISTINCO $DOMA MRNF MUNF MPNF MGAMNF MYNF PGAS ;
  482.  
  483.     VAL1 = 'EXTRAIRE' CHPRES1 PDUAL MOTDUA ;                                  
  484.  
  485.  'FINPROC' ((VAL1 '-' VAL0) '/' DELTATOT) ;
  486. *****************************************************
  487. *****************************************************
  488. ******** FIN PROCEDURES *****************************
  489. *****************************************************
  490. *****************************************************
  491. *****************************************************
  492.  
  493.  
  494. *****************************************************
  495. ******* TEST1 ***************************************
  496. *****************************************************
  497. *
  498. * On compare le jacobien et la variation du residu
  499. * en $DOM1 'CENTRE' par rapport à une variation
  500. * infinitésimal en $DOM1 'CENTRE'
  501. *
  502.  
  503.  PCEN9 = 'POIN' 1 ('DOMA' $DOM1 'CENTRE') ;
  504. ****************************************************************
  505. ****************************************************************
  506.   DRR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMDEN ;
  507.   DGXR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMMOX ;
  508.   DGYR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMMOY ;
  509.   DGZR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMMOZ ;
  510.   DRETR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMRET ;
  511.   DYH2R = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMYH2 ;
  512.   DYO2R = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMYO2 ;
  513.   DYHOR = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMDEN NOMYHO ;
  514. ***************************************
  515.   DRGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMDEN ;
  516.   DGXGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMMOX ;
  517.   DGYGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMMOY ;
  518.   DGZGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMMOZ ;
  519.   DRETGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMRET ;
  520.   DYH2GX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMYH2 ;
  521.   DYO2GX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMYO2 ;
  522.   DYHOGX = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOX NOMYHO ;
  523. ***************************************
  524.   DRGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMDEN ;
  525.   DGXGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMMOX ;
  526.   DGYGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMMOY ;
  527.   DGZGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMMOZ ; 
  528.   DRETGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMRET ;
  529.   DYH2GY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMYH2 ;
  530.   DYO2GY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMYO2 ;
  531.   DYHOGY = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOY NOMYHO ;
  532. ***************************************
  533.   DRGZ = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOZ NOMDEN ;
  534.   DGXGZ = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOZ NOMMOX ;
  535.   DGYGZ = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOZ NOMMOY ;
  536.   DGZGZ = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOZ NOMMOZ ; 
  537.   DRETGZ = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOZ NOMRET ;
  538.   DYH2GZ = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOZ NOMYH2 ;
  539.   DYO2GZ = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOZ NOMYO2 ;
  540.   DYHOGZ = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMMOZ NOMYHO ;    
  541. ***************************************
  542.   DRRET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMDEN ;
  543.   DGXRET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMMOX ;
  544.   DGYRET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMMOY ;
  545.   DGZRET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMMOZ ;
  546.   DRETRET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMRET ;
  547.   DYH2RET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMYH2 ;
  548.   DYO2RET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMYO2 ;
  549.   DYHORET = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMRET NOMYHO ;
  550. ***************************************
  551.   DRYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMDEN ;
  552.   DGXYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMMOX ;
  553.   DGYYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMMOY ;
  554.   DGZYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMMOZ ;
  555.   DRETYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMRET ;
  556.   DYH2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMYH2 ;
  557.   DYO2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMYO2 ;
  558.   DYHOYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYH2 NOMYHO ;
  559. ***************************************
  560.   DRYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMDEN ;
  561.   DGXYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMMOX ;
  562.   DGYYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMMOY ;
  563.   DGZYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMMOZ ;
  564.   DRETYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMRET ;
  565.   DYH2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMYH2 ;
  566.   DYO2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMYO2 ;
  567.   DYHOYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYO2 NOMYHO ;
  568. ****************************************
  569.   DRYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMDEN ;
  570.   DGXYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMMOX ;
  571.   DGYYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMMOY ;
  572.   DGZYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMMOZ ;
  573.   DRETYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMRET ;
  574.   DYH2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMYH2 ;
  575.   DYO2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMYO2 ;
  576.   DYHOYHO = JAEXAC RMAT $DOMTOT LINCO PCEN9 PCEN9 NOMYHO NOMYHO ;
  577.  
  578. *****************************************************************
  579. *****************************************************************
  580. *****************************************************************  
  581.   DRRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  582.       PCEN9 NOMDEN NOMDEN DELTA ;  
  583.   DGXRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  584.       PCEN9 NOMDEN NOMMOX DELTA ;      
  585.   DGYRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  586.       PCEN9 NOMDEN NOMMOY DELTA ;
  587.   DGZRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  588.       PCEN9 NOMDEN NOMMOZ DELTA ;    
  589.   DRETRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  590.       PCEN9 NOMDEN NOMRET DELTA ;
  591.   DYH2RN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  592.       PCEN9 NOMDEN NOMYH2 DELTA ;
  593.   DYO2RN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  594.       PCEN9 NOMDEN NOMYO2 DELTA ;
  595.   DYHORN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  596.       PCEN9 NOMDEN NOMYHO DELTA ;
  597. **************************************
  598.   DRGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  599.       PCEN9 NOMMOX NOMDEN DELTA ;  
  600.   DGXGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  601.       PCEN9 NOMMOX NOMMOX DELTA ;      
  602.   DGYGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  603.       PCEN9 NOMMOX NOMMOY DELTA ;
  604.   DGZGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  605.       PCEN9 NOMMOX NOMMOZ DELTA ;    
  606.   DRETGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  607.       PCEN9 NOMMOX NOMRET DELTA ;
  608.   DYH2GXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  609.       PCEN9 NOMMOX NOMYH2 DELTA ;
  610.   DYO2GXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  611.       PCEN9 NOMMOX NOMYO2 DELTA ;
  612.   DYHOGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  613.       PCEN9 NOMMOX NOMYHO DELTA ;
  614. **************************************
  615.   DRGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  616.       PCEN9 NOMMOY NOMDEN DELTA ;  
  617.   DGXGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  618.       PCEN9 NOMMOY NOMMOX DELTA ;      
  619.   DGYGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  620.       PCEN9 NOMMOY NOMMOY DELTA ;
  621.   DGZGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  622.       PCEN9 NOMMOY NOMMOZ DELTA ;    
  623.   DRETGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  624.       PCEN9 NOMMOY NOMRET DELTA ;
  625.   DYH2GYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  626.       PCEN9 NOMMOY NOMYH2 DELTA ;
  627.   DYO2GYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  628.       PCEN9 NOMMOY NOMYO2 DELTA ;
  629.   DYHOGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  630.       PCEN9 NOMMOY NOMYHO DELTA ;
  631. **************************************
  632.   DRGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  633.       PCEN9 NOMMOZ NOMDEN DELTA ;  
  634.   DGXGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  635.       PCEN9 NOMMOZ NOMMOX DELTA ;      
  636.   DGYGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  637.       PCEN9 NOMMOZ NOMMOY DELTA ;
  638.   DGZGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  639.       PCEN9 NOMMOZ NOMMOZ DELTA ;    
  640.   DRETGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  641.       PCEN9 NOMMOZ NOMRET DELTA ;
  642.   DYH2GZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  643.       PCEN9 NOMMOZ NOMYH2 DELTA ;
  644.   DYO2GZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  645.       PCEN9 NOMMOZ NOMYO2 DELTA ;
  646.   DYHOGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  647.       PCEN9 NOMMOZ NOMYHO DELTA ;      
  648. **************************************
  649.   DRRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  650.       PCEN9 NOMRET NOMDEN DELTA ;  
  651.   DGXRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  652.       PCEN9 NOMRET NOMMOX DELTA ;      
  653.   DGYRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  654.       PCEN9 NOMRET NOMMOY DELTA ;
  655.   DGZRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  656.       PCEN9 NOMRET NOMMOZ DELTA ;    
  657.   DRETRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  658.       PCEN9 NOMRET NOMRET DELTA ;
  659.   DYH2RETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  660.       PCEN9 NOMRET NOMYH2 DELTA ;
  661.   DYO2RETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  662.       PCEN9 NOMRET NOMYO2 DELTA ;
  663.   DYHORETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  664.       PCEN9 NOMRET NOMYHO DELTA ;
  665. **************************************
  666.   DRYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  667.       PCEN9 NOMYH2 NOMDEN DELTA ;  
  668.   DGXYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  669.       PCEN9 NOMYH2 NOMMOX DELTA ;      
  670.   DGYYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  671.       PCEN9 NOMYH2 NOMMOY DELTA ;
  672.   DGZYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  673.       PCEN9 NOMYH2 NOMMOZ DELTA ;    
  674.   DRETYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  675.       PCEN9 NOMYH2 NOMRET DELTA ;
  676.   DYH2YH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  677.       PCEN9 NOMYH2 NOMYH2 DELTA ;
  678.   DYO2YH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  679.       PCEN9 NOMYH2 NOMYO2 DELTA ;
  680.   DYHOYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  681.       PCEN9 NOMYH2 NOMYHO DELTA ;
  682. ***************************************
  683.   DRYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  684.       PCEN9 NOMYO2 NOMDEN DELTA ;  
  685.   DGXYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  686.       PCEN9 NOMYO2 NOMMOX DELTA ;      
  687.   DGYYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  688.       PCEN9 NOMYO2 NOMMOY DELTA ;
  689.   DGZYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  690.       PCEN9 NOMYO2 NOMMOZ DELTA ;    
  691.   DRETYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  692.       PCEN9 NOMYO2 NOMRET DELTA ;
  693.   DYH2YO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  694.       PCEN9 NOMYO2 NOMYH2 DELTA ;
  695.   DYO2YO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  696.       PCEN9 NOMYO2 NOMYO2 DELTA ;
  697.   DYHOYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  698.       PCEN9 NOMYO2 NOMYHO DELTA ;
  699. ***************************************
  700.   DRYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  701.       PCEN9 NOMYHO NOMDEN DELTA ;  
  702.   DGXYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  703.       PCEN9 NOMYHO NOMMOX DELTA ;      
  704.   DGYYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  705.       PCEN9 NOMYHO NOMMOY DELTA ;
  706.   DGZYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  707.       PCEN9 NOMYHO NOMMOZ DELTA ;    
  708.   DRETYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  709.       PCEN9 NOMYHO NOMRET DELTA ;
  710.   DYH2YHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  711.       PCEN9 NOMYHO NOMYH2 DELTA ;
  712.   DYO2YHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  713.       PCEN9 NOMYHO NOMYO2 DELTA ;
  714.   DYHOYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN9
  715.       PCEN9 NOMYHO NOMYHO DELTA ;  
  716.  
  717. * Test des comparaisons jacobien exact-jacobien numerique
  718. * Test des comparaisons jacobien exact-jacobien numerique
  719. **************************************************************
  720.  'SI' (('ABS'(DRR)) < DELTA) ;
  721.     ERR1 = ('ABS' (DRRN '-' DRR)) ;
  722.  'SINON' ;  
  723.     ERR1 = ('ABS' (DRRN '-' DRR)) '/' ('ABS' (DRR)) ;
  724.  'FINSI' ; 
  725. ***********************
  726.  'SI' (ERR1 > ERRTOL) ;
  727.       'MESSAGE' 'Problem 1 T1';
  728.       'ERREUR' 5 ;
  729.  'FINSI' ;
  730. **************************************************************
  731.  'SI' (('ABS'(DGXR)) < DELTA) ;
  732.     ERR1 = ('ABS' (DGXRN '-' DGXR)) ;
  733.  'SINON' ;  
  734.     ERR1 = ('ABS' (DGXRN '-' DGXR)) '/' ('ABS' (DGXR)) ;
  735.  'FINSI' ; 
  736. *********************** 
  737.  'SI' (ERR1 > ERRTOL) ;
  738.       'MESSAGE' 'Problem 2 T1';
  739.       'ERREUR' 5 ;
  740.  'FINSI' ;
  741. **************************************************************
  742.  'SI' (('ABS'(DGYR)) < DELTA) ;
  743.     ERR1 = ('ABS' (DGYRN '-' DGYR)) ;
  744.  'SINON' ;  
  745.     ERR1 = ('ABS' (DGYRN '-' DGYR)) '/' ('ABS' (DGYR)) ;
  746.  'FINSI' ; 
  747. ************************
  748.  'SI' (ERR1 > ERRTOL) ;
  749.       'MESSAGE' 'Problem 3 T1';
  750.        'ERREUR' 5 ;
  751.  'FINSI' ;
  752. **************************************************************
  753.  'SI' (('ABS'(DGZR)) < DELTA) ;
  754.     ERR1 = ('ABS' (DGZRN '-' DGZR)) ;
  755.  'SINON' ;  
  756.     ERR1 = ('ABS' (DGZRN '-' DGZR)) '/' ('ABS' (DGZR)) ;
  757.  'FINSI' ;
  758. ************************
  759.  'SI' (ERR1 > ERRTOL) ;
  760.       'MESSAGE' 'Problem 4 T1';
  761.        'ERREUR' 5 ;
  762.  'FINSI' ;
  763. ***************************************************************
  764.  'SI' (('ABS'(DRETR)) < DELTA) ;
  765.     ERR1 = ('ABS' (DRETRN '-' DRETR)) ;
  766.  'SINON' ;  
  767.     ERR1 = ('ABS' (DRETRN '-' DRETR)) '/' ('ABS' (DRETR)) ;
  768.  'FINSI' ;
  769. **************************
  770.  'SI' (ERR1 > ERRTOL) ;
  771.       'MESSAGE' 'Problem 5 T1';
  772.       'ERREUR' 5 ;
  773.  'FINSI' ;
  774. **************************************************************
  775. **************************************************************
  776.  'SI' (('ABS'(DYH2R)) < DELTA) ;
  777.     ERR1 = ('ABS' (DYH2RN '-' DYH2R)) ;
  778.  'SINON' ;  
  779.     ERR1 = ('ABS' (DYH2RN '-' DYH2R)) '/' ('ABS' (DYH2R)) ;
  780.  'FINSI' ;
  781. **************************
  782.  'SI' (ERR1 > ERRTOL) ;
  783.       'MESSAGE' 'Problem 6 T1';
  784.       'ERREUR' 5 ;
  785.  'FINSI' ;
  786. **************************************************************
  787.    'SI' (('ABS'(DYO2R)) < DELTA) ;
  788.     ERR1 = ('ABS' (DYO2RN '-' DYO2R)) ;
  789.  'SINON' ;  
  790.     ERR1 = ('ABS' (DYO2RN '-' DYO2R)) '/' ('ABS' (DYO2R)) ;
  791.  'FINSI' ;
  792. **************************
  793.  'SI' (ERR1 > ERRTOL) ;
  794.       'MESSAGE' 'Problem 7 T1';
  795.       'ERREUR' 5 ;
  796.  'FINSI' ;
  797. **************************************************************
  798.   'SI' (('ABS'(DYHOR)) < DELTA) ;
  799.     ERR1 = ('ABS' (DYHORN '-' DYHOR)) ;
  800.  'SINON' ;  
  801.     ERR1 = ('ABS' (DYHORN '-' DYHOR)) '/' ('ABS' (DYHOR)) ;
  802.  'FINSI' ;
  803. **************************
  804.  'SI' (ERR1 > ERRTOL) ;
  805.       'MESSAGE' 'Problem 8 T1';
  806.       'ERREUR' 5 ;
  807.  'FINSI' ; 
  808. **************************************************************
  809. **************************************************************
  810. **************************************************************
  811. *****    Second line !!!!!!!!!!!!!!!!!!!!!!!    **************
  812. **************************************************************
  813. ************************************************************** 
  814.  'SI' (('ABS'(DRGX)) < DELTA) ;
  815.     ERR1 = ('ABS' (DRGXN '-' DRGX)) ;
  816.  'SINON' ;   
  817.     ERR1 = ('ABS' (DRGXN '-' DRGX)) '/' ('ABS' (DRGX))  ;
  818.  'FINSI' ;
  819. *************************** 
  820.  'SI' (ERR1 > ERRTOL) ;
  821.       'MESSAGE' 'Problem 9 T1';
  822.       'ERREUR' 5 ;
  823.  'FINSI' ;
  824. ***************************************************************
  825.  'SI' (('ABS'(DGXGX)) < DELTA) ;
  826.    ERR1 = ('ABS' (DGXGXN '-' DGXGX)) ;
  827.  'SINON' ;
  828.    ERR1 = ('ABS' (DGXGXN '-' DGXGX)) '/' ('ABS' (DGXGX)) ;
  829.  'FINSI' ;  
  830. *****************************   
  831.  'SI' (ERR1 > ERRTOL) ;
  832.         'MESSAGE' 'Problem 10 T1';
  833.         'ERREUR' 5 ;
  834.  'FINSI' ;
  835. **************************************************************
  836.  'SI' (('ABS'(DGYGX)) < DELTA) ;
  837.     ERR1 = ('ABS' (DGYGXN '-' DGYGX)) ;
  838.  'SINON' ;
  839.     ERR1 = ('ABS' (DGYGXN '-' DGYGX)) '/' ('ABS' (DGYGX)) ;
  840.  'FINSI' ;
  841. ******************************     
  842.  'SI' (ERR1 > ERRTOL) ;
  843.         'MESSAGE' 'Problem 11 T1';
  844.         'ERREUR' 5 ;
  845.  'FINSI' ;
  846. **************************************************************
  847.  'SI' (('ABS'(DGZGX)) < DELTA) ;
  848.     ERR1 = ('ABS' (DGZGXN '-' DGZGX)) ;
  849.  'SINON' ;
  850.     ERR1 = ('ABS' (DGZGXN '-' DGZGX)) '/' ('ABS' (DGZGX)) ;
  851.  'FINSI' ;
  852. ******************************    
  853.  'SI' (ERR1 > ERRTOL) ;
  854.         'MESSAGE' 'Problem 12 T1';
  855.         'ERREUR' 5 ;
  856.  'FINSI' ;
  857. **************************************************************
  858.  'SI' (('ABS'(DRETGX)) < DELTA) ;
  859.     ERR1 = ('ABS' (DRETGXN '-' DRETGX)) ;
  860.  'SINON' ;
  861.     ERR1 = ('ABS' (DRETGXN '-' DRETGX)) '/' ('ABS' (DRETGX)) ;
  862.  'FINSI' ;
  863. *******************************    
  864.  'SI' (ERR1 > ERRTOL) ;
  865.         'MESSAGE' 'Problem 13 T1';
  866.         'ERREUR' 5 ;
  867.  'FINSI' ;
  868. **************************************************************
  869.  'SI' (('ABS'(DYH2GX)) < DELTA) ;
  870.     ERR1 = ('ABS' (DYH2GXN '-' DYH2GX)) ;
  871.  'SINON' ;  
  872.     ERR1 = ('ABS' (DYH2GXN '-' DYH2GX)) '/' ('ABS' (DYH2GX)) ;
  873.  'FINSI' ;
  874. **************************
  875.  'SI' (ERR1 > ERRTOL) ;
  876.       'MESSAGE' 'Problem 14 T1';
  877.       'ERREUR' 5 ;
  878.  'FINSI' ;
  879. **************************************************************
  880.  'SI' (('ABS'(DYO2GX)) < DELTA) ;
  881.     ERR1 = ('ABS' (DYO2GXN '-' DYO2GX)) ;
  882.  'SINON' ;  
  883.     ERR1 = ('ABS' (DYO2GXN '-' DYO2GX)) '/' ('ABS' (DYO2GX)) ;
  884.  'FINSI' ;
  885. **************************
  886.  'SI' (ERR1 > ERRTOL) ;
  887.       'MESSAGE' 'Problem 15 T1';
  888.       'ERREUR' 5 ;
  889.  'FINSI' ;
  890. **************************************************************
  891.  'SI' (('ABS'(DYHOGX)) < DELTA) ;
  892.     ERR1 = ('ABS' (DYHOGXN '-' DYHOGX)) ;
  893.  'SINON' ;  
  894.     ERR1 = ('ABS' (DYHOGXN '-' DYHOGX)) '/' ('ABS' (DYHOGX)) ;
  895.  'FINSI' ;
  896. **************************
  897.  'SI' (ERR1 > ERRTOL) ;
  898.       'MESSAGE' 'Problem 16 T1';
  899.       'ERREUR' 5 ;
  900.  'FINSI' ;
  901. **************************************************************
  902. *******        Third line!!!!!!!!!!!!!!!!      ***************
  903. **************************************************************
  904.  'SI' (('ABS'(DRGY)) < DELTA) ;
  905.     ERR1 = ('ABS' (DRGYN '-' DRGY)) ;
  906.  'SINON' ;
  907.     ERR1 = ('ABS' (DRGYN '-' DRGY)) '/' ('ABS' (DRGY)) ;
  908.  'FINSI' ;
  909. *******************************       
  910.  'SI' (ERR1 > ERRTOL) ;
  911.         'MESSAGE' 'Problem 17 T1';
  912.         'ERREUR' 5 ;
  913.  'FINSI' ;
  914. **************************************************************
  915.  'SI' (('ABS'(DGXGY)) < DELTA) ;
  916.     ERR1 = ('ABS' (DGXGYN '-' DGXGY)) ;
  917.  'SINON' ;
  918.     ERR1 = ('ABS' (DGXGYN '-' DGXGY))  '/' ('ABS' (DGXGY)) ;
  919.  'FINSI' ;
  920. ********************************      
  921.  'SI' (ERR1 > ERRTOL) ;
  922.         'MESSAGE' 'Problem 18 T1';
  923.         'ERREUR' 5 ;
  924.  'FINSI' ;
  925. **************************************************************
  926.  'SI' (('ABS'(DGYGY)) < DELTA) ;
  927.     ERR1 = ('ABS' (DGYGYN '-' DGYGY)) ;
  928.  'SINON' ;
  929.     ERR1 = ('ABS' (DGYGYN '-' DGYGY)) '/' ('ABS' (DGYGY)) ;
  930.  'FINSI' ;
  931. ***************
  932.  'SI' (ERR1 > ERRTOL) ;
  933.         'MESSAGE' 'Problem 19 T1';
  934.         'ERREUR' 5 ;
  935.  'FINSI' ;
  936. **************************************************************
  937.  'SI' (('ABS'(DGZGY)) < DELTA) ;
  938.     ERR1 = ('ABS' (DGZGYN '-' DGZGY)) ;
  939.  'SINON' ;
  940.     ERR1 = ('ABS' (DGZGYN '-' DGZGY)) '/' ('ABS' (DGZGY)) ;
  941.  'FINSI' ;
  942. ****************        
  943.  'SI' (ERR1 > ERRTOL) ;
  944.         'MESSAGE' 'Problem 20 T1';
  945.         'ERREUR' 5 ;
  946.  'FINSI' ;
  947. **************************************************************
  948.  'SI' (('ABS'(DRETGY)) < DELTA) ;
  949.     ERR1 = ('ABS' (DRETGYN '-' DRETGY)) ;
  950.  'SINON' ;
  951.     ERR1 = ('ABS' (DRETGYN '-' DRETGY)) '/' ('ABS' (DRETGY)) ;
  952.  'FINSI' ;
  953. ****************     
  954.  'SI' (ERR1 > ERRTOL) ;
  955.         'MESSAGE' 'Problem 21 T1';
  956.         'ERREUR' 5 ;
  957.  'FINSI' ;
  958. **************************************************************
  959.  'SI' (('ABS'(DYH2GY)) < DELTA) ;
  960.     ERR1 = ('ABS' (DYH2GYN '-' DYH2GY)) ;
  961.  'SINON' ;  
  962.     ERR1 = ('ABS' (DYH2GYN '-' DYH2GY)) '/' ('ABS' (DYH2GY)) ;
  963.  'FINSI' ;
  964. **************************
  965.  'SI' (ERR1 > ERRTOL) ;
  966.       'MESSAGE' 'Problem 22 T1';
  967.       'ERREUR' 5 ;
  968.  'FINSI' ;
  969. **************************************************************
  970.  'SI' (('ABS'(DYO2GY)) < DELTA) ;
  971.     ERR1 = ('ABS' (DYO2GYN '-' DYO2GY)) ;
  972.  'SINON' ;  
  973.     ERR1 = ('ABS' (DYO2GYN '-' DYO2GY)) '/' ('ABS' (DYO2GY)) ;
  974.  'FINSI' ;
  975. **************************
  976.  'SI' (ERR1 > ERRTOL) ;
  977.       'MESSAGE' 'Problem 23 T1';
  978.       'ERREUR' 5 ;
  979.  'FINSI' ;
  980. **************************************************************
  981.  'SI' (('ABS'(DYHOGY)) < DELTA) ;
  982.     ERR1 = ('ABS' (DYHOGYN '-' DYHOGY)) ;
  983.  'SINON' ;  
  984.     ERR1 = ('ABS' (DYHOGYN '-' DYHOGY)) '/' ('ABS' (DYHOGY)) ;
  985.  'FINSI' ;
  986. **************************
  987.  'SI' (ERR1 > ERRTOL) ;
  988.       'MESSAGE' 'Problem 24 T1';
  989.       'ERREUR' 5 ;
  990.  'FINSI' ;
  991. **************************************************************
  992. *******     Fouth line!!!!!!!!!!!!!!!!!!!   ******************
  993. **************************************************************
  994.  'SI' (('ABS'(DRGZ)) < DELTA) ;
  995.     ERR1 = ('ABS' (DRGZN '-' DRGZ))  ;
  996.  'SINON' ;
  997.     ERR1 = ('ABS' (DRGZN '-' DRGZ)) '/' ('ABS' (DRGZ))  ;
  998.  'FINSI' ;
  999. ***************      
  1000.  'SI' (ERR1 > ERRTOL) ;
  1001.         'MESSAGE' 'Problem 25 T1';
  1002.         'ERREUR' 5 ;
  1003.  'FINSI' ;
  1004. **************************************************************
  1005.  'SI' (('ABS'(DGXGZ)) < DELTA) ;
  1006.     ERR1 = ('ABS' (DGXGZN '-' DGXGZ)) ;
  1007.  'SINON' ;
  1008.     ERR1 = ('ABS' (DGXGZN '-' DGXGZ))  '/' ('ABS' (DGXGZ)) ;
  1009.  'FINSI' ;
  1010. *****************       
  1011.  'SI' (ERR1 > ERRTOL) ;
  1012.         'MESSAGE' 'Problem 26 T1';
  1013.         'ERREUR' 5 ;
  1014.  'FINSI' ;
  1015. **************************************************************
  1016.  'SI' (('ABS'(DGYGZ)) < DELTA) ;
  1017.     ERR1 = ('ABS' (DGYGZN '-' DGYGZ)) ;
  1018.  'SINON' ;
  1019.     ERR1 = ('ABS' (DGYGZN '-' DGYGZ)) '/' ('ABS' (DGYGZ)) ;
  1020.  'FINSI' ;
  1021. *****************       
  1022.  'SI' (ERR1 > ERRTOL) ;
  1023.         'MESSAGE' 'Problem 27 T1';
  1024.         'ERREUR' 5 ;
  1025.  'FINSI' ;
  1026. **************************************************************
  1027.  'SI' (('ABS'(DGZGZ)) < DELTA) ;
  1028.     ERR1 = ('ABS' (DGZGZN '-' DGZGZ)) ;
  1029.  'SINON' ;
  1030.     ERR1 = ('ABS' (DGZGZN '-' DGZGZ)) '/' ('ABS' (DGZGZ)) ;
  1031.  'FINSI' ;
  1032. *****************     
  1033.  'SI' (ERR1 > ERRTOL) ;
  1034.         'MESSAGE' 'Problem 28 T1';
  1035.         'ERREUR' 5 ;
  1036.  'FINSI' ;
  1037. **************************************************************
  1038.  'SI' (('ABS'(DRETGZ)) < DELTA) ;
  1039.     ERR1 = ('ABS' (DRETGZN '-' DRETGZ)) ;
  1040.  'SINON' ;
  1041.     ERR1 = ('ABS' (DRETGZN '-' DRETGZ)) '/' ('ABS' (DRETGZ)) ;
  1042.  'FINSI' ;     
  1043. *****************
  1044.  'SI' (ERR1 > ERRTOL) ;
  1045.         'MESSAGE' 'Problem 29 T1';
  1046.         'ERREUR' 5 ;
  1047.  'FINSI' ;
  1048. **************************************************************
  1049.  'SI' (('ABS'(DYH2GZ)) < DELTA) ;
  1050.     ERR1 = ('ABS' (DYH2GZN '-' DYH2GZ)) ;
  1051.  'SINON' ;  
  1052.     ERR1 = ('ABS' (DYH2GZN '-' DYH2GZ)) '/' ('ABS' (DYH2GZ)) ;
  1053.  'FINSI' ;
  1054. **************************
  1055.  'SI' (ERR1 > ERRTOL) ;
  1056.       'MESSAGE' 'Problem 30 T1';
  1057.       'ERREUR' 5 ;
  1058.  'FINSI' ;
  1059. **************************************************************
  1060.  'SI' (('ABS'(DYO2GZ)) < DELTA) ;
  1061.     ERR1 = ('ABS' (DYO2GZN '-' DYO2GZ)) ;
  1062.  'SINON' ;  
  1063.     ERR1 = ('ABS' (DYO2GZN '-' DYO2GZ)) '/' ('ABS' (DYO2GZ)) ;
  1064.  'FINSI' ;
  1065. **************************
  1066.  'SI' (ERR1 > ERRTOL) ;
  1067.       'MESSAGE' 'Problem 31 T1';
  1068.       'ERREUR' 5 ;
  1069.  'FINSI' ;
  1070. **************************************************************
  1071.  'SI' (('ABS'(DYHOGZ)) < DELTA) ;
  1072.     ERR1 = ('ABS' (DYHOGZN '-' DYHOGZ)) ;
  1073.  'SINON' ;  
  1074.     ERR1 = ('ABS' (DYHOGZN '-' DYHOGZ)) '/' ('ABS' (DYHOGZ)) ;
  1075.  'FINSI' ;
  1076. **************************
  1077.  'SI' (ERR1 > ERRTOL) ;
  1078.       'MESSAGE' 'Problem 32 T1';
  1079.       'ERREUR' 5 ;
  1080.  'FINSI' ;
  1081. **************************************************************
  1082. *******     Fith line!!!!!!!!!!!!!!!!!!!   ******************
  1083. **************************************************************
  1084.  'SI' (('ABS'(DRRET)) < DELTA) ;
  1085.     ERR1 = ('ABS' (DRRETN '-' DRRET)) ;
  1086.  'SINON' ;
  1087.     ERR1 = ('ABS' (DRRETN '-' DRRET)) '/' ('ABS' (DRRET)) ;
  1088.  'FINSI' ;
  1089. *****************       
  1090.  'SI' (ERR1 > ERRTOL) ;
  1091.         'MESSAGE' 'Problem 33 T1';
  1092.         'ERREUR' 5 ;
  1093.  'FINSI' ;
  1094. ***************************************************************
  1095.  'SI' (('ABS'(DGXRET)) < DELTA) ;
  1096.     ERR1 = ('ABS' (DGXRETN '-' DGXRET)) ;
  1097.  'SINON' ;
  1098.     ERR1 = ('ABS' (DGXRETN '-' DGXRET)) '/' ('ABS' (DGXRET)) ;
  1099.  'FINSI' ;
  1100. *****************    
  1101.  'SI' (ERR1 > ERRTOL) ;
  1102.         'MESSAGE' 'Problem 34 T1';
  1103.         'ERREUR' 5 ;
  1104.  'FINSI' ;
  1105. ****************************************************************
  1106.  'SI' (('ABS'(DGYRET)) < DELTA) ;
  1107.     ERR1 = ('ABS' (DGYRETN '-' DGYRET)) ; 
  1108.  'SINON' ;
  1109.     ERR1 = ('ABS' (DGYRETN '-' DGYRET)) '/' ('ABS' (DGYRET)) ;
  1110.  'FINSI' ;
  1111. *****************    
  1112.  'SI' (ERR1 > ERRTOL) ;
  1113.         'MESSAGE' 'Problem 35 T1';
  1114.         'ERREUR' 5 ;
  1115.  'FINSI' ;
  1116. ****************************************************************
  1117.  'SI' (('ABS'(DGZRET)) < DELTA) ;
  1118.     ERR1 = ('ABS' (DGZRETN '-' DGZRET)) ;
  1119.  'SINON' ;    
  1120.     ERR1 = ('ABS' (DGZRETN '-' DGZRET)) '/' ('ABS' (DGZRET)) ;
  1121.  'FINSI' ;
  1122. *****************    
  1123.  'SI' (ERR1 > ERRTOL) ;
  1124.         'MESSAGE' 'Problem 36 T1';
  1125.         'ERREUR' 5 ;
  1126.  'FINSI' ;
  1127. ****************************************************************
  1128.  'SI' (('ABS'(DRETRET)) < DELTA) ;
  1129.     ERR1 = ('ABS' (DRETRETN '-' DRETRET)) ;
  1130.  'SINON' ;
  1131.     ERR1 = ('ABS' (DRETRETN '-' DRETRET)) '/' ('ABS' (DRETRET)) ;
  1132.  'FINSI' ;
  1133. *****************       
  1134.  'SI' (ERR1 > ERRTOL) ;
  1135.         'MESSAGE' 'Problem 37 T1';
  1136.         'ERREUR' 5 ;
  1137.  'FINSI' ;
  1138.  'SI' (('ABS'(DYH2RET)) < DELTA) ;
  1139.     ERR1 = ('ABS' (DYH2RETN '-' DYH2RET)) ;
  1140.  'SINON' ;  
  1141.     ERR1 = ('ABS' (DYH2RETN '-' DYH2RET)) '/' ('ABS' (DYH2RET)) ;
  1142.  'FINSI' ;
  1143. **************************
  1144.  'SI' (ERR1 > ERRTOL) ;
  1145.       'MESSAGE' 'Problem 38 T1';
  1146.       'ERREUR' 5 ;
  1147.  'FINSI' ;
  1148. **************************************************************
  1149.  'SI' (('ABS'(DYO2RET)) < DELTA) ;
  1150.     ERR1 = ('ABS' (DYO2RETN '-' DYO2RET)) ;
  1151.  'SINON' ;  
  1152.     ERR1 = ('ABS' (DYO2RETN '-' DYO2RET)) '/' ('ABS' (DYO2RET)) ;
  1153.  'FINSI' ;
  1154. **************************
  1155.  'SI' (ERR1 > ERRTOL) ;
  1156.       'MESSAGE' 'Problem 39 T1';
  1157.       'ERREUR' 5 ;
  1158.  'FINSI' ;
  1159. **************************************************************
  1160.  'SI' (('ABS'(DYHORET)) < DELTA) ;
  1161.     ERR1 = ('ABS' (DYHORETN '-' DYHORET)) ;
  1162.  'SINON' ;  
  1163.     ERR1 = ('ABS' (DYHORETN '-' DYHORET)) '/' ('ABS' (DYHORET)) ;
  1164.  'FINSI' ;
  1165. **************************
  1166.  'SI' (ERR1 > ERRTOL) ;
  1167.       'MESSAGE' 'Problem 40 T1';
  1168.       'ERREUR' 5 ;
  1169.  'FINSI' ;
  1170. **************************************************************
  1171. *******     Sixth line!!!!!!!!!!!!!!!!!!!   ******************
  1172. **************************************************************
  1173.  'SI' (('ABS'(DRYH2)) < DELTA) ;
  1174.     ERR1 = ('ABS' (DRYH2N '-' DRYH2)) ;
  1175.  'SINON' ;
  1176.     ERR1 = ('ABS' (DRYH2N '-' DRYH2)) '/' ('ABS' (DRYH2)) ;
  1177.  'FINSI' ;
  1178. *****************
  1179.  'SI' (ERR1 > ERRTOL) ;
  1180.         'MESSAGE' 'Problem 41 T1';
  1181.         'ERREUR' 5 ;
  1182.  'FINSI' ;
  1183. **************************************************************
  1184.  'SI' (('ABS'(DGXYH2)) < DELTA) ;
  1185.     ERR1 = ('ABS' (DGXYH2N '-' DGXYH2)) ;
  1186.  'SINON' ;
  1187.     ERR1 = ('ABS' (DGXYH2N '-' DGXYH2)) '/' ('ABS' (DGXYH2)) ;
  1188.  'FINSI' ;
  1189. *****************
  1190.  'SI' (ERR1 > ERRTOL) ;
  1191.         'MESSAGE' 'Problem 42 T1';
  1192.         'ERREUR' 5 ;
  1193.  'FINSI' ;
  1194. *************************************************************
  1195.  'SI' (('ABS'(DGYYH2)) < DELTA) ;
  1196.     ERR1 = ('ABS' (DGYYH2N '-' DGYYH2)) ; 
  1197.  'SINON' ;
  1198.     ERR1 = ('ABS' (DGYYH2N '-' DGYYH2)) '/' ('ABS' (DGYYH2)) ;
  1199.  'FINSI' ;
  1200. *****************    
  1201.  'SI' (ERR1 > ERRTOL) ;
  1202.         'MESSAGE' 'Problem 43 T1';
  1203.         'ERREUR' 5 ;
  1204.  'FINSI' ;
  1205. **************************************************************
  1206.  'SI' (('ABS'(DGZYH2)) < DELTA) ;
  1207.     ERR1 = ('ABS' (DGZYH2N '-' DGZYH2)) ; 
  1208.  'SINON' ;
  1209.     ERR1 = ('ABS' (DGZYH2N '-' DGZYH2)) '/' ('ABS' (DGZYH2)) ;
  1210.  'FINSI' ;
  1211. *****************    
  1212.  'SI' (ERR1 > ERRTOL) ;
  1213.         'MESSAGE' 'Problem 44 T1';
  1214.         'ERREUR' 5 ;
  1215.  'FINSI' ; 
  1216. **************************************************************
  1217.  'SI' (('ABS'(DRETYH2)) < DELTA) ;
  1218.     ERR1 = ('ABS' (DRETYH2N '-' DRETYH2)) ;
  1219.  'SINON' ;
  1220.     ERR1 = ('ABS' (DRETYH2N '-' DRETYH2)) '/' ('ABS' (DRETYH2)) ;
  1221.  'FINSI' ;
  1222. *****************
  1223.  'SI' (ERR1 > ERRTOL) ;
  1224.         'MESSAGE' 'Problem 45 T1';
  1225.         'ERREUR' 5 ;
  1226.  'FINSI' ;
  1227. **************************************************************
  1228.  'SI' (('ABS'(DYH2YH2)) < DELTA) ;
  1229.     ERR1 = ('ABS' (DYH2YH2N '-' DYH2YH2)) ;
  1230.  'SINON' ;  
  1231.     ERR1 = ('ABS' (DYH2YH2N '-' DYH2YH2)) '/' ('ABS' (DYH2YH2)) ;
  1232.  'FINSI' ;
  1233. **************************
  1234.  'SI' (ERR1 > ERRTOL) ;
  1235.       'MESSAGE' 'Problem 46 T1';
  1236.       'ERREUR' 5 ;
  1237.  'FINSI' ;
  1238. **************************************************************
  1239.  'SI' (('ABS'(DYO2YH2)) < DELTA) ;
  1240.     ERR1 = ('ABS' (DYO2YH2N '-' DYO2YH2)) ;
  1241.  'SINON' ;  
  1242.     ERR1 = ('ABS' (DYO2YH2N '-' DYO2YH2)) '/' ('ABS' (DYO2YH2)) ;
  1243.  'FINSI' ;
  1244. **************************
  1245.  'SI' (ERR1 > ERRTOL) ;
  1246.       'MESSAGE' 'Problem 47 T1';
  1247.       'ERREUR' 5 ;
  1248.  'FINSI' ;
  1249. **************************************************************
  1250.  'SI' (('ABS'(DYHOYH2)) < DELTA) ;
  1251.     ERR1 = ('ABS' (DYHOYH2N '-' DYHOYH2)) ;
  1252.  'SINON' ;  
  1253.     ERR1 = ('ABS' (DYHOYH2N '-' DYHOYH2)) '/' ('ABS' (DYHOYH2)) ;
  1254.  'FINSI' ;
  1255. **************************
  1256.  'SI' (ERR1 > ERRTOL) ;
  1257.       'MESSAGE' 'Problem 48 T1';
  1258.       'ERREUR' 5 ;
  1259.  'FINSI' ;
  1260. **************************************************************
  1261. *********  Seventh line !!!!!!!!!!!!!!!!!!!!!!   *************
  1262. **************************************************************
  1263.  'SI' (('ABS'(DRYO2)) < DELTA) ;
  1264.     ERR1 = ('ABS' (DRYO2N '-' DRYO2)) ;
  1265.  'SINON' ;
  1266.     ERR1 = ('ABS' (DRYO2N '-' DRYO2)) '/' ('ABS' (DRYO2)) ;
  1267.  'FINSI' ;
  1268. *****************
  1269.  'SI' (ERR1 > ERRTOL) ;
  1270.         'MESSAGE' 'Problem 49 T1';
  1271.         'ERREUR' 5 ;
  1272.  'FINSI' ;
  1273. **************************************************************
  1274.  'SI' (('ABS'(DGXYO2)) < DELTA) ;
  1275.     ERR1 = ('ABS' (DGXYO2N '-' DGXYO2)) ;
  1276.  'SINON' ;
  1277.     ERR1 = ('ABS' (DGXYO2N '-' DGXYO2)) '/' ('ABS' (DGXYO2)) ;
  1278.  'FINSI' ;
  1279. *****************
  1280.  'SI' (ERR1 > ERRTOL) ;
  1281.         'MESSAGE' 'Problem 50 T1';
  1282.         'ERREUR' 5 ;
  1283.  'FINSI' ;
  1284. *************************************************************
  1285.  'SI' (('ABS'(DGYYO2)) < DELTA) ;
  1286.     ERR1 = ('ABS' (DGYYO2N '-' DGYYO2)) ; 
  1287.  'SINON' ;
  1288.     ERR1 = ('ABS' (DGYYO2N '-' DGYYO2)) '/' ('ABS' (DGYYO2)) ;
  1289.  'FINSI' ;
  1290. *****************    
  1291.  'SI' (ERR1 > ERRTOL) ;
  1292.         'MESSAGE' 'Problem 51 T1';
  1293.         'ERREUR' 5 ;
  1294.  'FINSI' ;
  1295. **************************************************************
  1296.  'SI' (('ABS'(DGZYO2)) < DELTA) ;
  1297.     ERR1 = ('ABS' (DGZYO2N '-' DGZYO2)) ; 
  1298.  'SINON' ;
  1299.     ERR1 = ('ABS' (DGZYO2N '-' DGZYO2)) '/' ('ABS' (DGZYO2)) ;
  1300.  'FINSI' ;
  1301. *****************    
  1302.  'SI' (ERR1 > ERRTOL) ;
  1303.         'MESSAGE' 'Problem 52 T1';
  1304.         'ERREUR' 5 ;
  1305.  'FINSI' ; 
  1306. **************************************************************
  1307.  'SI' (('ABS'(DRETYO2)) < DELTA) ;
  1308.     ERR1 = ('ABS' (DRETYO2N '-' DRETYO2)) ;
  1309.  'SINON' ;
  1310.     ERR1 = ('ABS' (DRETYO2N '-' DRETYO2)) '/' ('ABS' (DRETYO2)) ;
  1311.  'FINSI' ;
  1312. *****************
  1313.  'SI' (ERR1 > ERRTOL) ;
  1314.         'MESSAGE' 'Problem  53 T1';
  1315.         'ERREUR' 5 ;
  1316.  'FINSI' ;
  1317. **************************************************************
  1318.  'SI' (('ABS'(DYH2YO2)) < DELTA) ;
  1319.     ERR1 = ('ABS' (DYH2YO2N '-' DYH2YO2)) ;
  1320.  'SINON' ;  
  1321.     ERR1 = ('ABS' (DYH2YO2N '-' DYH2YO2)) '/' ('ABS' (DYH2YO2)) ;
  1322.  'FINSI' ;
  1323. **************************
  1324.  'SI' (ERR1 > ERRTOL) ;
  1325.       'MESSAGE' 'Problem 54 T1';
  1326.       'ERREUR' 5 ;
  1327.  'FINSI' ;
  1328. **************************************************************
  1329.  'SI' (('ABS'(DYO2YO2)) < DELTA) ;
  1330.     ERR1 = ('ABS' (DYO2YO2N '-' DYO2YO2)) ;
  1331.  'SINON' ;  
  1332.     ERR1 = ('ABS' (DYO2YO2N '-' DYO2YO2)) '/' ('ABS' (DYO2YO2)) ;
  1333.  'FINSI' ;
  1334. **************************
  1335.  'SI' (ERR1 > ERRTOL) ;
  1336.       'MESSAGE' 'Problem 55 T1';
  1337.       'ERREUR' 5 ;
  1338.  'FINSI' ;
  1339. **************************************************************
  1340.  'SI' (('ABS'(DYHOYO2)) < DELTA) ;
  1341.     ERR1 = ('ABS' (DYHOYO2N '-' DYHOYO2)) ;
  1342.  'SINON' ;  
  1343.     ERR1 = ('ABS' (DYHOYO2N '-' DYHOYO2)) '/' ('ABS' (DYHOYO2)) ;
  1344.  'FINSI' ;
  1345. **************************
  1346.  'SI' (ERR1 > ERRTOL) ;
  1347.       'MESSAGE' 'Problem 56 T1';
  1348.       'ERREUR' 5 ;
  1349.  'FINSI' ;
  1350. **************************************************************
  1351. *********  Eigth line !!!!!!!!!!!!!!!!!!!!!!   ***************
  1352. **************************************************************
  1353.  'SI' (('ABS'(DRYHO)) < DELTA) ;
  1354.     ERR1 = ('ABS' (DRYHON '-' DRYHO)) ;
  1355.  'SINON' ;
  1356.     ERR1 = ('ABS' (DRYHON '-' DRYHO)) '/' ('ABS' (DRYHO)) ;
  1357.  'FINSI' ;
  1358. *****************
  1359.  'SI' (ERR1 > ERRTOL) ;
  1360.         'MESSAGE' 'Problem 57 T1';
  1361.         'ERREUR' 5 ;
  1362.  'FINSI' ;
  1363. **************************************************************
  1364.  'SI' (('ABS'(DGXYHO)) < DELTA) ;
  1365.     ERR1 = ('ABS' (DGXYHON '-' DGXYHO)) ;
  1366.  'SINON' ;
  1367.     ERR1 = ('ABS' (DGXYHON '-' DGXYHO)) '/' ('ABS' (DGXYHO)) ;
  1368.  'FINSI' ;
  1369. *****************
  1370.  'SI' (ERR1 > ERRTOL) ;
  1371.         'MESSAGE' 'Problem 58 T1';
  1372.         'ERREUR' 5 ;
  1373.  'FINSI' ;
  1374. *************************************************************
  1375.  'SI' (('ABS'(DGYYHO)) < DELTA) ;
  1376.     ERR1 = ('ABS' (DGYYHON '-' DGYYHO)) ; 
  1377.  'SINON' ;
  1378.     ERR1 = ('ABS' (DGYYHON '-' DGYYHO)) '/' ('ABS' (DGYYHO)) ;
  1379.  'FINSI' ;
  1380. *****************    
  1381.  'SI' (ERR1 > ERRTOL) ;
  1382.         'MESSAGE' 'Problem 59 T1';
  1383.         'ERREUR' 5 ;
  1384.  'FINSI' ;
  1385. **************************************************************
  1386.  'SI' (('ABS'(DGZYHO)) < DELTA) ;
  1387.     ERR1 = ('ABS' (DGZYHON '-' DGZYHO)) ; 
  1388.  'SINON' ;
  1389.     ERR1 = ('ABS' (DGZYHON '-' DGZYHO)) '/' ('ABS' (DGZYHO)) ;
  1390.  'FINSI' ;
  1391. *****************    
  1392.  'SI' (ERR1 > ERRTOL) ;
  1393.         'MESSAGE' 'Problem 60 T1';
  1394.         'ERREUR' 5 ;
  1395.  'FINSI' ; 
  1396. **************************************************************
  1397.  'SI' (('ABS'(DRETYHO)) < DELTA) ;
  1398.     ERR1 = ('ABS' (DRETYHON '-' DRETYHO)) ;
  1399.  'SINON' ;
  1400.     ERR1 = ('ABS' (DRETYHON '-' DRETYHO)) '/' ('ABS' (DRETYHO)) ;
  1401.  'FINSI' ;
  1402. *****************
  1403.  'SI' (ERR1 > ERRTOL) ;
  1404.         'MESSAGE' 'Problem 61 T1';
  1405.         'ERREUR' 5 ;
  1406.  'FINSI' ;
  1407. **************************************************************
  1408.  'SI' (('ABS'(DYH2YHO)) < DELTA) ;
  1409.     ERR1 = ('ABS' (DYH2YHON '-' DYH2YHO)) ;
  1410.  'SINON' ;  
  1411.     ERR1 = ('ABS' (DYH2YHON '-' DYH2YHO)) '/' ('ABS' (DYH2YHO)) ;
  1412.  'FINSI' ;
  1413. **************************
  1414.  'SI' (ERR1 > ERRTOL) ;
  1415.       'MESSAGE' 'Problem 62 T1';
  1416.       'ERREUR' 5 ;
  1417.  'FINSI' ;
  1418. **************************************************************
  1419.  'SI' (('ABS'(DYO2YHO)) < DELTA) ;
  1420.     ERR1 = ('ABS' (DYO2YHON '-' DYO2YHO)) ;
  1421.  'SINON' ;  
  1422.     ERR1 = ('ABS' (DYO2YHON '-' DYO2YHO)) '/' ('ABS' (DYO2YHO)) ;
  1423.  'FINSI' ;
  1424. **************************
  1425.  'SI' (ERR1 > ERRTOL) ;
  1426.       'MESSAGE' 'Problem 63 T1';
  1427.       'ERREUR' 5 ;
  1428.  'FINSI' ;
  1429. **************************************************************
  1430.  'SI' (('ABS'(DYHOYHO)) < DELTA) ;
  1431.     ERR1 = ('ABS' (DYHOYHON '-' DYHOYHO)) ;
  1432.  'SINON' ;  
  1433.     ERR1 = ('ABS' (DYHOYHON '-' DYHOYHO)) '/' ('ABS' (DYHOYHO)) ;
  1434.  'FINSI' ;
  1435. **************************
  1436.  'SI' (ERR1 > ERRTOL) ;
  1437.       'MESSAGE' 'Problem 64 T1';
  1438.       'ERREUR' 5 ;
  1439.  'FINSI' ;
  1440.  
  1441. *****************************************************
  1442. *****************************************************
  1443. ******* TEST2 ***************************************
  1444. *****************************************************
  1445. *****************************************************
  1446. *
  1447. * On compare le jacobien et la variation du residu
  1448. * en $DOM1 . 'CENTRE' par rapport à une variation
  1449. * infinitésimal en $DOM27 . 'CENTRE'
  1450. *
  1451.  
  1452.  PCEN27 = ('DOMA '$DOM27 'CENTRE') 'POIN' 1 ;
  1453.  
  1454. ****************************************************************
  1455. ****************************************************************
  1456.   DRR = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMDEN ;
  1457.   DGXR = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMMOX ;
  1458.   DGYR = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMMOY ;
  1459.   DGZR = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMMOZ ;
  1460.   DRETR = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMRET ;
  1461.   DYH2R = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMYH2 ;
  1462.   DYO2R = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMYO2 ;
  1463.   DYHOR = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMYHO ;
  1464. ***************************************
  1465.   DRGX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMDEN ;
  1466.   DGXGX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMMOX ;
  1467.   DGYGX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMMOY ;
  1468.   DGZGX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMMOZ ;
  1469.   DRETGX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMRET ;
  1470.   DYH2GX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMYH2 ;
  1471.   DYO2GX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMYO2 ;
  1472.   DYHOGX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMYHO ;
  1473. ***************************************
  1474.   DRGY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMDEN ;
  1475.   DGXGY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMMOX ;
  1476.   DGYGY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMMOY ;
  1477.   DGZGY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMMOZ ; 
  1478.   DRETGY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMRET ;
  1479.   DYH2GY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMYH2 ;
  1480.   DYO2GY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMYO2 ;
  1481.   DYHOGY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMYHO ;
  1482. ***************************************
  1483.   DRGZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMDEN ;
  1484.   DGXGZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMMOX ;
  1485.   DGYGZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMMOY ;
  1486.   DGZGZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMMOZ ; 
  1487.   DRETGZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMRET ;
  1488.   DYH2GZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMYH2 ;
  1489.   DYO2GZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMYO2 ;
  1490.   DYHOGZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMYHO ;    
  1491. ***************************************
  1492.   DRRET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMDEN ;
  1493.   DGXRET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMMOX ;
  1494.   DGYRET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMMOY ;
  1495.   DGZRET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMMOZ ;
  1496.   DRETRET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMRET ;
  1497.   DYH2RET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMYH2 ;
  1498.   DYO2RET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMYO2 ;
  1499.   DYHORET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMYHO ;
  1500. ***************************************
  1501.   DRYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMDEN ;
  1502.   DGXYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMMOX ;
  1503.   DGYYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMMOY ;
  1504.   DGZYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMMOZ ;
  1505.   DRETYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMRET ;
  1506.   DYH2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMYH2 ;
  1507.   DYO2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMYO2 ;
  1508.   DYHOYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMYHO ;
  1509. ***************************************
  1510.   DRYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMDEN ;
  1511.   DGXYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMMOX ;
  1512.   DGYYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMMOY ;
  1513.   DGZYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMMOZ ;
  1514.   DRETYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMRET ;
  1515.   DYH2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMYH2 ;
  1516.   DYO2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMYO2 ;
  1517.   DYHOYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMYHO ;
  1518. ****************************************
  1519.   DRYHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMDEN ;
  1520.   DGXYHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMMOX ;
  1521.   DGYYHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMMOY ;
  1522.   DGZYHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMMOZ ;
  1523.   DRETYHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMRET ;
  1524.   DYH2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMYH2 ;
  1525.   DYO2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMYO2 ;
  1526.   DYHOYHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMYHO ;
  1527.  
  1528.  
  1529.  
  1530. *****************************************************************
  1531. *****************************************************************
  1532. *****************************************************************  
  1533.   DRRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1534.       PCEN9 NOMDEN NOMDEN DELTA ;  
  1535.   DGXRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1536.       PCEN9 NOMDEN NOMMOX DELTA ;      
  1537.   DGYRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1538.       PCEN9 NOMDEN NOMMOY DELTA ;
  1539.   DGZRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1540.       PCEN9 NOMDEN NOMMOZ DELTA ;    
  1541.   DRETRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1542.       PCEN9 NOMDEN NOMRET DELTA ;
  1543.   DYH2RN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1544.       PCEN9 NOMDEN NOMYH2 DELTA ;
  1545.   DYO2RN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1546.       PCEN9 NOMDEN NOMYO2 DELTA ;
  1547.   DYHORN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1548.       PCEN9 NOMDEN NOMYHO DELTA ;
  1549. **************************************
  1550.   DRGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1551.       PCEN9 NOMMOX NOMDEN DELTA ;  
  1552.   DGXGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1553.       PCEN9 NOMMOX NOMMOX DELTA ;      
  1554.   DGYGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1555.       PCEN9 NOMMOX NOMMOY DELTA ;
  1556.   DGZGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1557.       PCEN9 NOMMOX NOMMOZ DELTA ;    
  1558.   DRETGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1559.       PCEN9 NOMMOX NOMRET DELTA ;
  1560.   DYH2GXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1561.       PCEN9 NOMMOX NOMYH2 DELTA ;
  1562.   DYO2GXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1563.       PCEN9 NOMMOX NOMYO2 DELTA ;
  1564.   DYHOGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1565.       PCEN9 NOMMOX NOMYHO DELTA ;
  1566. **************************************
  1567.   DRGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1568.       PCEN9 NOMMOY NOMDEN DELTA ;  
  1569.   DGXGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1570.       PCEN9 NOMMOY NOMMOX DELTA ;      
  1571.   DGYGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1572.       PCEN9 NOMMOY NOMMOY DELTA ;
  1573.   DGZGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1574.       PCEN9 NOMMOY NOMMOZ DELTA ;    
  1575.   DRETGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1576.       PCEN9 NOMMOY NOMRET DELTA ;
  1577.   DYH2GYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1578.       PCEN9 NOMMOY NOMYH2 DELTA ;
  1579.   DYO2GYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1580.       PCEN9 NOMMOY NOMYO2 DELTA ;
  1581.   DYHOGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1582.       PCEN9 NOMMOY NOMYHO DELTA ;
  1583. **************************************
  1584.   DRGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1585.       PCEN9 NOMMOZ NOMDEN DELTA ;  
  1586.   DGXGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1587.       PCEN9 NOMMOZ NOMMOX DELTA ;      
  1588.   DGYGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1589.       PCEN9 NOMMOZ NOMMOY DELTA ;
  1590.   DGZGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1591.       PCEN9 NOMMOZ NOMMOZ DELTA ;    
  1592.   DRETGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1593.       PCEN9 NOMMOZ NOMRET DELTA ;
  1594.   DYH2GZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1595.       PCEN9 NOMMOZ NOMYH2 DELTA ;
  1596.   DYO2GZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1597.       PCEN9 NOMMOZ NOMYO2 DELTA ;
  1598.   DYHOGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1599.       PCEN9 NOMMOZ NOMYHO DELTA ;      
  1600. **************************************
  1601.   DRRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1602.       PCEN9 NOMRET NOMDEN DELTA ;  
  1603.   DGXRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1604.       PCEN9 NOMRET NOMMOX DELTA ;      
  1605.   DGYRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1606.       PCEN9 NOMRET NOMMOY DELTA ;
  1607.   DGZRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1608.       PCEN9 NOMRET NOMMOZ DELTA ;    
  1609.   DRETRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1610.       PCEN9 NOMRET NOMRET DELTA ;
  1611.   DYH2RETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1612.       PCEN9 NOMRET NOMYH2 DELTA ;
  1613.   DYO2RETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1614.       PCEN9 NOMRET NOMYO2 DELTA ;
  1615.   DYHORETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1616.       PCEN9 NOMRET NOMYHO DELTA ;
  1617. **************************************
  1618.   DRYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1619.       PCEN9 NOMYH2 NOMDEN DELTA ;  
  1620.   DGXYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1621.       PCEN9 NOMYH2 NOMMOX DELTA ;      
  1622.   DGYYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1623.       PCEN9 NOMYH2 NOMMOY DELTA ;
  1624.   DGZYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1625.       PCEN9 NOMYH2 NOMMOZ DELTA ;    
  1626.   DRETYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1627.       PCEN9 NOMYH2 NOMRET DELTA ;
  1628.   DYH2YH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1629.       PCEN9 NOMYH2 NOMYH2 DELTA ;
  1630.   DYO2YH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1631.       PCEN9 NOMYH2 NOMYO2 DELTA ;
  1632.   DYHOYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1633.       PCEN9 NOMYH2 NOMYHO DELTA ;
  1634. ***************************************
  1635.   DRYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1636.       PCEN9 NOMYO2 NOMDEN DELTA ;  
  1637.   DGXYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1638.       PCEN9 NOMYO2 NOMMOX DELTA ;      
  1639.   DGYYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1640.       PCEN9 NOMYO2 NOMMOY DELTA ;
  1641.   DGZYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1642.       PCEN9 NOMYO2 NOMMOZ DELTA ;    
  1643.   DRETYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1644.       PCEN9 NOMYO2 NOMRET DELTA ;
  1645.   DYH2YO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1646.       PCEN9 NOMYO2 NOMYH2 DELTA ;
  1647.   DYO2YO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1648.       PCEN9 NOMYO2 NOMYO2 DELTA ;
  1649.   DYHOYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1650.       PCEN9 NOMYO2 NOMYHO DELTA ;
  1651. ***************************************
  1652.   DRYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1653.       PCEN9 NOMYHO NOMDEN DELTA ;  
  1654.   DGXYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1655.       PCEN9 NOMYHO NOMMOX DELTA ;      
  1656.   DGYYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1657.       PCEN9 NOMYHO NOMMOY DELTA ;
  1658.   DGZYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1659.       PCEN9 NOMYHO NOMMOZ DELTA ;    
  1660.   DRETYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1661.       PCEN9 NOMYHO NOMRET DELTA ;
  1662.   DYH2YHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1663.       PCEN9 NOMYHO NOMYH2 DELTA ;
  1664.   DYO2YHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1665.       PCEN9 NOMYHO NOMYO2 DELTA ;
  1666.   DYHOYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  1667.       PCEN9 NOMYHO NOMYHO DELTA ;   
  1668.  
  1669. * Test des comparaisons jacobien exact-jacobien numerique
  1670. * Test des comparaisons jacobien exact-jacobien numerique
  1671. **************************************************************
  1672.  'SI' (('ABS'(DRR)) < DELTA) ;
  1673.     ERR1 = ('ABS' (DRRN '-' DRR)) ;
  1674.  'SINON' ;  
  1675.     ERR1 = ('ABS' (DRRN '-' DRR)) '/' ('ABS' (DRR)) ;
  1676.  'FINSI' ; 
  1677. ***********************
  1678.  'SI' (ERR1 > ERRTOL) ;
  1679.       'MESSAGE' 'Problem 1 T2';
  1680.       'ERREUR' 5 ;
  1681.  'FINSI' ;
  1682. **************************************************************
  1683.  'SI' (('ABS'(DGXR)) < DELTA) ;
  1684.     ERR1 = ('ABS' (DGXRN '-' DGXR)) ;
  1685.  'SINON' ;  
  1686.     ERR1 = ('ABS' (DGXRN '-' DGXR)) '/' ('ABS' (DGXR)) ;
  1687.  'FINSI' ; 
  1688. *********************** 
  1689.  'SI' (ERR1 > ERRTOL) ;
  1690.       'MESSAGE' 'Problem 2 T2';
  1691.       'ERREUR' 5 ;
  1692.  'FINSI' ;
  1693. **************************************************************
  1694.  'SI' (('ABS'(DGYR)) < DELTA) ;
  1695.     ERR1 = ('ABS' (DGYRN '-' DGYR)) ;
  1696.  'SINON' ;  
  1697.     ERR1 = ('ABS' (DGYRN '-' DGYR)) '/' ('ABS' (DGYR)) ;
  1698.  'FINSI' ; 
  1699. ************************
  1700.  'SI' (ERR1 > ERRTOL) ;
  1701.       'MESSAGE' 'Problem 3 T2';
  1702.        'ERREUR' 5 ;
  1703.  'FINSI' ;
  1704. **************************************************************
  1705.  'SI' (('ABS'(DGZR)) < DELTA) ;
  1706.     ERR1 = ('ABS' (DGZRN '-' DGZR)) ;
  1707.  'SINON' ;  
  1708.     ERR1 = ('ABS' (DGZRN '-' DGZR)) '/' ('ABS' (DGZR)) ;
  1709.  'FINSI' ;
  1710. ************************
  1711.  'SI' (ERR1 > ERRTOL) ;
  1712.       'MESSAGE' 'Problem 4 T2';
  1713.        'ERREUR' 5 ;
  1714.  'FINSI' ;
  1715. ***************************************************************
  1716.  'SI' (('ABS'(DRETR)) < DELTA) ;
  1717.     ERR1 = ('ABS' (DRETRN '-' DRETR)) ;
  1718.  'SINON' ;  
  1719.     ERR1 = ('ABS' (DRETRN '-' DRETR)) '/' ('ABS' (DRETR)) ;
  1720.  'FINSI' ;
  1721. **************************
  1722.  'SI' (ERR1 > ERRTOL) ;
  1723.       'MESSAGE' 'Problem 5 T2';
  1724.       'ERREUR' 5 ;
  1725.  'FINSI' ;
  1726. **************************************************************
  1727. **************************************************************
  1728.  'SI' (('ABS'(DYH2R)) < DELTA) ;
  1729.     ERR1 = ('ABS' (DYH2RN '-' DYH2R)) ;
  1730.  'SINON' ;  
  1731.     ERR1 = ('ABS' (DYH2RN '-' DYH2R)) '/' ('ABS' (DYH2R)) ;
  1732.  'FINSI' ;
  1733. **************************
  1734.  'SI' (ERR1 > ERRTOL) ;
  1735.       'MESSAGE' 'Problem 6 T2';
  1736.       'ERREUR' 5 ;
  1737.  'FINSI' ;
  1738. **************************************************************
  1739.    'SI' (('ABS'(DYO2R)) < DELTA) ;
  1740.     ERR1 = ('ABS' (DYO2RN '-' DYO2R)) ;
  1741.  'SINON' ;  
  1742.     ERR1 = ('ABS' (DYO2RN '-' DYO2R)) '/' ('ABS' (DYO2R)) ;
  1743.  'FINSI' ;
  1744. **************************
  1745.  'SI' (ERR1 > ERRTOL) ;
  1746.       'MESSAGE' 'Problem 7 T2';
  1747.       'ERREUR' 5 ;
  1748.  'FINSI' ;
  1749. **************************************************************
  1750.   'SI' (('ABS'(DYHOR)) < DELTA) ;
  1751.     ERR1 = ('ABS' (DYHORN '-' DYHOR)) ;
  1752.  'SINON' ;  
  1753.     ERR1 = ('ABS' (DYHORN '-' DYHOR)) '/' ('ABS' (DYHOR)) ;
  1754.  'FINSI' ;
  1755. **************************
  1756.  'SI' (ERR1 > ERRTOL) ;
  1757.       'MESSAGE' 'Problem 8 T2';
  1758.       'ERREUR' 5 ;
  1759.  'FINSI' ; 
  1760. **************************************************************
  1761. **************************************************************
  1762. **************************************************************
  1763. *****    Second line !!!!!!!!!!!!!!!!!!!!!!!    **************
  1764. **************************************************************
  1765. ************************************************************** 
  1766.  'SI' (('ABS'(DRGX)) < DELTA) ;
  1767.     ERR1 = ('ABS' (DRGXN '-' DRGX)) ;
  1768.  'SINON' ;   
  1769.     ERR1 = ('ABS' (DRGXN '-' DRGX)) '/' ('ABS' (DRGX))  ;
  1770.  'FINSI' ;
  1771. *************************** 
  1772.  'SI' (ERR1 > ERRTOL) ;
  1773.       'MESSAGE' 'Problem 9 T2';
  1774.       'ERREUR' 5 ;
  1775.  'FINSI' ;
  1776. ***************************************************************
  1777.  'SI' (('ABS'(DGXGX)) < DELTA) ;
  1778.    ERR1 = ('ABS' (DGXGXN '-' DGXGX)) ;
  1779.  'SINON' ;
  1780.    ERR1 = ('ABS' (DGXGXN '-' DGXGX)) '/' ('ABS' (DGXGX)) ;
  1781.  'FINSI' ;  
  1782. *****************************   
  1783.  'SI' (ERR1 > ERRTOL) ;
  1784.         'MESSAGE' 'Problem 10 T2';
  1785.         'ERREUR' 5 ;
  1786.  'FINSI' ;
  1787. **************************************************************
  1788.  'SI' (('ABS'(DGYGX)) < DELTA) ;
  1789.     ERR1 = ('ABS' (DGYGXN '-' DGYGX)) ;
  1790.  'SINON' ;
  1791.     ERR1 = ('ABS' (DGYGXN '-' DGYGX)) '/' ('ABS' (DGYGX)) ;
  1792.  'FINSI' ;
  1793. ******************************     
  1794.  'SI' (ERR1 > ERRTOL) ;
  1795.         'MESSAGE' 'Problem 11 T2';
  1796.         'ERREUR' 5 ;
  1797.  'FINSI' ;
  1798. **************************************************************
  1799.  'SI' (('ABS'(DGZGX)) < DELTA) ;
  1800.     ERR1 = ('ABS' (DGZGXN '-' DGZGX)) ;
  1801.  'SINON' ;
  1802.     ERR1 = ('ABS' (DGZGXN '-' DGZGX)) '/' ('ABS' (DGZGX)) ;
  1803.  'FINSI' ;
  1804. ******************************    
  1805.  'SI' (ERR1 > ERRTOL) ;
  1806.         'MESSAGE' 'Problem 12 T2';
  1807.         'ERREUR' 5 ;
  1808.  'FINSI' ;
  1809. **************************************************************
  1810.  'SI' (('ABS'(DRETGX)) < DELTA) ;
  1811.     ERR1 = ('ABS' (DRETGXN '-' DRETGX)) ;
  1812.  'SINON' ;
  1813.     ERR1 = ('ABS' (DRETGXN '-' DRETGX)) '/' ('ABS' (DRETGX)) ;
  1814.  'FINSI' ;
  1815. *******************************    
  1816.  'SI' (ERR1 > ERRTOL) ;
  1817.         'MESSAGE' 'Problem 13 T2';
  1818.         'ERREUR' 5 ;
  1819.  'FINSI' ;
  1820. **************************************************************
  1821.  'SI' (('ABS'(DYH2GX)) < DELTA) ;
  1822.     ERR1 = ('ABS' (DYH2GXN '-' DYH2GX)) ;
  1823.  'SINON' ;  
  1824.     ERR1 = ('ABS' (DYH2GXN '-' DYH2GX)) '/' ('ABS' (DYH2GX)) ;
  1825.  'FINSI' ;
  1826. **************************
  1827.  'SI' (ERR1 > ERRTOL) ;
  1828.       'MESSAGE' 'Problem 14 T2';
  1829.       'ERREUR' 5 ;
  1830.  'FINSI' ;
  1831. **************************************************************
  1832.  'SI' (('ABS'(DYO2GX)) < DELTA) ;
  1833.     ERR1 = ('ABS' (DYO2GXN '-' DYO2GX)) ;
  1834.  'SINON' ;  
  1835.     ERR1 = ('ABS' (DYO2GXN '-' DYO2GX)) '/' ('ABS' (DYO2GX)) ;
  1836.  'FINSI' ;
  1837. **************************
  1838.  'SI' (ERR1 > ERRTOL) ;
  1839.       'MESSAGE' 'Problem 15 T2';
  1840.       'ERREUR' 5 ;
  1841.  'FINSI' ;
  1842. **************************************************************
  1843.  'SI' (('ABS'(DYHOGX)) < DELTA) ;
  1844.     ERR1 = ('ABS' (DYHOGXN '-' DYHOGX)) ;
  1845.  'SINON' ;  
  1846.     ERR1 = ('ABS' (DYHOGXN '-' DYHOGX)) '/' ('ABS' (DYHOGX)) ;
  1847.  'FINSI' ;
  1848. **************************
  1849.  'SI' (ERR1 > ERRTOL) ;
  1850.       'MESSAGE' 'Problem 16 T2';
  1851.       'ERREUR' 5 ;
  1852.  'FINSI' ;
  1853. **************************************************************
  1854. *******        Third line!!!!!!!!!!!!!!!!      ***************
  1855. **************************************************************
  1856.  'SI' (('ABS'(DRGY)) < DELTA) ;
  1857.     ERR1 = ('ABS' (DRGYN '-' DRGY)) ;
  1858.  'SINON' ;
  1859.     ERR1 = ('ABS' (DRGYN '-' DRGY)) '/' ('ABS' (DRGY)) ;
  1860.  'FINSI' ;
  1861. *******************************       
  1862.  'SI' (ERR1 > ERRTOL) ;
  1863.         'MESSAGE' 'Problem 17 T2';
  1864.         'ERREUR' 5 ;
  1865.  'FINSI' ;
  1866. **************************************************************
  1867.  'SI' (('ABS'(DGXGY)) < DELTA) ;
  1868.     ERR1 = ('ABS' (DGXGYN '-' DGXGY)) ;
  1869.  'SINON' ;
  1870.     ERR1 = ('ABS' (DGXGYN '-' DGXGY))  '/' ('ABS' (DGXGY)) ;
  1871.  'FINSI' ;
  1872. ********************************      
  1873.  'SI' (ERR1 > ERRTOL) ;
  1874.         'MESSAGE' 'Problem 18 T2';
  1875.         'ERREUR' 5 ;
  1876.  'FINSI' ;
  1877. **************************************************************
  1878.  'SI' (('ABS'(DGYGY)) < DELTA) ;
  1879.     ERR1 = ('ABS' (DGYGYN '-' DGYGY)) ;
  1880.  'SINON' ;
  1881.     ERR1 = ('ABS' (DGYGYN '-' DGYGY)) '/' ('ABS' (DGYGY)) ;
  1882.  'FINSI' ;
  1883. ***************
  1884.  'SI' (ERR1 > ERRTOL) ;
  1885.         'MESSAGE' 'Problem 19 T2';
  1886.         'ERREUR' 5 ;
  1887.  'FINSI' ;
  1888. **************************************************************
  1889.  'SI' (('ABS'(DGZGY)) < DELTA) ;
  1890.     ERR1 = ('ABS' (DGZGYN '-' DGZGY)) ;
  1891.  'SINON' ;
  1892.     ERR1 = ('ABS' (DGZGYN '-' DGZGY)) '/' ('ABS' (DGZGY)) ;
  1893.  'FINSI' ;
  1894. ****************        
  1895.  'SI' (ERR1 > ERRTOL) ;
  1896.         'MESSAGE' 'Problem 20 T2';
  1897.         'ERREUR' 5 ;
  1898.  'FINSI' ;
  1899. **************************************************************
  1900.  'SI' (('ABS'(DRETGY)) < DELTA) ;
  1901.     ERR1 = ('ABS' (DRETGYN '-' DRETGY)) ;
  1902.  'SINON' ;
  1903.     ERR1 = ('ABS' (DRETGYN '-' DRETGY)) '/' ('ABS' (DRETGY)) ;
  1904.  'FINSI' ;
  1905. ****************     
  1906.  'SI' (ERR1 > ERRTOL) ;
  1907.         'MESSAGE' 'Problem 21 T2';
  1908.         'ERREUR' 5 ;
  1909.  'FINSI' ;
  1910. **************************************************************
  1911.  'SI' (('ABS'(DYH2GY)) < DELTA) ;
  1912.     ERR1 = ('ABS' (DYH2GYN '-' DYH2GY)) ;
  1913.  'SINON' ;  
  1914.     ERR1 = ('ABS' (DYH2GYN '-' DYH2GY)) '/' ('ABS' (DYH2GY)) ;
  1915.  'FINSI' ;
  1916. **************************
  1917.  'SI' (ERR1 > ERRTOL) ;
  1918.       'MESSAGE' 'Problem 22 T2';
  1919.       'ERREUR' 5 ;
  1920.  'FINSI' ;
  1921. **************************************************************
  1922.  'SI' (('ABS'(DYO2GY)) < DELTA) ;
  1923.     ERR1 = ('ABS' (DYO2GYN '-' DYO2GY)) ;
  1924.  'SINON' ;  
  1925.     ERR1 = ('ABS' (DYO2GYN '-' DYO2GY)) '/' ('ABS' (DYO2GY)) ;
  1926.  'FINSI' ;
  1927. **************************
  1928.  'SI' (ERR1 > ERRTOL) ;
  1929.       'MESSAGE' 'Problem 23 T2';
  1930.       'ERREUR' 5 ;
  1931.  'FINSI' ;
  1932. **************************************************************
  1933.  'SI' (('ABS'(DYHOGY)) < DELTA) ;
  1934.     ERR1 = ('ABS' (DYHOGYN '-' DYHOGY)) ;
  1935.  'SINON' ;  
  1936.     ERR1 = ('ABS' (DYHOGYN '-' DYHOGY)) '/' ('ABS' (DYHOGY)) ;
  1937.  'FINSI' ;
  1938. **************************
  1939.  'SI' (ERR1 > ERRTOL) ;
  1940.       'MESSAGE' 'Problem 24 T2';
  1941.       'ERREUR' 5 ;
  1942.  'FINSI' ;
  1943. **************************************************************
  1944. *******     Fouth line!!!!!!!!!!!!!!!!!!!   ******************
  1945. **************************************************************
  1946.  'SI' (('ABS'(DRGZ)) < DELTA) ;
  1947.     ERR1 = ('ABS' (DRGZN '-' DRGZ))  ;
  1948.  'SINON' ;
  1949.     ERR1 = ('ABS' (DRGZN '-' DRGZ)) '/' ('ABS' (DRGZ))  ;
  1950.  'FINSI' ;
  1951. ***************      
  1952.  'SI' (ERR1 > ERRTOL) ;
  1953.         'MESSAGE' 'Problem 25 T2';
  1954.         'ERREUR' 5 ;
  1955.  'FINSI' ;
  1956. **************************************************************
  1957.  'SI' (('ABS'(DGXGZ)) < DELTA) ;
  1958.     ERR1 = ('ABS' (DGXGZN '-' DGXGZ)) ;
  1959.  'SINON' ;
  1960.     ERR1 = ('ABS' (DGXGZN '-' DGXGZ))  '/' ('ABS' (DGXGZ)) ;
  1961.  'FINSI' ;
  1962. *****************       
  1963.  'SI' (ERR1 > ERRTOL) ;
  1964.         'MESSAGE' 'Problem 26 T2';
  1965.         'ERREUR' 5 ;
  1966.  'FINSI' ;
  1967. **************************************************************
  1968.  'SI' (('ABS'(DGYGZ)) < DELTA) ;
  1969.     ERR1 = ('ABS' (DGYGZN '-' DGYGZ)) ;
  1970.  'SINON' ;
  1971.     ERR1 = ('ABS' (DGYGZN '-' DGYGZ)) '/' ('ABS' (DGYGZ)) ;
  1972.  'FINSI' ;
  1973. *****************       
  1974.  'SI' (ERR1 > ERRTOL) ;
  1975.         'MESSAGE' 'Problem 27 T2';
  1976.         'ERREUR' 5 ;
  1977.  'FINSI' ;
  1978. **************************************************************
  1979.  'SI' (('ABS'(DGZGZ)) < DELTA) ;
  1980.     ERR1 = ('ABS' (DGZGZN '-' DGZGZ)) ;
  1981.  'SINON' ;
  1982.     ERR1 = ('ABS' (DGZGZN '-' DGZGZ)) '/' ('ABS' (DGZGZ)) ;
  1983.  'FINSI' ;
  1984. *****************     
  1985.  'SI' (ERR1 > ERRTOL) ;
  1986.         'MESSAGE' 'Problem 28 T2';
  1987.         'ERREUR' 5 ;
  1988.  'FINSI' ;
  1989. **************************************************************
  1990.  'SI' (('ABS'(DRETGZ)) < DELTA) ;
  1991.     ERR1 = ('ABS' (DRETGZN '-' DRETGZ)) ;
  1992.  'SINON' ;
  1993.     ERR1 = ('ABS' (DRETGZN '-' DRETGZ)) '/' ('ABS' (DRETGZ)) ;
  1994.  'FINSI' ;     
  1995. *****************
  1996.  'SI' (ERR1 > ERRTOL) ;
  1997.         'MESSAGE' 'Problem 29 T2';
  1998.         'ERREUR' 5 ;
  1999.  'FINSI' ;
  2000. **************************************************************
  2001.  'SI' (('ABS'(DYH2GZ)) < DELTA) ;
  2002.     ERR1 = ('ABS' (DYH2GZN '-' DYH2GZ)) ;
  2003.  'SINON' ;  
  2004.     ERR1 = ('ABS' (DYH2GZN '-' DYH2GZ)) '/' ('ABS' (DYH2GZ)) ;
  2005.  'FINSI' ;
  2006. **************************
  2007.  'SI' (ERR1 > ERRTOL) ;
  2008.       'MESSAGE' 'Problem 30 T2';
  2009.       'ERREUR' 5 ;
  2010.  'FINSI' ;
  2011. **************************************************************
  2012.  'SI' (('ABS'(DYO2GZ)) < DELTA) ;
  2013.     ERR1 = ('ABS' (DYO2GZN '-' DYO2GZ)) ;
  2014.  'SINON' ;  
  2015.     ERR1 = ('ABS' (DYO2GZN '-' DYO2GZ)) '/' ('ABS' (DYO2GZ)) ;
  2016.  'FINSI' ;
  2017. **************************
  2018.  'SI' (ERR1 > ERRTOL) ;
  2019.       'MESSAGE' 'Problem 31 T2';
  2020.       'ERREUR' 5 ;
  2021.  'FINSI' ;
  2022. **************************************************************
  2023.  'SI' (('ABS'(DYHOGZ)) < DELTA) ;
  2024.     ERR1 = ('ABS' (DYHOGZN '-' DYHOGZ)) ;
  2025.  'SINON' ;  
  2026.     ERR1 = ('ABS' (DYHOGZN '-' DYHOGZ)) '/' ('ABS' (DYHOGZ)) ;
  2027.  'FINSI' ;
  2028. **************************
  2029.  'SI' (ERR1 > ERRTOL) ;
  2030.       'MESSAGE' 'Problem 32 T2';
  2031.       'ERREUR' 5 ;
  2032.  'FINSI' ;
  2033. **************************************************************
  2034. *******     Fith line!!!!!!!!!!!!!!!!!!!   ******************
  2035. **************************************************************
  2036.  'SI' (('ABS'(DRRET)) < DELTA) ;
  2037.     ERR1 = ('ABS' (DRRETN '-' DRRET)) ;
  2038.  'SINON' ;
  2039.     ERR1 = ('ABS' (DRRETN '-' DRRET)) '/' ('ABS' (DRRET)) ;
  2040.  'FINSI' ;
  2041. *****************       
  2042.  'SI' (ERR1 > ERRTOL) ;
  2043.         'MESSAGE' 'Problem 33 T2';
  2044.         'ERREUR' 5 ;
  2045.  'FINSI' ;
  2046. ***************************************************************
  2047.  'SI' (('ABS'(DGXRET)) < DELTA) ;
  2048.     ERR1 = ('ABS' (DGXRETN '-' DGXRET)) ;
  2049.  'SINON' ;
  2050.     ERR1 = ('ABS' (DGXRETN '-' DGXRET)) '/' ('ABS' (DGXRET)) ;
  2051.  'FINSI' ;
  2052. *****************    
  2053.  'SI' (ERR1 > ERRTOL) ;
  2054.         'MESSAGE' 'Problem 34 T2';
  2055.         'ERREUR' 5 ;
  2056.  'FINSI' ;
  2057. ****************************************************************
  2058.  'SI' (('ABS'(DGYRET)) < DELTA) ;
  2059.     ERR1 = ('ABS' (DGYRETN '-' DGYRET)) ; 
  2060.  'SINON' ;
  2061.     ERR1 = ('ABS' (DGYRETN '-' DGYRET)) '/' ('ABS' (DGYRET)) ;
  2062.  'FINSI' ;
  2063. *****************    
  2064.  'SI' (ERR1 > ERRTOL) ;
  2065.         'MESSAGE' 'Problem 35 T2';
  2066.         'ERREUR' 5 ;
  2067.  'FINSI' ;
  2068. ****************************************************************
  2069.  'SI' (('ABS'(DGZRET)) < DELTA) ;
  2070.     ERR1 = ('ABS' (DGZRETN '-' DGZRET)) ;
  2071.  'SINON' ;    
  2072.     ERR1 = ('ABS' (DGZRETN '-' DGZRET)) '/' ('ABS' (DGZRET)) ;
  2073.  'FINSI' ;
  2074. *****************    
  2075.  'SI' (ERR1 > ERRTOL) ;
  2076.         'MESSAGE' 'Problem 36 T2';
  2077.         'ERREUR' 5 ;
  2078.  'FINSI' ;
  2079. ****************************************************************
  2080.  'SI' (('ABS'(DRETRET)) < DELTA) ;
  2081.     ERR1 = ('ABS' (DRETRETN '-' DRETRET)) ;
  2082.  'SINON' ;
  2083.     ERR1 = ('ABS' (DRETRETN '-' DRETRET)) '/' ('ABS' (DRETRET)) ;
  2084.  'FINSI' ;
  2085. *****************       
  2086.  'SI' (ERR1 > ERRTOL) ;
  2087.         'MESSAGE' 'Problem 37 T2';
  2088.         'ERREUR' 5 ;
  2089.  'FINSI' ;
  2090.  'SI' (('ABS'(DYH2RET)) < DELTA) ;
  2091.     ERR1 = ('ABS' (DYH2RETN '-' DYH2RET)) ;
  2092.  'SINON' ;  
  2093.     ERR1 = ('ABS' (DYH2RETN '-' DYH2RET)) '/' ('ABS' (DYH2RET)) ;
  2094.  'FINSI' ;
  2095. **************************
  2096.  'SI' (ERR1 > ERRTOL) ;
  2097.       'MESSAGE' 'Problem 38 T2';
  2098.       'ERREUR' 5 ;
  2099.  'FINSI' ;
  2100. **************************************************************
  2101.  'SI' (('ABS'(DYO2RET)) < DELTA) ;
  2102.     ERR1 = ('ABS' (DYO2RETN '-' DYO2RET)) ;
  2103.  'SINON' ;  
  2104.     ERR1 = ('ABS' (DYO2RETN '-' DYO2RET)) '/' ('ABS' (DYO2RET)) ;
  2105.  'FINSI' ;
  2106. **************************
  2107.  'SI' (ERR1 > ERRTOL) ;
  2108.       'MESSAGE' 'Problem 39 T2';
  2109.       'ERREUR' 5 ;
  2110.  'FINSI' ;
  2111. **************************************************************
  2112.  'SI' (('ABS'(DYHORET)) < DELTA) ;
  2113.     ERR1 = ('ABS' (DYHORETN '-' DYHORET)) ;
  2114.  'SINON' ;  
  2115.     ERR1 = ('ABS' (DYHORETN '-' DYHORET)) '/' ('ABS' (DYHORET)) ;
  2116.  'FINSI' ;
  2117. **************************
  2118.  'SI' (ERR1 > ERRTOL) ;
  2119.       'MESSAGE' 'Problem 40 T2';
  2120.       'ERREUR' 5 ;
  2121.  'FINSI' ;
  2122. **************************************************************
  2123. *******     Sixth line!!!!!!!!!!!!!!!!!!!   ******************
  2124. **************************************************************
  2125.  'SI' (('ABS'(DRYH2)) < DELTA) ;
  2126.     ERR1 = ('ABS' (DRYH2N '-' DRYH2)) ;
  2127.  'SINON' ;
  2128.     ERR1 = ('ABS' (DRYH2N '-' DRYH2)) '/' ('ABS' (DRYH2)) ;
  2129.  'FINSI' ;
  2130. *****************
  2131.  'SI' (ERR1 > ERRTOL) ;
  2132.         'MESSAGE' 'Problem 41 T2';
  2133.         'ERREUR' 5 ;
  2134.  'FINSI' ;
  2135. **************************************************************
  2136.  'SI' (('ABS'(DGXYH2)) < DELTA) ;
  2137.     ERR1 = ('ABS' (DGXYH2N '-' DGXYH2)) ;
  2138.  'SINON' ;
  2139.     ERR1 = ('ABS' (DGXYH2N '-' DGXYH2)) '/' ('ABS' (DGXYH2)) ;
  2140.  'FINSI' ;
  2141. *****************
  2142.  'SI' (ERR1 > ERRTOL) ;
  2143.         'MESSAGE' 'Problem 42 T2';
  2144.         'ERREUR' 5 ;
  2145.  'FINSI' ;
  2146. *************************************************************
  2147.  'SI' (('ABS'(DGYYH2)) < DELTA) ;
  2148.     ERR1 = ('ABS' (DGYYH2N '-' DGYYH2)) ; 
  2149.  'SINON' ;
  2150.     ERR1 = ('ABS' (DGYYH2N '-' DGYYH2)) '/' ('ABS' (DGYYH2)) ;
  2151.  'FINSI' ;
  2152. *****************    
  2153.  'SI' (ERR1 > ERRTOL) ;
  2154.         'MESSAGE' 'Problem 43 T2';
  2155.         'ERREUR' 5 ;
  2156.  'FINSI' ;
  2157. **************************************************************
  2158.  'SI' (('ABS'(DGZYH2)) < DELTA) ;
  2159.     ERR1 = ('ABS' (DGZYH2N '-' DGZYH2)) ; 
  2160.  'SINON' ;
  2161.     ERR1 = ('ABS' (DGZYH2N '-' DGZYH2)) '/' ('ABS' (DGZYH2)) ;
  2162.  'FINSI' ;
  2163. *****************    
  2164.  'SI' (ERR1 > ERRTOL) ;
  2165.         'MESSAGE' 'Problem 44 T2';
  2166.         'ERREUR' 5 ;
  2167.  'FINSI' ; 
  2168. **************************************************************
  2169.  'SI' (('ABS'(DRETYH2)) < DELTA) ;
  2170.     ERR1 = ('ABS' (DRETYH2N '-' DRETYH2)) ;
  2171.  'SINON' ;
  2172.     ERR1 = ('ABS' (DRETYH2N '-' DRETYH2)) '/' ('ABS' (DRETYH2)) ;
  2173.  'FINSI' ;
  2174. *****************
  2175.  'SI' (ERR1 > ERRTOL) ;
  2176.         'MESSAGE' 'Problem 45 T2';
  2177.         'ERREUR' 5 ;
  2178.  'FINSI' ;
  2179. **************************************************************
  2180.  'SI' (('ABS'(DYH2YH2)) < DELTA) ;
  2181.     ERR1 = ('ABS' (DYH2YH2N '-' DYH2YH2)) ;
  2182.  'SINON' ;  
  2183.     ERR1 = ('ABS' (DYH2YH2N '-' DYH2YH2)) '/' ('ABS' (DYH2YH2)) ;
  2184.  'FINSI' ;
  2185. **************************
  2186.  'SI' (ERR1 > ERRTOL) ;
  2187.       'MESSAGE' 'Problem 46 T2';
  2188.       'ERREUR' 5 ;
  2189.  'FINSI' ;
  2190. **************************************************************
  2191.  'SI' (('ABS'(DYO2YH2)) < DELTA) ;
  2192.     ERR1 = ('ABS' (DYO2YH2N '-' DYO2YH2)) ;
  2193.  'SINON' ;  
  2194.     ERR1 = ('ABS' (DYO2YH2N '-' DYO2YH2)) '/' ('ABS' (DYO2YH2)) ;
  2195.  'FINSI' ;
  2196. **************************
  2197.  'SI' (ERR1 > ERRTOL) ;
  2198.       'MESSAGE' 'Problem 47 T2';
  2199.       'ERREUR' 5 ;
  2200.  'FINSI' ;
  2201. **************************************************************
  2202.  'SI' (('ABS'(DYHOYH2)) < DELTA) ;
  2203.     ERR1 = ('ABS' (DYHOYH2N '-' DYHOYH2)) ;
  2204.  'SINON' ;  
  2205.     ERR1 = ('ABS' (DYHOYH2N '-' DYHOYH2)) '/' ('ABS' (DYHOYH2)) ;
  2206.  'FINSI' ;
  2207. **************************
  2208.  'SI' (ERR1 > ERRTOL) ;
  2209.       'MESSAGE' 'Problem 48 T2';
  2210.       'ERREUR' 5 ;
  2211.  'FINSI' ;
  2212. **************************************************************
  2213. *********  Seventh line !!!!!!!!!!!!!!!!!!!!!!   *************
  2214. **************************************************************
  2215.  'SI' (('ABS'(DRYO2)) < DELTA) ;
  2216.     ERR1 = ('ABS' (DRYO2N '-' DRYO2)) ;
  2217.  'SINON' ;
  2218.     ERR1 = ('ABS' (DRYO2N '-' DRYO2)) '/' ('ABS' (DRYO2)) ;
  2219.  'FINSI' ;
  2220. *****************
  2221.  'SI' (ERR1 > ERRTOL) ;
  2222.         'MESSAGE' 'Problem 49 T2';
  2223.         'ERREUR' 5 ;
  2224.  'FINSI' ;
  2225. **************************************************************
  2226.  'SI' (('ABS'(DGXYO2)) < DELTA) ;
  2227.     ERR1 = ('ABS' (DGXYO2N '-' DGXYO2)) ;
  2228.  'SINON' ;
  2229.     ERR1 = ('ABS' (DGXYO2N '-' DGXYO2)) '/' ('ABS' (DGXYO2)) ;
  2230.  'FINSI' ;
  2231. *****************
  2232.  'SI' (ERR1 > ERRTOL) ;
  2233.         'MESSAGE' 'Problem 50 T2';
  2234.         'ERREUR' 5 ;
  2235.  'FINSI' ;
  2236. *************************************************************
  2237.  'SI' (('ABS'(DGYYO2)) < DELTA) ;
  2238.     ERR1 = ('ABS' (DGYYO2N '-' DGYYO2)) ; 
  2239.  'SINON' ;
  2240.     ERR1 = ('ABS' (DGYYO2N '-' DGYYO2)) '/' ('ABS' (DGYYO2)) ;
  2241.  'FINSI' ;
  2242. *****************    
  2243.  'SI' (ERR1 > ERRTOL) ;
  2244.         'MESSAGE' 'Problem 51 T2';
  2245.         'ERREUR' 5 ;
  2246.  'FINSI' ;
  2247. **************************************************************
  2248.  'SI' (('ABS'(DGZYO2)) < DELTA) ;
  2249.     ERR1 = ('ABS' (DGZYO2N '-' DGZYO2)) ; 
  2250.  'SINON' ;
  2251.     ERR1 = ('ABS' (DGZYO2N '-' DGZYO2)) '/' ('ABS' (DGZYO2)) ;
  2252.  'FINSI' ;
  2253. *****************    
  2254.  'SI' (ERR1 > ERRTOL) ;
  2255.         'MESSAGE' 'Problem 52 T2';
  2256.         'ERREUR' 5 ;
  2257.  'FINSI' ; 
  2258. **************************************************************
  2259.  'SI' (('ABS'(DRETYO2)) < DELTA) ;
  2260.     ERR1 = ('ABS' (DRETYO2N '-' DRETYO2)) ;
  2261.  'SINON' ;
  2262.     ERR1 = ('ABS' (DRETYO2N '-' DRETYO2)) '/' ('ABS' (DRETYO2)) ;
  2263.  'FINSI' ;
  2264. *****************
  2265.  'SI' (ERR1 > ERRTOL) ;
  2266.         'MESSAGE' 'Problem  53 T2';
  2267.         'ERREUR' 5 ;
  2268.  'FINSI' ;
  2269. **************************************************************
  2270.  'SI' (('ABS'(DYH2YO2)) < DELTA) ;
  2271.     ERR1 = ('ABS' (DYH2YO2N '-' DYH2YO2)) ;
  2272.  'SINON' ;  
  2273.     ERR1 = ('ABS' (DYH2YO2N '-' DYH2YO2)) '/' ('ABS' (DYH2YO2)) ;
  2274.  'FINSI' ;
  2275. **************************
  2276.  'SI' (ERR1 > ERRTOL) ;
  2277.       'MESSAGE' 'Problem 54 T2';
  2278.       'ERREUR' 5 ;
  2279.  'FINSI' ;
  2280. **************************************************************
  2281.  'SI' (('ABS'(DYO2YO2)) < DELTA) ;
  2282.     ERR1 = ('ABS' (DYO2YO2N '-' DYO2YO2)) ;
  2283.  'SINON' ;  
  2284.     ERR1 = ('ABS' (DYO2YO2N '-' DYO2YO2)) '/' ('ABS' (DYO2YO2)) ;
  2285.  'FINSI' ;
  2286. **************************
  2287.  'SI' (ERR1 > ERRTOL) ;
  2288.       'MESSAGE' 'Problem 55 T2';
  2289.       'ERREUR' 5 ;
  2290.  'FINSI' ;
  2291. **************************************************************
  2292.  'SI' (('ABS'(DYHOYO2)) < DELTA) ;
  2293.     ERR1 = ('ABS' (DYHOYO2N '-' DYHOYO2)) ;
  2294.  'SINON' ;  
  2295.     ERR1 = ('ABS' (DYHOYO2N '-' DYHOYO2)) '/' ('ABS' (DYHOYO2)) ;
  2296.  'FINSI' ;
  2297. **************************
  2298.  'SI' (ERR1 > ERRTOL) ;
  2299.       'MESSAGE' 'Problem 56 T2';
  2300.       'ERREUR' 5 ;
  2301.  'FINSI' ;
  2302. **************************************************************
  2303. *********  Eigth line !!!!!!!!!!!!!!!!!!!!!!   ***************
  2304. **************************************************************
  2305.  'SI' (('ABS'(DRYHO)) < DELTA) ;
  2306.     ERR1 = ('ABS' (DRYHON '-' DRYHO)) ;
  2307.  'SINON' ;
  2308.     ERR1 = ('ABS' (DRYHON '-' DRYHO)) '/' ('ABS' (DRYHO)) ;
  2309.  'FINSI' ;
  2310. *****************
  2311.  'SI' (ERR1 > ERRTOL) ;
  2312.         'MESSAGE' 'Problem 57 T2';
  2313.         'ERREUR' 5 ;
  2314.  'FINSI' ;
  2315. **************************************************************
  2316.  'SI' (('ABS'(DGXYHO)) < DELTA) ;
  2317.     ERR1 = ('ABS' (DGXYHON '-' DGXYHO)) ;
  2318.  'SINON' ;
  2319.     ERR1 = ('ABS' (DGXYHON '-' DGXYHO)) '/' ('ABS' (DGXYHO)) ;
  2320.  'FINSI' ;
  2321. *****************
  2322.  'SI' (ERR1 > ERRTOL) ;
  2323.         'MESSAGE' 'Problem 58 T2';
  2324.         'ERREUR' 5 ;
  2325.  'FINSI' ;
  2326. *************************************************************
  2327.  'SI' (('ABS'(DGYYHO)) < DELTA) ;
  2328.     ERR1 = ('ABS' (DGYYHON '-' DGYYHO)) ; 
  2329.  'SINON' ;
  2330.     ERR1 = ('ABS' (DGYYHON '-' DGYYHO)) '/' ('ABS' (DGYYHO)) ;
  2331.  'FINSI' ;
  2332. *****************    
  2333.  'SI' (ERR1 > ERRTOL) ;
  2334.         'MESSAGE' 'Problem 59 T2';
  2335.         'ERREUR' 5 ;
  2336.  'FINSI' ;
  2337. **************************************************************
  2338.  'SI' (('ABS'(DGZYHO)) < DELTA) ;
  2339.     ERR1 = ('ABS' (DGZYHON '-' DGZYHO)) ; 
  2340.  'SINON' ;
  2341.     ERR1 = ('ABS' (DGZYHON '-' DGZYHO)) '/' ('ABS' (DGZYHO)) ;
  2342.  'FINSI' ;
  2343. *****************    
  2344.  'SI' (ERR1 > ERRTOL) ;
  2345.         'MESSAGE' 'Problem 60 T2';
  2346.         'ERREUR' 5 ;
  2347.  'FINSI' ; 
  2348. **************************************************************
  2349.  'SI' (('ABS'(DRETYHO)) < DELTA) ;
  2350.     ERR1 = ('ABS' (DRETYHON '-' DRETYHO)) ;
  2351.  'SINON' ;
  2352.     ERR1 = ('ABS' (DRETYHON '-' DRETYHO)) '/' ('ABS' (DRETYHO)) ;
  2353.  'FINSI' ;
  2354. *****************
  2355.  'SI' (ERR1 > ERRTOL) ;
  2356.         'MESSAGE' 'Problem 61 T2';
  2357.         'ERREUR' 5 ;
  2358.  'FINSI' ;
  2359. **************************************************************
  2360.  'SI' (('ABS'(DYH2YHO)) < DELTA) ;
  2361.     ERR1 = ('ABS' (DYH2YHON '-' DYH2YHO)) ;
  2362.  'SINON' ;  
  2363.     ERR1 = ('ABS' (DYH2YHON '-' DYH2YHO)) '/' ('ABS' (DYH2YHO)) ;
  2364.  'FINSI' ;
  2365. **************************
  2366.  'SI' (ERR1 > ERRTOL) ;
  2367.       'MESSAGE' 'Problem 62 T2';
  2368.       'ERREUR' 5 ;
  2369.  'FINSI' ;
  2370. **************************************************************
  2371.  'SI' (('ABS'(DYO2YHO)) < DELTA) ;
  2372.     ERR1 = ('ABS' (DYO2YHON '-' DYO2YHO)) ;
  2373.  'SINON' ;  
  2374.     ERR1 = ('ABS' (DYO2YHON '-' DYO2YHO)) '/' ('ABS' (DYO2YHO)) ;
  2375.  'FINSI' ;
  2376. **************************
  2377.  'SI' (ERR1 > ERRTOL) ;
  2378.       'MESSAGE' 'Problem 63 T2';
  2379.       'ERREUR' 5 ;
  2380.  'FINSI' ;
  2381. **************************************************************
  2382.  'SI' (('ABS'(DYHOYHO)) < DELTA) ;
  2383.     ERR1 = ('ABS' (DYHOYHON '-' DYHOYHO)) ;
  2384.  'SINON' ;  
  2385.     ERR1 = ('ABS' (DYHOYHON '-' DYHOYHO)) '/' ('ABS' (DYHOYHO)) ;
  2386.  'FINSI' ;
  2387. **************************
  2388.  'SI' (ERR1 > ERRTOL) ;
  2389.       'MESSAGE' 'Problem 64 T2';
  2390.       'ERREUR' 5 ;
  2391.  'FINSI' ;
  2392.  
  2393. *****************************************************
  2394. *****************************************************
  2395. ******* TEST3 ***************************************
  2396. *****************************************************
  2397. *****************************************************
  2398. *
  2399. * On observe la variation du residu en $DOM27 . 'CENTRE'
  2400. * par rapport à une variation infinitésimal en
  2401. * $DOM27 . 'CENTRE' (NB : DOM27 est sur le bord!).
  2402. *
  2403.  
  2404.  
  2405. * Les grandeurs pour adimesionner les erreurs
  2406.  
  2407.   PCEN9 = ('DOMA '$DOM27 'CENTRE') 'POIN' 1 ;
  2408. ****************************************************************
  2409. ****************************************************************
  2410.   DRR = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMDEN ;
  2411.   DGXR = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMMOX ;
  2412.   DGYR = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMMOY ;
  2413.   DGZR = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMMOZ ;
  2414.   DRETR = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMRET ;
  2415.   DYH2R = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMYH2 ;
  2416.   DYO2R = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMYO2 ;
  2417.   DYHOR = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMDEN NOMYHO ;
  2418. ***************************************
  2419.   DRGX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMDEN ;
  2420.   DGXGX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMMOX ;
  2421.   DGYGX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMMOY ;
  2422.   DGZGX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMMOZ ;
  2423.   DRETGX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMRET ;
  2424.   DYH2GX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMYH2 ;
  2425.   DYO2GX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMYO2 ;
  2426.   DYHOGX = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOX NOMYHO ;
  2427. ***************************************
  2428.   DRGY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMDEN ;
  2429.   DGXGY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMMOX ;
  2430.   DGYGY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMMOY ;
  2431.   DGZGY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMMOZ ; 
  2432.   DRETGY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMRET ;
  2433.   DYH2GY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMYH2 ;
  2434.   DYO2GY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMYO2 ;
  2435.   DYHOGY = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOY NOMYHO ;
  2436. ***************************************
  2437.   DRGZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMDEN ;
  2438.   DGXGZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMMOX ;
  2439.   DGYGZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMMOY ;
  2440.   DGZGZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMMOZ ; 
  2441.   DRETGZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMRET ;
  2442.   DYH2GZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMYH2 ;
  2443.   DYO2GZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMYO2 ;
  2444.   DYHOGZ = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMMOZ NOMYHO ;    
  2445. ***************************************
  2446.   DRRET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMDEN ;
  2447.   DGXRET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMMOX ;
  2448.   DGYRET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMMOY ;
  2449.   DGZRET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMMOZ ;
  2450.   DRETRET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMRET ;
  2451.   DYH2RET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMYH2 ;
  2452.   DYO2RET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMYO2 ;
  2453.   DYHORET = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMRET NOMYHO ;
  2454. ***************************************
  2455.   DRYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMDEN ;
  2456.   DGXYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMMOX ;
  2457.   DGYYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMMOY ;
  2458.   DGZYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMMOZ ;
  2459.   DRETYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMRET ;
  2460.   DYH2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMYH2 ;
  2461.   DYO2YH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMYO2 ;
  2462.   DYHOYH2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYH2 NOMYHO ;
  2463. ***************************************
  2464.   DRYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMDEN ;
  2465.   DGXYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMMOX ;
  2466.   DGYYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMMOY ;
  2467.   DGZYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMMOZ ;
  2468.   DRETYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMRET ;
  2469.   DYH2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMYH2 ;
  2470.   DYO2YO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMYO2 ;
  2471.   DYHOYO2 = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYO2 NOMYHO ;
  2472. ****************************************
  2473.   DRYHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMDEN ;
  2474.   DGXYHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMMOX ;
  2475.   DGYYHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMMOY ;
  2476.   DGZYHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMMOZ ;
  2477.   DRETYHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMRET ;
  2478.   DYH2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMYH2 ;
  2479.   DYO2YHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMYO2 ;
  2480.   DYHOYHO = JAEXAC RMAT $DOMTOT LINCO PCEN27 PCEN9 NOMYHO NOMYHO ;
  2481.  
  2482.  
  2483.  
  2484. *****************************************************************
  2485. *****************************************************************
  2486. *****************************************************************  
  2487.   DRRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2488.       PCEN9 NOMDEN NOMDEN DELTA ;  
  2489.   DGXRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2490.       PCEN9 NOMDEN NOMMOX DELTA ;      
  2491.   DGYRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2492.       PCEN9 NOMDEN NOMMOY DELTA ;
  2493.   DGZRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2494.       PCEN9 NOMDEN NOMMOZ DELTA ;    
  2495.   DRETRN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2496.       PCEN9 NOMDEN NOMRET DELTA ;
  2497.   DYH2RN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2498.       PCEN9 NOMDEN NOMYH2 DELTA ;
  2499.   DYO2RN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2500.       PCEN9 NOMDEN NOMYO2 DELTA ;
  2501.   DYHORN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2502.       PCEN9 NOMDEN NOMYHO DELTA ;
  2503. **************************************
  2504.   DRGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2505.       PCEN9 NOMMOX NOMDEN DELTA ;  
  2506.   DGXGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2507.       PCEN9 NOMMOX NOMMOX DELTA ;      
  2508.   DGYGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2509.       PCEN9 NOMMOX NOMMOY DELTA ;
  2510.   DGZGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2511.       PCEN9 NOMMOX NOMMOZ DELTA ;    
  2512.   DRETGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2513.       PCEN9 NOMMOX NOMRET DELTA ;
  2514.   DYH2GXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2515.       PCEN9 NOMMOX NOMYH2 DELTA ;
  2516.   DYO2GXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2517.       PCEN9 NOMMOX NOMYO2 DELTA ;
  2518.   DYHOGXN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2519.       PCEN9 NOMMOX NOMYHO DELTA ;
  2520. **************************************
  2521.   DRGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2522.       PCEN9 NOMMOY NOMDEN DELTA ;  
  2523.   DGXGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2524.       PCEN9 NOMMOY NOMMOX DELTA ;      
  2525.   DGYGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2526.       PCEN9 NOMMOY NOMMOY DELTA ;
  2527.   DGZGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2528.       PCEN9 NOMMOY NOMMOZ DELTA ;    
  2529.   DRETGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2530.       PCEN9 NOMMOY NOMRET DELTA ;
  2531.   DYH2GYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2532.       PCEN9 NOMMOY NOMYH2 DELTA ;
  2533.   DYO2GYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2534.       PCEN9 NOMMOY NOMYO2 DELTA ;
  2535.   DYHOGYN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2536.       PCEN9 NOMMOY NOMYHO DELTA ;
  2537. **************************************
  2538.   DRGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2539.       PCEN9 NOMMOZ NOMDEN DELTA ;  
  2540.   DGXGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2541.       PCEN9 NOMMOZ NOMMOX DELTA ;      
  2542.   DGYGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2543.       PCEN9 NOMMOZ NOMMOY DELTA ;
  2544.   DGZGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2545.       PCEN9 NOMMOZ NOMMOZ DELTA ;    
  2546.   DRETGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2547.       PCEN9 NOMMOZ NOMRET DELTA ;
  2548.   DYH2GZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2549.       PCEN9 NOMMOZ NOMYH2 DELTA ;
  2550.   DYO2GZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2551.       PCEN9 NOMMOZ NOMYO2 DELTA ;
  2552.   DYHOGZN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2553.       PCEN9 NOMMOZ NOMYHO DELTA ;      
  2554. **************************************
  2555.   DRRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2556.       PCEN9 NOMRET NOMDEN DELTA ;  
  2557.   DGXRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2558.       PCEN9 NOMRET NOMMOX DELTA ;      
  2559.   DGYRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2560.       PCEN9 NOMRET NOMMOY DELTA ;
  2561.   DGZRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2562.       PCEN9 NOMRET NOMMOZ DELTA ;    
  2563.   DRETRETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2564.       PCEN9 NOMRET NOMRET DELTA ;
  2565.   DYH2RETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2566.       PCEN9 NOMRET NOMYH2 DELTA ;
  2567.   DYO2RETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2568.       PCEN9 NOMRET NOMYO2 DELTA ;
  2569.   DYHORETN = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2570.       PCEN9 NOMRET NOMYHO DELTA ;
  2571. **************************************
  2572.   DRYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2573.       PCEN9 NOMYH2 NOMDEN DELTA ;  
  2574.   DGXYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2575.       PCEN9 NOMYH2 NOMMOX DELTA ;      
  2576.   DGYYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2577.       PCEN9 NOMYH2 NOMMOY DELTA ;
  2578.   DGZYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2579.       PCEN9 NOMYH2 NOMMOZ DELTA ;    
  2580.   DRETYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2581.       PCEN9 NOMYH2 NOMRET DELTA ;
  2582.   DYH2YH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2583.       PCEN9 NOMYH2 NOMYH2 DELTA ;
  2584.   DYO2YH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2585.       PCEN9 NOMYH2 NOMYO2 DELTA ;
  2586.   DYHOYH2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2587.       PCEN9 NOMYH2 NOMYHO DELTA ;
  2588. ***************************************
  2589.   DRYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2590.       PCEN9 NOMYO2 NOMDEN DELTA ;  
  2591.   DGXYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2592.       PCEN9 NOMYO2 NOMMOX DELTA ;      
  2593.   DGYYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2594.       PCEN9 NOMYO2 NOMMOY DELTA ;
  2595.   DGZYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2596.       PCEN9 NOMYO2 NOMMOZ DELTA ;    
  2597.   DRETYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2598.       PCEN9 NOMYO2 NOMRET DELTA ;
  2599.   DYH2YO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2600.       PCEN9 NOMYO2 NOMYH2 DELTA ;
  2601.   DYO2YO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2602.       PCEN9 NOMYO2 NOMYO2 DELTA ;
  2603.   DYHOYO2N = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2604.       PCEN9 NOMYO2 NOMYHO DELTA ;
  2605. ***************************************
  2606.   DRYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2607.       PCEN9 NOMYHO NOMDEN DELTA ;  
  2608.   DGXYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2609.       PCEN9 NOMYHO NOMMOX DELTA ;      
  2610.   DGYYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2611.       PCEN9 NOMYHO NOMMOY DELTA ;
  2612.   DGZYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2613.       PCEN9 NOMYHO NOMMOZ DELTA ;    
  2614.   DRETYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2615.       PCEN9 NOMYHO NOMRET DELTA ;
  2616.   DYH2YHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2617.       PCEN9 NOMYHO NOMYH2 DELTA ;
  2618.   DYO2YHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2619.       PCEN9 NOMYHO NOMYO2 DELTA ;
  2620.   DYHOYHON = JACNUM $DOMTOT METO RN GN RETN RYN GAMN PGAS LINCO PCEN27
  2621.       PCEN9 NOMYHO NOMYHO DELTA ;   
  2622.  
  2623. * Test des comparaisons jacobien exact-jacobien numerique
  2624. * Test des comparaisons jacobien exact-jacobien numerique
  2625. **************************************************************
  2626. 'SI' (('ABS'(DRR)) < DELTA) ;
  2627.     ERR1 = ('ABS' (DRRN '-' DRR)) ;
  2628.  'SINON' ;  
  2629.     ERR1 = ('ABS' (DRRN '-' DRR)) '/' ('ABS' (DRR)) ;
  2630.  'FINSI' ; 
  2631. ***********************
  2632.  'SI' (ERR1 > ERRTOL) ;
  2633.       'MESSAGE' 'Problem 1 T3';
  2634.       'ERREUR' 5 ;
  2635.  'FINSI' ;
  2636. **************************************************************
  2637.  'SI' (('ABS'(DGXR)) < DELTA) ;
  2638.     ERR1 = ('ABS' (DGXRN '-' DGXR)) ;
  2639.  'SINON' ;  
  2640.     ERR1 = ('ABS' (DGXRN '-' DGXR)) '/' ('ABS' (DGXR)) ;
  2641.  'FINSI' ; 
  2642. *********************** 
  2643.  'SI' (ERR1 > ERRTOL) ;
  2644.       'MESSAGE' 'Problem 2 T3';
  2645.       'ERREUR' 5 ;
  2646.  'FINSI' ;
  2647. **************************************************************
  2648.  'SI' (('ABS'(DGYR)) < DELTA) ;
  2649.     ERR1 = ('ABS' (DGYRN '-' DGYR)) ;
  2650.  'SINON' ;  
  2651.     ERR1 = ('ABS' (DGYRN '-' DGYR)) '/' ('ABS' (DGYR)) ;
  2652.  'FINSI' ; 
  2653. ************************
  2654.  'SI' (ERR1 > ERRTOL) ;
  2655.       'MESSAGE' 'Problem 3 T3';
  2656.        'ERREUR' 5 ;
  2657.  'FINSI' ;
  2658. **************************************************************
  2659.  'SI' (('ABS'(DGZR)) < DELTA) ;
  2660.     ERR1 = ('ABS' (DGZRN '-' DGZR)) ;
  2661.  'SINON' ;  
  2662.     ERR1 = ('ABS' (DGZRN '-' DGZR)) '/' ('ABS' (DGZR)) ;
  2663.  'FINSI' ;
  2664. ************************
  2665.  'SI' (ERR1 > ERRTOL) ;
  2666.       'MESSAGE' 'Problem 4 T3';
  2667.        'ERREUR' 5 ;
  2668.  'FINSI' ;
  2669. ***************************************************************
  2670.  'SI' (('ABS'(DRETR)) < DELTA) ;
  2671.     ERR1 = ('ABS' (DRETRN '-' DRETR)) ;
  2672.  'SINON' ;  
  2673.     ERR1 = ('ABS' (DRETRN '-' DRETR)) '/' ('ABS' (DRETR)) ;
  2674.  'FINSI' ;
  2675. **************************
  2676.  'SI' (ERR1 > ERRTOL) ;
  2677.       'MESSAGE' 'Problem 5 T3';
  2678.       'ERREUR' 5 ;
  2679.  'FINSI' ;
  2680. **************************************************************
  2681. **************************************************************
  2682.  'SI' (('ABS'(DYH2R)) < DELTA) ;
  2683.     ERR1 = ('ABS' (DYH2RN '-' DYH2R)) ;
  2684.  'SINON' ;  
  2685.     ERR1 = ('ABS' (DYH2RN '-' DYH2R)) '/' ('ABS' (DYH2R)) ;
  2686.  'FINSI' ;
  2687. **************************
  2688.  'SI' (ERR1 > ERRTOL) ;
  2689.       'MESSAGE' 'Problem 6 T3';
  2690.       'ERREUR' 5 ;
  2691.  'FINSI' ;
  2692. **************************************************************
  2693.    'SI' (('ABS'(DYO2R)) < DELTA) ;
  2694.     ERR1 = ('ABS' (DYO2RN '-' DYO2R)) ;
  2695.  'SINON' ;  
  2696.     ERR1 = ('ABS' (DYO2RN '-' DYO2R)) '/' ('ABS' (DYO2R)) ;
  2697.  'FINSI' ;
  2698. **************************
  2699.  'SI' (ERR1 > ERRTOL) ;
  2700.       'MESSAGE' 'Problem 7 T3';
  2701.       'ERREUR' 5 ;
  2702.  'FINSI' ;
  2703. **************************************************************
  2704.   'SI' (('ABS'(DYHOR)) < DELTA) ;
  2705.     ERR1 = ('ABS' (DYHORN '-' DYHOR)) ;
  2706.  'SINON' ;  
  2707.     ERR1 = ('ABS' (DYHORN '-' DYHOR)) '/' ('ABS' (DYHOR)) ;
  2708.  'FINSI' ;
  2709. **************************
  2710.  'SI' (ERR1 > ERRTOL) ;
  2711.       'MESSAGE' 'Problem 8 T3';
  2712.       'ERREUR' 5 ;
  2713.  'FINSI' ; 
  2714. **************************************************************
  2715. **************************************************************
  2716. **************************************************************
  2717. *****    Second line !!!!!!!!!!!!!!!!!!!!!!!    **************
  2718. **************************************************************
  2719. ************************************************************** 
  2720.  'SI' (('ABS'(DRGX)) < DELTA) ;
  2721.     ERR1 = ('ABS' (DRGXN '-' DRGX)) ;
  2722.  'SINON' ;   
  2723.     ERR1 = ('ABS' (DRGXN '-' DRGX)) '/' ('ABS' (DRGX))  ;
  2724.  'FINSI' ;
  2725. *************************** 
  2726.  'SI' (ERR1 > ERRTOL) ;
  2727.       'MESSAGE' 'Problem 9 T3';
  2728.       'ERREUR' 5 ;
  2729.  'FINSI' ;
  2730. ***************************************************************
  2731.  'SI' (('ABS'(DGXGX)) < DELTA) ;
  2732.    ERR1 = ('ABS' (DGXGXN '-' DGXGX)) ;
  2733.  'SINON' ;
  2734.    ERR1 = ('ABS' (DGXGXN '-' DGXGX)) '/' ('ABS' (DGXGX)) ;
  2735.  'FINSI' ;  
  2736. *****************************   
  2737.  'SI' (ERR1 > ERRTOL) ;
  2738.         'MESSAGE' 'Problem 10 T3';
  2739.         'ERREUR' 5 ;
  2740.  'FINSI' ;
  2741. **************************************************************
  2742.  'SI' (('ABS'(DGYGX)) < DELTA) ;
  2743.     ERR1 = ('ABS' (DGYGXN '-' DGYGX)) ;
  2744.  'SINON' ;
  2745.     ERR1 = ('ABS' (DGYGXN '-' DGYGX)) '/' ('ABS' (DGYGX)) ;
  2746.  'FINSI' ;
  2747. ******************************     
  2748.  'SI' (ERR1 > ERRTOL) ;
  2749.         'MESSAGE' 'Problem 11 T3';
  2750.         'ERREUR' 5 ;
  2751.  'FINSI' ;
  2752. **************************************************************
  2753.  'SI' (('ABS'(DGZGX)) < DELTA) ;
  2754.     ERR1 = ('ABS' (DGZGXN '-' DGZGX)) ;
  2755.  'SINON' ;
  2756.     ERR1 = ('ABS' (DGZGXN '-' DGZGX)) '/' ('ABS' (DGZGX)) ;
  2757.  'FINSI' ;
  2758. ******************************    
  2759.  'SI' (ERR1 > ERRTOL) ;
  2760.         'MESSAGE' 'Problem 12 T3';
  2761.         'ERREUR' 5 ;
  2762.  'FINSI' ;
  2763. **************************************************************
  2764.  'SI' (('ABS'(DRETGX)) < DELTA) ;
  2765.     ERR1 = ('ABS' (DRETGXN '-' DRETGX)) ;
  2766.  'SINON' ;
  2767.     ERR1 = ('ABS' (DRETGXN '-' DRETGX)) '/' ('ABS' (DRETGX)) ;
  2768.  'FINSI' ;
  2769. *******************************    
  2770.  'SI' (ERR1 > ERRTOL) ;
  2771.         'MESSAGE' 'Problem 13 T3';
  2772.         'ERREUR' 5 ;
  2773.  'FINSI' ;
  2774. **************************************************************
  2775.  'SI' (('ABS'(DYH2GX)) < DELTA) ;
  2776.     ERR1 = ('ABS' (DYH2GXN '-' DYH2GX)) ;
  2777.  'SINON' ;  
  2778.     ERR1 = ('ABS' (DYH2GXN '-' DYH2GX)) '/' ('ABS' (DYH2GX)) ;
  2779.  'FINSI' ;
  2780. **************************
  2781.  'SI' (ERR1 > ERRTOL) ;
  2782.       'MESSAGE' 'Problem 14 T3';
  2783.       'ERREUR' 5 ;
  2784.  'FINSI' ;
  2785. **************************************************************
  2786.  'SI' (('ABS'(DYO2GX)) < DELTA) ;
  2787.     ERR1 = ('ABS' (DYO2GXN '-' DYO2GX)) ;
  2788.  'SINON' ;  
  2789.     ERR1 = ('ABS' (DYO2GXN '-' DYO2GX)) '/' ('ABS' (DYO2GX)) ;
  2790.  'FINSI' ;
  2791. **************************
  2792.  'SI' (ERR1 > ERRTOL) ;
  2793.       'MESSAGE' 'Problem 15 T3';
  2794.       'ERREUR' 5 ;
  2795.  'FINSI' ;
  2796. **************************************************************
  2797.  'SI' (('ABS'(DYHOGX)) < DELTA) ;
  2798.     ERR1 = ('ABS' (DYHOGXN '-' DYHOGX)) ;
  2799.  'SINON' ;  
  2800.     ERR1 = ('ABS' (DYHOGXN '-' DYHOGX)) '/' ('ABS' (DYHOGX)) ;
  2801.  'FINSI' ;
  2802. **************************
  2803.  'SI' (ERR1 > ERRTOL) ;
  2804.       'MESSAGE' 'Problem 16 T3';
  2805.       'ERREUR' 5 ;
  2806.  'FINSI' ;
  2807. **************************************************************
  2808. *******        Third line!!!!!!!!!!!!!!!!      ***************
  2809. **************************************************************
  2810.  'SI' (('ABS'(DRGY)) < DELTA) ;
  2811.     ERR1 = ('ABS' (DRGYN '-' DRGY)) ;
  2812.  'SINON' ;
  2813.     ERR1 = ('ABS' (DRGYN '-' DRGY)) '/' ('ABS' (DRGY)) ;
  2814.  'FINSI' ;
  2815. *******************************       
  2816.  'SI' (ERR1 > ERRTOL) ;
  2817.         'MESSAGE' 'Problem 17 T3';
  2818.         'ERREUR' 5 ;
  2819.  'FINSI' ;
  2820. **************************************************************
  2821.  'SI' (('ABS'(DGXGY)) < DELTA) ;
  2822.     ERR1 = ('ABS' (DGXGYN '-' DGXGY)) ;
  2823.  'SINON' ;
  2824.     ERR1 = ('ABS' (DGXGYN '-' DGXGY))  '/' ('ABS' (DGXGY)) ;
  2825.  'FINSI' ;
  2826. ********************************      
  2827.  'SI' (ERR1 > ERRTOL) ;
  2828.         'MESSAGE' 'Problem 18 T3';
  2829.         'ERREUR' 5 ;
  2830.  'FINSI' ;
  2831. **************************************************************
  2832.  'SI' (('ABS'(DGYGY)) < DELTA) ;
  2833.     ERR1 = ('ABS' (DGYGYN '-' DGYGY)) ;
  2834.  'SINON' ;
  2835.     ERR1 = ('ABS' (DGYGYN '-' DGYGY)) '/' ('ABS' (DGYGY)) ;
  2836.  'FINSI' ;
  2837. ***************
  2838.  'SI' (ERR1 > ERRTOL) ;
  2839.         'MESSAGE' 'Problem 19 T3';
  2840.         'ERREUR' 5 ;
  2841.  'FINSI' ;
  2842. **************************************************************
  2843.  'SI' (('ABS'(DGZGY)) < DELTA) ;
  2844.     ERR1 = ('ABS' (DGZGYN '-' DGZGY)) ;
  2845.  'SINON' ;
  2846.     ERR1 = ('ABS' (DGZGYN '-' DGZGY)) '/' ('ABS' (DGZGY)) ;
  2847.  'FINSI' ;
  2848. ****************        
  2849.  'SI' (ERR1 > ERRTOL) ;
  2850.         'MESSAGE' 'Problem 20 T3';
  2851.         'ERREUR' 5 ;
  2852.  'FINSI' ;
  2853. **************************************************************
  2854.  'SI' (('ABS'(DRETGY)) < DELTA) ;
  2855.     ERR1 = ('ABS' (DRETGYN '-' DRETGY)) ;
  2856.  'SINON' ;
  2857.     ERR1 = ('ABS' (DRETGYN '-' DRETGY)) '/' ('ABS' (DRETGY)) ;
  2858.  'FINSI' ;
  2859. ****************     
  2860.  'SI' (ERR1 > ERRTOL) ;
  2861.         'MESSAGE' 'Problem 21 T3';
  2862.         'ERREUR' 5 ;
  2863.  'FINSI' ;
  2864. **************************************************************
  2865.  'SI' (('ABS'(DYH2GY)) < DELTA) ;
  2866.     ERR1 = ('ABS' (DYH2GYN '-' DYH2GY)) ;
  2867.  'SINON' ;  
  2868.     ERR1 = ('ABS' (DYH2GYN '-' DYH2GY)) '/' ('ABS' (DYH2GY)) ;
  2869.  'FINSI' ;
  2870. **************************
  2871.  'SI' (ERR1 > ERRTOL) ;
  2872.       'MESSAGE' 'Problem 22 T3';
  2873.       'ERREUR' 5 ;
  2874.  'FINSI' ;
  2875. **************************************************************
  2876.  'SI' (('ABS'(DYO2GY)) < DELTA) ;
  2877.     ERR1 = ('ABS' (DYO2GYN '-' DYO2GY)) ;
  2878.  'SINON' ;  
  2879.     ERR1 = ('ABS' (DYO2GYN '-' DYO2GY)) '/' ('ABS' (DYO2GY)) ;
  2880.  'FINSI' ;
  2881. **************************
  2882.  'SI' (ERR1 > ERRTOL) ;
  2883.       'MESSAGE' 'Problem 23 T3';
  2884.       'ERREUR' 5 ;
  2885.  'FINSI' ;
  2886. **************************************************************
  2887.  'SI' (('ABS'(DYHOGY)) < DELTA) ;
  2888.     ERR1 = ('ABS' (DYHOGYN '-' DYHOGY)) ;
  2889.  'SINON' ;  
  2890.     ERR1 = ('ABS' (DYHOGYN '-' DYHOGY)) '/' ('ABS' (DYHOGY)) ;
  2891.  'FINSI' ;
  2892. **************************
  2893.  'SI' (ERR1 > ERRTOL) ;
  2894.       'MESSAGE' 'Problem 24 T3';
  2895.       'ERREUR' 5 ;
  2896.  'FINSI' ;
  2897. **************************************************************
  2898. *******     Fouth line!!!!!!!!!!!!!!!!!!!   ******************
  2899. **************************************************************
  2900.  'SI' (('ABS'(DRGZ)) < DELTA) ;
  2901.     ERR1 = ('ABS' (DRGZN '-' DRGZ))  ;
  2902.  'SINON' ;
  2903.     ERR1 = ('ABS' (DRGZN '-' DRGZ)) '/' ('ABS' (DRGZ))  ;
  2904.  'FINSI' ;
  2905. ***************      
  2906.  'SI' (ERR1 > ERRTOL) ;
  2907.         'MESSAGE' 'Problem 25 T3';
  2908.         'ERREUR' 5 ;
  2909.  'FINSI' ;
  2910. **************************************************************
  2911.  'SI' (('ABS'(DGXGZ)) < DELTA) ;
  2912.     ERR1 = ('ABS' (DGXGZN '-' DGXGZ)) ;
  2913.  'SINON' ;
  2914.     ERR1 = ('ABS' (DGXGZN '-' DGXGZ))  '/' ('ABS' (DGXGZ)) ;
  2915.  'FINSI' ;
  2916. *****************       
  2917.  'SI' (ERR1 > ERRTOL) ;
  2918.         'MESSAGE' 'Problem 26 T3';
  2919.         'ERREUR' 5 ;
  2920.  'FINSI' ;
  2921. **************************************************************
  2922.  'SI' (('ABS'(DGYGZ)) < DELTA) ;
  2923.     ERR1 = ('ABS' (DGYGZN '-' DGYGZ)) ;
  2924.  'SINON' ;
  2925.     ERR1 = ('ABS' (DGYGZN '-' DGYGZ)) '/' ('ABS' (DGYGZ)) ;
  2926.  'FINSI' ;
  2927. *****************       
  2928.  'SI' (ERR1 > ERRTOL) ;
  2929.         'MESSAGE' 'Problem 27 T3';
  2930.         'ERREUR' 5 ;
  2931.  'FINSI' ;
  2932. **************************************************************
  2933.  'SI' (('ABS'(DGZGZ)) < DELTA) ;
  2934.     ERR1 = ('ABS' (DGZGZN '-' DGZGZ)) ;
  2935.  'SINON' ;
  2936.     ERR1 = ('ABS' (DGZGZN '-' DGZGZ)) '/' ('ABS' (DGZGZ)) ;
  2937.  'FINSI' ;
  2938. *****************     
  2939.  'SI' (ERR1 > ERRTOL) ;
  2940.         'MESSAGE' 'Problem 28 T3';
  2941.         'ERREUR' 5 ;
  2942.  'FINSI' ;
  2943. **************************************************************
  2944.  'SI' (('ABS'(DRETGZ)) < DELTA) ;
  2945.     ERR1 = ('ABS' (DRETGZN '-' DRETGZ)) ;
  2946.  'SINON' ;
  2947.     ERR1 = ('ABS' (DRETGZN '-' DRETGZ)) '/' ('ABS' (DRETGZ)) ;
  2948.  'FINSI' ;     
  2949. *****************
  2950.  'SI' (ERR1 > ERRTOL) ;
  2951.         'MESSAGE' 'Problem 29 T3';
  2952.         'ERREUR' 5 ;
  2953.  'FINSI' ;
  2954. **************************************************************
  2955.  'SI' (('ABS'(DYH2GZ)) < DELTA) ;
  2956.     ERR1 = ('ABS' (DYH2GZN '-' DYH2GZ)) ;
  2957.  'SINON' ;  
  2958.     ERR1 = ('ABS' (DYH2GZN '-' DYH2GZ)) '/' ('ABS' (DYH2GZ)) ;
  2959.  'FINSI' ;
  2960. **************************
  2961.  'SI' (ERR1 > ERRTOL) ;
  2962.       'MESSAGE' 'Problem 30 T3';
  2963.       'ERREUR' 5 ;
  2964.  'FINSI' ;
  2965. **************************************************************
  2966.  'SI' (('ABS'(DYO2GZ)) < DELTA) ;
  2967.     ERR1 = ('ABS' (DYO2GZN '-' DYO2GZ)) ;
  2968.  'SINON' ;  
  2969.     ERR1 = ('ABS' (DYO2GZN '-' DYO2GZ)) '/' ('ABS' (DYO2GZ)) ;
  2970.  'FINSI' ;
  2971. **************************
  2972.  'SI' (ERR1 > ERRTOL) ;
  2973.       'MESSAGE' 'Problem 31 T3';
  2974.       'ERREUR' 5 ;
  2975.  'FINSI' ;
  2976. **************************************************************
  2977.  'SI' (('ABS'(DYHOGZ)) < DELTA) ;
  2978.     ERR1 = ('ABS' (DYHOGZN '-' DYHOGZ)) ;
  2979.  'SINON' ;  
  2980.     ERR1 = ('ABS' (DYHOGZN '-' DYHOGZ)) '/' ('ABS' (DYHOGZ)) ;
  2981.  'FINSI' ;
  2982. **************************
  2983.  'SI' (ERR1 > ERRTOL) ;
  2984.       'MESSAGE' 'Problem 32 T3';
  2985.       'ERREUR' 5 ;
  2986.  'FINSI' ;
  2987. **************************************************************
  2988. *******     Fith line!!!!!!!!!!!!!!!!!!!   ******************
  2989. **************************************************************
  2990.  'SI' (('ABS'(DRRET)) < DELTA) ;
  2991.     ERR1 = ('ABS' (DRRETN '-' DRRET)) ;
  2992.  'SINON' ;
  2993.     ERR1 = ('ABS' (DRRETN '-' DRRET)) '/' ('ABS' (DRRET)) ;
  2994.  'FINSI' ;
  2995. *****************       
  2996.  'SI' (ERR1 > ERRTOL) ;
  2997.         'MESSAGE' 'Problem 33 T3';
  2998.         'ERREUR' 5 ;
  2999.  'FINSI' ;
  3000. ***************************************************************
  3001.  'SI' (('ABS'(DGXRET)) < DELTA) ;
  3002.     ERR1 = ('ABS' (DGXRETN '-' DGXRET)) ;
  3003.  'SINON' ;
  3004.     ERR1 = ('ABS' (DGXRETN '-' DGXRET)) '/' ('ABS' (DGXRET)) ;
  3005.  'FINSI' ;
  3006. *****************    
  3007.  'SI' (ERR1 > ERRTOL) ;
  3008.         'MESSAGE' 'Problem 34 T3';
  3009.         'ERREUR' 5 ;
  3010.  'FINSI' ;
  3011. ****************************************************************
  3012.  'SI' (('ABS'(DGYRET)) < DELTA) ;
  3013.     ERR1 = ('ABS' (DGYRETN '-' DGYRET)) ; 
  3014.  'SINON' ;
  3015.     ERR1 = ('ABS' (DGYRETN '-' DGYRET)) '/' ('ABS' (DGYRET)) ;
  3016.  'FINSI' ;
  3017. *****************    
  3018.  'SI' (ERR1 > ERRTOL) ;
  3019.         'MESSAGE' 'Problem 35 T3';
  3020.         'ERREUR' 5 ;
  3021.  'FINSI' ;
  3022. ****************************************************************
  3023.  'SI' (('ABS'(DGZRET)) < DELTA) ;
  3024.     ERR1 = ('ABS' (DGZRETN '-' DGZRET)) ;
  3025.  'SINON' ;    
  3026.     ERR1 = ('ABS' (DGZRETN '-' DGZRET)) '/' ('ABS' (DGZRET)) ;
  3027.  'FINSI' ;
  3028. *****************    
  3029.  'SI' (ERR1 > ERRTOL) ;
  3030.         'MESSAGE' 'Problem 36 T3';
  3031.         'ERREUR' 5 ;
  3032.  'FINSI' ;
  3033. ****************************************************************
  3034.  'SI' (('ABS'(DRETRET)) < DELTA) ;
  3035.     ERR1 = ('ABS' (DRETRETN '-' DRETRET)) ;
  3036.  'SINON' ;
  3037.     ERR1 = ('ABS' (DRETRETN '-' DRETRET)) '/' ('ABS' (DRETRET)) ;
  3038.  'FINSI' ;
  3039. *****************       
  3040.  'SI' (ERR1 > ERRTOL) ;
  3041.         'MESSAGE' 'Problem 37 T3';
  3042.         'ERREUR' 5 ;
  3043.  'FINSI' ;
  3044.  'SI' (('ABS'(DYH2RET)) < DELTA) ;
  3045.     ERR1 = ('ABS' (DYH2RETN '-' DYH2RET)) ;
  3046.  'SINON' ;  
  3047.     ERR1 = ('ABS' (DYH2RETN '-' DYH2RET)) '/' ('ABS' (DYH2RET)) ;
  3048.  'FINSI' ;
  3049. **************************
  3050.  'SI' (ERR1 > ERRTOL) ;
  3051.       'MESSAGE' 'Problem 38 T3';
  3052.       'ERREUR' 5 ;
  3053.  'FINSI' ;
  3054. **************************************************************
  3055.  'SI' (('ABS'(DYO2RET)) < DELTA) ;
  3056.     ERR1 = ('ABS' (DYO2RETN '-' DYO2RET)) ;
  3057.  'SINON' ;  
  3058.     ERR1 = ('ABS' (DYO2RETN '-' DYO2RET)) '/' ('ABS' (DYO2RET)) ;
  3059.  'FINSI' ;
  3060. **************************
  3061.  'SI' (ERR1 > ERRTOL) ;
  3062.       'MESSAGE' 'Problem 39 T3';
  3063.       'ERREUR' 5 ;
  3064.  'FINSI' ;
  3065. **************************************************************
  3066.  'SI' (('ABS'(DYHORET)) < DELTA) ;
  3067.     ERR1 = ('ABS' (DYHORETN '-' DYHORET)) ;
  3068.  'SINON' ;  
  3069.     ERR1 = ('ABS' (DYHORETN '-' DYHORET)) '/' ('ABS' (DYHORET)) ;
  3070.  'FINSI' ;
  3071. **************************
  3072.  'SI' (ERR1 > ERRTOL) ;
  3073.       'MESSAGE' 'Problem 40 T3';
  3074.       'ERREUR' 5 ;
  3075.  'FINSI' ;
  3076. **************************************************************
  3077. *******     Sixth line!!!!!!!!!!!!!!!!!!!   ******************
  3078. **************************************************************
  3079.  'SI' (('ABS'(DRYH2)) < DELTA) ;
  3080.     ERR1 = ('ABS' (DRYH2N '-' DRYH2)) ;
  3081.  'SINON' ;
  3082.     ERR1 = ('ABS' (DRYH2N '-' DRYH2)) '/' ('ABS' (DRYH2)) ;
  3083.  'FINSI' ;
  3084. *****************
  3085.  'SI' (ERR1 > ERRTOL) ;
  3086.         'MESSAGE' 'Problem 41 T3';
  3087.         'ERREUR' 5 ;
  3088.  'FINSI' ;
  3089. **************************************************************
  3090.  'SI' (('ABS'(DGXYH2)) < DELTA) ;
  3091.     ERR1 = ('ABS' (DGXYH2N '-' DGXYH2)) ;
  3092.  'SINON' ;
  3093.     ERR1 = ('ABS' (DGXYH2N '-' DGXYH2)) '/' ('ABS' (DGXYH2)) ;
  3094.  'FINSI' ;
  3095. *****************
  3096.  'SI' (ERR1 > ERRTOL) ;
  3097.         'MESSAGE' 'Problem 42 T3';
  3098.         'ERREUR' 5 ;
  3099.  'FINSI' ;
  3100. *************************************************************
  3101.  'SI' (('ABS'(DGYYH2)) < DELTA) ;
  3102.     ERR1 = ('ABS' (DGYYH2N '-' DGYYH2)) ; 
  3103.  'SINON' ;
  3104.     ERR1 = ('ABS' (DGYYH2N '-' DGYYH2)) '/' ('ABS' (DGYYH2)) ;
  3105.  'FINSI' ;
  3106. *****************    
  3107.  'SI' (ERR1 > ERRTOL) ;
  3108.         'MESSAGE' 'Problem 43 T3';
  3109.         'ERREUR' 5 ;
  3110.  'FINSI' ;
  3111. **************************************************************
  3112.  'SI' (('ABS'(DGZYH2)) < DELTA) ;
  3113.     ERR1 = ('ABS' (DGZYH2N '-' DGZYH2)) ; 
  3114.  'SINON' ;
  3115.     ERR1 = ('ABS' (DGZYH2N '-' DGZYH2)) '/' ('ABS' (DGZYH2)) ;
  3116.  'FINSI' ;
  3117. *****************    
  3118.  'SI' (ERR1 > ERRTOL) ;
  3119.         'MESSAGE' 'Problem 44 T3';
  3120.         'ERREUR' 5 ;
  3121.  'FINSI' ; 
  3122. **************************************************************
  3123.  'SI' (('ABS'(DRETYH2)) < DELTA) ;
  3124.     ERR1 = ('ABS' (DRETYH2N '-' DRETYH2)) ;
  3125.  'SINON' ;
  3126.     ERR1 = ('ABS' (DRETYH2N '-' DRETYH2)) '/' ('ABS' (DRETYH2)) ;
  3127.  'FINSI' ;
  3128. *****************
  3129.  'SI' (ERR1 > ERRTOL) ;
  3130.         'MESSAGE' 'Problem 45 T3';
  3131.         'ERREUR' 5 ;
  3132.  'FINSI' ;
  3133. **************************************************************
  3134.  'SI' (('ABS'(DYH2YH2)) < DELTA) ;
  3135.     ERR1 = ('ABS' (DYH2YH2N '-' DYH2YH2)) ;
  3136.  'SINON' ;  
  3137.     ERR1 = ('ABS' (DYH2YH2N '-' DYH2YH2)) '/' ('ABS' (DYH2YH2)) ;
  3138.  'FINSI' ;
  3139. **************************
  3140.  'SI' (ERR1 > ERRTOL) ;
  3141.       'MESSAGE' 'Problem 46 T3';
  3142.       'ERREUR' 5 ;
  3143.  'FINSI' ;
  3144. **************************************************************
  3145.  'SI' (('ABS'(DYO2YH2)) < DELTA) ;
  3146.     ERR1 = ('ABS' (DYO2YH2N '-' DYO2YH2)) ;
  3147.  'SINON' ;  
  3148.     ERR1 = ('ABS' (DYO2YH2N '-' DYO2YH2)) '/' ('ABS' (DYO2YH2)) ;
  3149.  'FINSI' ;
  3150. **************************
  3151.  'SI' (ERR1 > ERRTOL) ;
  3152.       'MESSAGE' 'Problem 47 T3';
  3153.       'ERREUR' 5 ;
  3154.  'FINSI' ;
  3155. **************************************************************
  3156.  'SI' (('ABS'(DYHOYH2)) < DELTA) ;
  3157.     ERR1 = ('ABS' (DYHOYH2N '-' DYHOYH2)) ;
  3158.  'SINON' ;  
  3159.     ERR1 = ('ABS' (DYHOYH2N '-' DYHOYH2)) '/' ('ABS' (DYHOYH2)) ;
  3160.  'FINSI' ;
  3161. **************************
  3162.  'SI' (ERR1 > ERRTOL) ;
  3163.       'MESSAGE' 'Problem 48 T3';
  3164.       'ERREUR' 5 ;
  3165.  'FINSI' ;
  3166. **************************************************************
  3167. *********  Seventh line !!!!!!!!!!!!!!!!!!!!!!   *************
  3168. **************************************************************
  3169.  'SI' (('ABS'(DRYO2)) < DELTA) ;
  3170.     ERR1 = ('ABS' (DRYO2N '-' DRYO2)) ;
  3171.  'SINON' ;
  3172.     ERR1 = ('ABS' (DRYO2N '-' DRYO2)) '/' ('ABS' (DRYO2)) ;
  3173.  'FINSI' ;
  3174. *****************
  3175.  'SI' (ERR1 > ERRTOL) ;
  3176.         'MESSAGE' 'Problem 49 T3';
  3177.         'ERREUR' 5 ;
  3178.  'FINSI' ;
  3179. **************************************************************
  3180.  'SI' (('ABS'(DGXYO2)) < DELTA) ;
  3181.     ERR1 = ('ABS' (DGXYO2N '-' DGXYO2)) ;
  3182.  'SINON' ;
  3183.     ERR1 = ('ABS' (DGXYO2N '-' DGXYO2)) '/' ('ABS' (DGXYO2)) ;
  3184.  'FINSI' ;
  3185. *****************
  3186.  'SI' (ERR1 > ERRTOL) ;
  3187.         'MESSAGE' 'Problem 50 T3';
  3188.         'ERREUR' 5 ;
  3189.  'FINSI' ;
  3190. *************************************************************
  3191.  'SI' (('ABS'(DGYYO2)) < DELTA) ;
  3192.     ERR1 = ('ABS' (DGYYO2N '-' DGYYO2)) ; 
  3193.  'SINON' ;
  3194.     ERR1 = ('ABS' (DGYYO2N '-' DGYYO2)) '/' ('ABS' (DGYYO2)) ;
  3195.  'FINSI' ;
  3196. *****************    
  3197.  'SI' (ERR1 > ERRTOL) ;
  3198.         'MESSAGE' 'Problem 51 T3';
  3199.         'ERREUR' 5 ;
  3200.  'FINSI' ;
  3201. **************************************************************
  3202.  'SI' (('ABS'(DGZYO2)) < DELTA) ;
  3203.     ERR1 = ('ABS' (DGZYO2N '-' DGZYO2)) ; 
  3204.  'SINON' ;
  3205.     ERR1 = ('ABS' (DGZYO2N '-' DGZYO2)) '/' ('ABS' (DGZYO2)) ;
  3206.  'FINSI' ;
  3207. *****************    
  3208.  'SI' (ERR1 > ERRTOL) ;
  3209.         'MESSAGE' 'Problem 52 T3';
  3210.         'ERREUR' 5 ;
  3211.  'FINSI' ; 
  3212. **************************************************************
  3213.  'SI' (('ABS'(DRETYO2)) < DELTA) ;
  3214.     ERR1 = ('ABS' (DRETYO2N '-' DRETYO2)) ;
  3215.  'SINON' ;
  3216.     ERR1 = ('ABS' (DRETYO2N '-' DRETYO2)) '/' ('ABS' (DRETYO2)) ;
  3217.  'FINSI' ;
  3218. *****************
  3219.  'SI' (ERR1 > ERRTOL) ;
  3220.         'MESSAGE' 'Problem  53 T3';
  3221.         'ERREUR' 5 ;
  3222.  'FINSI' ;
  3223. **************************************************************
  3224.  'SI' (('ABS'(DYH2YO2)) < DELTA) ;
  3225.     ERR1 = ('ABS' (DYH2YO2N '-' DYH2YO2)) ;
  3226.  'SINON' ;  
  3227.     ERR1 = ('ABS' (DYH2YO2N '-' DYH2YO2)) '/' ('ABS' (DYH2YO2)) ;
  3228.  'FINSI' ;
  3229. **************************
  3230.  'SI' (ERR1 > ERRTOL) ;
  3231.       'MESSAGE' 'Problem 54 T3';
  3232.       'ERREUR' 5 ;
  3233.  'FINSI' ;
  3234. **************************************************************
  3235.  'SI' (('ABS'(DYO2YO2)) < DELTA) ;
  3236.     ERR1 = ('ABS' (DYO2YO2N '-' DYO2YO2)) ;
  3237.  'SINON' ;  
  3238.     ERR1 = ('ABS' (DYO2YO2N '-' DYO2YO2)) '/' ('ABS' (DYO2YO2)) ;
  3239.  'FINSI' ;
  3240. **************************
  3241.  'SI' (ERR1 > ERRTOL) ;
  3242.       'MESSAGE' 'Problem 55 T3';
  3243.       'ERREUR' 5 ;
  3244.  'FINSI' ;
  3245. **************************************************************
  3246.  'SI' (('ABS'(DYHOYO2)) < DELTA) ;
  3247.     ERR1 = ('ABS' (DYHOYO2N '-' DYHOYO2)) ;
  3248.  'SINON' ;  
  3249.     ERR1 = ('ABS' (DYHOYO2N '-' DYHOYO2)) '/' ('ABS' (DYHOYO2)) ;
  3250.  'FINSI' ;
  3251. **************************
  3252.  'SI' (ERR1 > ERRTOL) ;
  3253.       'MESSAGE' 'Problem 56 T3';
  3254.       'ERREUR' 5 ;
  3255.  'FINSI' ;
  3256. **************************************************************
  3257. *********  Eigth line !!!!!!!!!!!!!!!!!!!!!!   ***************
  3258. **************************************************************
  3259.  'SI' (('ABS'(DRYHO)) < DELTA) ;
  3260.     ERR1 = ('ABS' (DRYHON '-' DRYHO)) ;
  3261.  'SINON' ;
  3262.     ERR1 = ('ABS' (DRYHON '-' DRYHO)) '/' ('ABS' (DRYHO)) ;
  3263.  'FINSI' ;
  3264. *****************
  3265.  'SI' (ERR1 > ERRTOL) ;
  3266.         'MESSAGE' 'Problem 57 T3';
  3267.         'ERREUR' 5 ;
  3268.  'FINSI' ;
  3269. **************************************************************
  3270.  'SI' (('ABS'(DGXYHO)) < DELTA) ;
  3271.     ERR1 = ('ABS' (DGXYHON '-' DGXYHO)) ;
  3272.  'SINON' ;
  3273.     ERR1 = ('ABS' (DGXYHON '-' DGXYHO)) '/' ('ABS' (DGXYHO)) ;
  3274.  'FINSI' ;
  3275. *****************
  3276.  'SI' (ERR1 > ERRTOL) ;
  3277.         'MESSAGE' 'Problem 58 T3';
  3278.         'ERREUR' 5 ;
  3279.  'FINSI' ;
  3280. *************************************************************
  3281.  'SI' (('ABS'(DGYYHO)) < DELTA) ;
  3282.     ERR1 = ('ABS' (DGYYHON '-' DGYYHO)) ; 
  3283.  'SINON' ;
  3284.     ERR1 = ('ABS' (DGYYHON '-' DGYYHO)) '/' ('ABS' (DGYYHO)) ;
  3285.  'FINSI' ;
  3286. *****************    
  3287.  'SI' (ERR1 > ERRTOL) ;
  3288.         'MESSAGE' 'Problem 59 T3';
  3289.         'ERREUR' 5 ;
  3290.  'FINSI' ;
  3291. **************************************************************
  3292.  'SI' (('ABS'(DGZYHO)) < DELTA) ;
  3293.     ERR1 = ('ABS' (DGZYHON '-' DGZYHO)) ; 
  3294.  'SINON' ;
  3295.     ERR1 = ('ABS' (DGZYHON '-' DGZYHO)) '/' ('ABS' (DGZYHO)) ;
  3296.  'FINSI' ;
  3297. *****************    
  3298.  'SI' (ERR1 > ERRTOL) ;
  3299.         'MESSAGE' 'Problem 60 T3';
  3300.         'ERREUR' 5 ;
  3301.  'FINSI' ; 
  3302. **************************************************************
  3303.  'SI' (('ABS'(DRETYHO)) < DELTA) ;
  3304.     ERR1 = ('ABS' (DRETYHON '-' DRETYHO)) ;
  3305.  'SINON' ;
  3306.     ERR1 = ('ABS' (DRETYHON '-' DRETYHO)) '/' ('ABS' (DRETYHO)) ;
  3307.  'FINSI' ;
  3308. *****************
  3309.  'SI' (ERR1 > ERRTOL) ;
  3310.         'MESSAGE' 'Problem 61 T3';
  3311.         'ERREUR' 5 ;
  3312.  'FINSI' ;
  3313. **************************************************************
  3314.  'SI' (('ABS'(DYH2YHO)) < DELTA) ;
  3315.     ERR1 = ('ABS' (DYH2YHON '-' DYH2YHO)) ;
  3316.  'SINON' ;  
  3317.     ERR1 = ('ABS' (DYH2YHON '-' DYH2YHO)) '/' ('ABS' (DYH2YHO)) ;
  3318.  'FINSI' ;
  3319. **************************
  3320.  'SI' (ERR1 > ERRTOL) ;
  3321.       'MESSAGE' 'Problem 62 T3';
  3322.       'ERREUR' 5 ;
  3323.  'FINSI' ;
  3324. **************************************************************
  3325.  'SI' (('ABS'(DYO2YHO)) < DELTA) ;
  3326.     ERR1 = ('ABS' (DYO2YHON '-' DYO2YHO)) ;
  3327.  'SINON' ;  
  3328.     ERR1 = ('ABS' (DYO2YHON '-' DYO2YHO)) '/' ('ABS' (DYO2YHO)) ;
  3329.  'FINSI' ;
  3330. **************************
  3331.  'SI' (ERR1 > ERRTOL) ;
  3332.       'MESSAGE' 'Problem 63 T3';
  3333.       'ERREUR' 5 ;
  3334.  'FINSI' ;
  3335. **************************************************************
  3336.  'SI' (('ABS'(DYHOYHO)) < DELTA) ;
  3337.     ERR1 = ('ABS' (DYHOYHON '-' DYHOYHO)) ;
  3338.  'SINON' ;  
  3339.     ERR1 = ('ABS' (DYHOYHON '-' DYHOYHO)) '/' ('ABS' (DYHOYHO)) ;
  3340.  'FINSI' ;
  3341. **************************
  3342.  'SI' (ERR1 > ERRTOL) ;
  3343.       'MESSAGE' 'Problem 64 T3';
  3344.       'ERREUR' 5 ;
  3345.  'FINSI' ;
  3346. *************************************************************  
  3347.  
  3348.  'FIN' ;
  3349.  
  3350.  
  3351.  
  3352.  
  3353.  
  3354.  
  3355.  
  3356.  
  3357.  
  3358.  
  3359.  
  3360.  
  3361.  

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