Télécharger la plaquette Cast3M


fr en

Télécharger kbmmsp_impl2D.dgibi

Retour à la liste

Numérotation des lignes : 

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

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