 problem: Rosenbrock Function                             


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to        94 (  0.1 percent)

 Integer Workspace declared as      1000 is used up to        54 (  5.4 percent)



 N =   2

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     2
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.348E+01       0.366E+01                  0         2   0.576E+01
         1      0.345E+01    *  0.362E+01      0.01000
         1      0.345E+01       0.356E+01                  0         2   0.560E+01
         2      0.116E+02    *  0.753E+00      0.55029
         2      0.116E+02       0.115E+01                  0         2   0.100E+01
         3      0.664E+01    *  0.523E+00      0.87723
         3      0.664E+01       0.818E+00                  0         2   0.352E+01
         4      0.102E+00    *  0.102E-01      1.00000
         4      0.102E+00       0.102E-01                  0         2   0.400E+01
         5      0.000E+00    *  0.000E+00      1.00000
         5      0.000E+00       0.000E+00                  0         2   0.400E+01
         6      0.000E+00    *  0.000E+00      1.00000



 Solution of nonlinear system of equations obtained within   6 iteration steps

 Achieved relative accuracy 0.000E+00

   Subcondition ( 1,   2)  0.400E+01
   Sensitivity ( 1,   2)  0.141E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       6  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       6  ***
   ***  Function eval.    :       7  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.0000E+00

 problem: Powell Singular Function                        


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       134 (  0.2 percent)

 Integer Workspace declared as      1000 is used up to        56 (  5.6 percent)



 N =   4

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     4
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.733E+01       0.677E+00                  0         4   0.203E+01
         1      0.726E+01    *  0.670E+00      0.01000
         1      0.726E+01       0.671E+00                  0         4   0.203E+01
         2      0.157E+01    *  0.616E-01      1.00000
         2      0.157E+01       0.160E+00                  0         4   0.203E+01
         3      0.393E+00    *  0.400E-01      1.00000
         3      0.393E+00       0.153E+00                  0         4   0.243E+01
         4      0.981E-01    *  0.381E-01      1.00000
         4      0.981E-01       0.763E-01                  0         4   0.243E+01
         5      0.245E-01    *  0.191E-01      1.00000
         5      0.245E-01       0.381E-01                  0         4   0.243E+01
         6      0.613E-02    *  0.953E-02      1.00000
         6      0.613E-02       0.191E-01                  0         4   0.243E+01
         7      0.153E-02    *  0.477E-02      1.00000
         7      0.153E-02       0.953E-02                  0         4   0.243E+01
         8      0.383E-03    *  0.238E-02      1.00000
         8      0.383E-03       0.477E-02                  0         4   0.243E+01
         9      0.958E-04    *  0.119E-02      1.00000
         9      0.958E-04       0.238E-02                  0         4   0.243E+01
        10      0.240E-04    *  0.596E-03      1.00000
        10      0.240E-04       0.119E-02                  0         4   0.243E+01
        11      0.599E-05    *  0.298E-03      1.00000
        11      0.599E-05       0.596E-03                  0         4   0.243E+01
        12      0.150E-05    *  0.149E-03      1.00000
        12      0.150E-05       0.298E-03                  0         4   0.243E+01
        13      0.374E-06    *  0.745E-04      1.00000
        13      0.374E-06       0.149E-03                  0         4   0.243E+01
        14      0.936E-07    *  0.372E-04      1.00000
        14      0.936E-07       0.745E-04                  0         4   0.243E+01
        15      0.234E-07    *  0.186E-04      1.00000
        15      0.234E-07       0.372E-04                  0         4   0.243E+01
        16      0.585E-08    *  0.931E-05      1.00000
        16      0.585E-08       0.186E-04                  0         4   0.243E+01
        17      0.146E-08    *  0.465E-05      1.00000
        17      0.146E-08       0.931E-05                  0         4   0.243E+01
        18      0.366E-09    *  0.233E-05      1.00000
        18      0.366E-09       0.465E-05                  0         4   0.243E+01
        19      0.914E-10    *  0.116E-05      1.00000
        19      0.914E-10       0.233E-05                  0         4   0.243E+01
        20      0.229E-10    *  0.582E-06      1.00000
        20      0.229E-10       0.116E-05                  0         4   0.243E+01
        21      0.571E-11    *  0.291E-06      1.00000
        21      0.571E-11       0.582E-06                  0         4   0.243E+01
        22      0.143E-11    *  0.145E-06      1.00000
        22      0.143E-11       0.291E-06                  0         4   0.243E+01
        23      0.357E-12    *  0.727E-07      1.00000
        23      0.357E-12       0.145E-06                  0         4   0.243E+01
        24      0.893E-13    *  0.364E-07      1.00000
        24      0.893E-13       0.727E-07                  0         4   0.243E+01
        25      0.223E-13    *  0.182E-07      1.00000
        25      0.223E-13       0.364E-07                  0         4   0.243E+01
        26      0.558E-14    *  0.909E-08      1.00000
        26      0.558E-14       0.182E-07                  0         4   0.243E+01
        27      0.139E-14    *  0.455E-08      1.00000
        27      0.139E-14       0.909E-08                  0         4   0.243E+01
        28      0.349E-15    *  0.227E-08      1.00000
        28      0.349E-15       0.455E-08                  0         4   0.243E+01
        29      0.872E-16    *  0.114E-08      1.00000
        29      0.872E-16       0.227E-08                  0         4   0.243E+01
        30      0.218E-16    *  0.568E-09      1.00000
        30      0.218E-16       0.114E-08                  0         4   0.243E+01
        31      0.545E-17    *  0.284E-09      1.00000
        31      0.545E-17       0.568E-09                  0         4   0.243E+01
        32      0.136E-17    *  0.142E-09      1.00000
        32      0.136E-17       0.284E-09                  0         4   0.243E+01
        33      0.340E-18    *  0.710E-10      1.00000



 Solution of nonlinear system of equations obtained within  33 iteration steps

 Achieved relative accuracy 0.710E-10

 Warning: No quadratic or superlinear convergence established yet
          your solution may perhaps may be less accurate 
          as indicated by the standard error estimate

   Subcondition ( 1,   4)  0.243E+01
   Sensitivity ( 1,   4)  0.141E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      33  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      33  ***
   ***  Function eval.    :      34  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1915E-18

 problem: Powell Badly Scaled Function                    


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to        94 (  0.1 percent)

 Integer Workspace declared as      1000 is used up to        54 (  5.4 percent)



 N =   2

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     2
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.753E+00       0.707E+00                  0         2   0.544E+01
         1      0.746E+00    *  0.700E+00      0.01000
         1      0.746E+00       0.703E+00                  0         2   0.546E+01
         2      0.682E+00    *  0.260E+00      0.99501
         2      0.682E+00       0.469E+00                  0         2   0.985E+01
         3      0.719E-01    *  0.261E+00      0.58428
         3      0.719E-01       0.307E+00                  0         2   0.116E+02
         4      0.639E-01    *  0.147E+00      0.73779
         4      0.639E-01       0.238E+00                  0         2   0.188E+02
         5      0.397E-01    *  0.121E+00      0.67876
         5      0.397E-01       0.192E+00                  0         2   0.299E+02
         6      0.277E-01    *  0.947E-01      0.70405
         6      0.277E-01       0.160E+00                  0         2   0.506E+02
         7      0.203E-01    *  0.794E-01      0.69948
         7      0.203E-01       0.137E+00                  0         2   0.871E+02
         8      0.154E-01    *  0.665E-01      0.71181
         8      0.154E-01       0.117E+00                  0         2   0.153E+03
         9      0.119E-01    *  0.555E-01      0.72874
         9      0.119E-01       0.991E-01                  0         2   0.273E+03
        10      0.907E-02    *  0.442E-01      0.76856
        10      0.907E-02       0.798E-01                  0         2   0.494E+03
        11      0.653E-02    *  0.302E-01      0.86073
        11      0.653E-02       0.552E-01                  0         2   0.901E+03
        12      0.343E-02    *  0.139E-01      1.00000
        12      0.343E-02       0.237E-01                  0         2   0.154E+04
        13      0.628E-03    *  0.309E-02      1.00000
        13      0.628E-03       0.391E-02                  0         2   0.194E+04
        14      0.175E-04    *  0.958E-04      1.00000
        14      0.175E-04       0.988E-04                  0         2   0.200E+04
        15      0.113E-07    *  0.632E-07      1.00000
        15      0.113E-07       0.631E-07                  0         2   0.200E+04
        16      0.457E-14    *  0.502E-16      1.00000



 Solution of nonlinear system of equations obtained within  16 iteration steps

 Achieved relative accuracy 0.502E-16

   Subcondition ( 1,   2)  0.200E+04
   Sensitivity ( 1,   2)  0.141E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      16  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      16  ***
   ***  Function eval.    :      17  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1491E-16

 problem: Wood Function                                   


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       134 (  0.2 percent)

 Integer Workspace declared as      1000 is used up to        56 (  5.6 percent)



 N =   4

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     4
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.428E+04       0.497E+01                  0         4   0.846E+02
         1      0.423E+04    *  0.492E+01      0.01000
         1      0.423E+04       0.490E+01                  0         4   0.848E+02
         2      0.418E+03    *  0.710E+00      1.00000
         2      0.418E+03       0.722E+00                  0         4   0.844E+02
         3      0.256E+03    *  0.161E+00      0.81339
         3      0.256E+03       0.240E+00                  0         4   0.307E+02
         4      0.976E+02    *  0.544E-01      1.00000
         4      0.976E+02       0.204E+00                  0         4   0.376E+02
         5      0.334E+02    *  0.493E-01      1.00000
         5      0.334E+02       0.155E+00                  0         4   0.537E+02
         6      0.966E+01    *  0.322E-01      1.00000
         6      0.966E+01       0.788E-01                  0         4   0.112E+03
         7      0.161E+01    *  0.110E-01      1.00000
         7      0.161E+01       0.192E-01                  0         4   0.605E+03
         8      0.703E-01    *  0.600E-02      1.00000
         8      0.703E-01       0.579E-01                  0         4   0.482E+04
         9      0.683E-01    *  0.560E-01      0.03209
         9      0.683E-01       0.548E-01                  0         4   0.477E+04
        10      0.237E+00    *  0.435E-01      1.00000
        10      0.237E+00       0.104E+00                  0         4   0.235E+05
        11      0.221E+00    *  0.811E-01      0.15641
        11      0.221E+00       0.517E-02                  0         4   0.170E+05
        11      0.263E-02    *  0.200E-01      1.00000
        12      0.192E+00    *  0.417E-02      0.12897
        12      0.192E+00       0.130E-02                  0         4   0.113E+05
        13      0.521E-01    *  0.119E-02      0.72945
        13      0.521E-01       0.671E-03                  0         4   0.424E+04
        14      0.542E-04    *  0.152E-03      1.00000
        14      0.542E-04       0.124E-03                  0         4   0.344E+04
        15      0.121E-05    *  0.200E-07      1.00000
        15      0.121E-05       0.205E-07                  0         4   0.344E+04
        16      0.306E-13    *  0.867E-13      1.00000



 Solution of nonlinear system of equations obtained within  16 iteration steps

 Achieved relative accuracy 0.867E-13

   Subcondition ( 1,   4)  0.344E+04
   Sensitivity ( 1,   4)  0.141E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      16  ***
   ***  Corrector steps   :       1  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      16  ***
   ***  Function eval.    :      18  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.2060E-13

 problem: Helical Valley Function                         


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       113 (  0.2 percent)

 Integer Workspace declared as      1000 is used up to        55 (  5.5 percent)



 N =   3

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     3
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.289E+02       0.181E+01                  0         3   0.139E+01
         1      0.286E+02    *  0.180E+01      0.01000
         1      0.286E+02       0.180E+01                  0         3   0.139E+01
         2      0.215E+02    *  0.136E+01      0.32137
         2      0.215E+02       0.191E+01                  0         3   0.148E+01
         3      0.155E+02    *  0.139E+01      0.32066
         3      0.155E+02       0.160E+01                  0         3   0.208E+01
         4      0.906E+01    *  0.843E+00      0.60661
         4      0.906E+01       0.706E+00                  0         3   0.202E+01
         5      0.518E+01    *  0.553E+00      1.00000
         5      0.518E+01       0.442E+00                  0         3   0.225E+01
         6      0.174E+01    *  0.174E+00      1.00000
         6      0.174E+01       0.125E+00                  0         3   0.148E+01
         7      0.184E+00    *  0.121E-01      1.00000
         7      0.184E+00       0.140E-01                  0         3   0.141E+01
         8      0.289E-02    *  0.192E-03      1.00000
         8      0.289E-02       0.191E-03                  0         3   0.139E+01
         9      0.456E-06    *  0.351E-07      1.00000
         9      0.456E-06       0.353E-07                  0         3   0.139E+01
        10      0.175E-13    *  0.114E-14      1.00000



 Solution of nonlinear system of equations obtained within  10 iteration steps

 Achieved relative accuracy 0.114E-14

   Subcondition ( 1,   3)  0.139E+01
   Sensitivity ( 1,   3)  0.118E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      10  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      10  ***
   ***  Function eval.    :      11  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.5323E-21

 problem: Watson Function                                 


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       302 (  0.4 percent)

 Integer Workspace declared as      1000 is used up to        62 (  6.2 percent)



 N =  10

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    10
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.300E+02       0.331E+02                  0        10   0.442E+11
         1      0.297E+02    *  0.327E+02      0.01000
         1      0.297E+02       0.325E+02                  0        10   0.441E+11
         2      0.548E+01    *  0.937E+01      1.00000
         2      0.548E+01       0.721E+00                  0        10   0.333E+11
         3      0.633E+00    *  0.987E-01      1.00000
         3      0.633E+00       0.341E+00                  0        10   0.303E+11
         4      0.567E+00    *  0.227E+00      0.33070
         4      0.567E+00       0.362E+00                  0        10   0.301E+11
         5      0.102E+01    *  0.179E+00      1.00000
         5      0.102E+01       0.789E-01                  0        10   0.259E+11
         6      0.228E+00    *  0.763E-01      1.00000
         6      0.228E+00       0.335E+00                  0        10   0.233E+11
         7      0.206E+00    *  0.277E+00      0.17791
         7      0.206E+00       0.336E+00                  0        10   0.221E+11
         8      0.530E+00    *  0.255E+00      1.00000
         8      0.530E+00       0.396E-01                  0        10   0.177E+11
         8      0.602E-01    *  0.538E-01      1.00000
         9      0.343E+00    *  0.310E-01      0.36867
         9      0.343E+00       0.935E-01                  0        10   0.146E+11
        10      0.293E+00    *  0.811E-01      0.15073
        10      0.293E+00       0.106E+00                  0        10   0.148E+11
        11      0.167E+00    *  0.696E-01      0.49780
        11      0.167E+00       0.133E+00                  0        10   0.150E+11
        12      0.106E+00    *  0.749E-01      0.51974
        12      0.106E+00       0.127E+00                  0        10   0.164E+11
        13      0.784E-01    *  0.142E-01      1.00000
        13      0.784E-01       0.556E-01                  0        10   0.178E+11
        14      0.241E-01    *  0.184E-01      0.90666
        14      0.241E-01       0.340E-01                  0        10   0.164E+11
        15      0.693E-02    *  0.395E-02      1.00000
        15      0.693E-02       0.622E-02                  0        10   0.161E+11
        16      0.265E-03    *  0.331E-03      1.00000
        16      0.265E-03       0.364E-03                  0        10   0.160E+11
        17      0.866E-06    *  0.811E-06      1.00000
        17      0.866E-06       0.815E-06                  0        10   0.161E+11
        18      0.457E-11    *  0.605E-11      1.00000



 Solution of nonlinear system of equations obtained within  18 iteration steps

 Achieved relative accuracy 0.605E-11

   Subcondition ( 1,  10)  0.161E+11
   Sensitivity ( 1,  10)  0.316E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      18  ***
   ***  Corrector steps   :       1  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      18  ***
   ***  Function eval.    :      20  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1603E-13

 problem: Chebyquad Function, n=7                         


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       209 (  0.3 percent)

 Integer Workspace declared as      1000 is used up to        59 (  5.9 percent)



 N =   7

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     7
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.695E-01       0.120E+00                  0         7   0.100E+02
         1      0.688E-01    *  0.119E+00      0.01000
         1      0.688E-01       0.116E+00                  0         7   0.995E+01
         2      0.312E-01    *  0.320E-01      0.42935
         2      0.312E-01       0.156E+00                  0         7   0.223E+02
         3      0.285E-01    *  0.139E+00      0.05767
         3      0.285E-01       0.507E-01                  0         7   0.827E+01
         4      0.193E-01    *  0.306E-01      0.27534
         4      0.193E-01       0.180E-01                  0         7   0.394E+01
         5      0.694E-02    *  0.403E-02      1.00000
         5      0.694E-02       0.294E-02                  0         7   0.243E+01
         6      0.135E-03    *  0.647E-04      1.00000
         6      0.135E-03       0.685E-04                  0         7   0.249E+01
         7      0.121E-06    *  0.575E-07      1.00000
         7      0.121E-06       0.576E-07                  0         7   0.249E+01
         8      0.785E-13    *  0.379E-13      1.00000



 Solution of nonlinear system of equations obtained within   8 iteration steps

 Achieved relative accuracy 0.379E-13

   Subcondition ( 1,   7)  0.249E+01
   Sensitivity ( 1,   7)  0.236E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       8  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       8  ***
   ***  Function eval.    :       9  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1358E-15

 problem: Chebyquad Function, n=9                         


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       269 (  0.4 percent)

 Integer Workspace declared as      1000 is used up to        61 (  6.1 percent)



 N =   9

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     9
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.566E-01       0.237E+00                  0         9   0.401E+02
         1      0.561E-01    *  0.235E+00      0.01000
         1      0.561E-01       0.224E+00                  0         9   0.394E+02
         2      0.655E-01    *  0.166E+00      0.19826
         2      0.655E-01       0.399E+00                  0         9   0.778E+02
         2      0.672E-01    *  0.424E+00      0.03445
         3      0.652E-01    *  0.398E+00      0.00615
         3      0.652E-01       0.229E+02                  0         9   0.442E+04
         3      0.653E-01    *  0.229E+02      0.00010
         3 Not accepted damping factor        0.00000     0         9
         3      0.652E-01       0.223E+02                  0         8   0.344E+04
         3 Not accepted damping factor        0.00000     0         8
         3      0.652E-01       0.365E-01                  0         7   0.168E+02
         4      0.451E-01    *  0.273E-01      0.21755
         4      0.451E-01       0.131E+02                  0         9   0.486E+04
         4      0.451E-01    *  0.131E+02      0.00010
         4 Not accepted damping factor        0.00000     0         9
         4      0.451E-01       0.162E+01                  0         8   0.324E+03
         4 Not accepted damping factor        0.00008     0         8
         4      0.451E-01       0.213E-01                  0         7   0.918E+01
         5      0.336E-01    *  0.111E-01      0.98445
         5      0.336E-01       0.389E+00                  0         9   0.647E+03
         6      0.335E-01    *  0.386E+00      0.00155
         6      0.335E-01       0.523E-01                  0         9   0.867E+02
         7      0.330E-01    *  0.512E-01      0.01329
         7      0.330E-01       0.266E-01                  0         9   0.437E+02
         8      0.313E-01    *  0.246E-01      0.05310
         8      0.313E-01       0.138E-01                  0         9   0.222E+02
         9      0.244E-01    *  0.983E-02      0.21583
         9      0.244E-01       0.680E-02                  0         9   0.114E+02
        10      0.996E-03    *  0.203E-02      1.00000
        10      0.996E-03       0.124E-02                  0         9   0.620E+01
        11      0.520E-04    *  0.572E-04      1.00000
        11      0.520E-04       0.631E-04                  0         9   0.677E+01
        12      0.145E-06    *  0.167E-06      1.00000
        12      0.145E-06       0.167E-06                  0         9   0.680E+01
        13      0.120E-11    *  0.124E-11      1.00000



 Solution of nonlinear system of equations obtained within  13 iteration steps

 Achieved relative accuracy 0.124E-11

   Subcondition ( 1,   9)  0.680E+01
   Sensitivity ( 1,   9)  0.273E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      13  ***
   ***  Corrector steps   :       3  ***
   ***  Rejected rk-1 st. :       4  ***
   ***  Jacobian eval.    :      13  ***
   ***  Function eval.    :      17  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.8562E-16

 problem: Brown Almost Linear Function                    


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       302 (  0.4 percent)

 Integer Workspace declared as      1000 is used up to        62 (  6.2 percent)



 N =  10

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    10
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.523E+01       0.167E+04                  0        10   0.233E+02
         0      0.138E+08    *  0.737E+11      0.01000
         0      0.523E+01    *  0.167E+04      0.00010
         0 Not accepted damping factor        0.00002     0        10
         0      0.523E+01       0.140E+02                  0         9   0.377E+01
         1      0.521E+01    *  0.140E+02      0.00010
         1      0.521E+01       0.163E+04                  0        10   0.233E+02
         1      0.521E+01    *  0.163E+04      0.00010
         1 Not accepted damping factor        0.00002     0        10
         1      0.521E+01       0.137E+02                  0         9   0.377E+01
         2      0.462E+01    *  0.136E+02      0.00419
         2      0.462E+01       0.563E+03                  0        10   0.217E+02
         2      0.462E+01    *  0.563E+03      0.00010
         2 Not accepted damping factor        0.00004     0        10
         2      0.462E+01       0.529E+01                  0         9   0.377E+01
         3      0.364E+01    *  0.517E+01      0.01771
         3      0.364E+01       0.120E+03                  0        10   0.201E+02
         4      0.364E+01    *  0.120E+03      0.00010
         4      0.364E+01       0.243E+03                  0        10   0.407E+02
         5      0.364E+01    *  0.243E+03      0.00010
         5      0.364E+01       0.298E+03                  0        10   0.504E+02
         5      0.364E+01    *  0.298E+03      0.00010
         5 Not accepted damping factor        0.00004     0        10
         5      0.364E+01       0.146E+01                  0         9   0.377E+01
         5      0.570E+03    *  0.221E+04      1.00000
         6      0.364E+01    *  0.146E+01      0.00033
         6      0.364E+01       0.297E+03                  0        10   0.506E+02
         6      0.363E+01    *  0.297E+03      0.00010
         6 Not accepted damping factor        0.00004     0        10
         6      0.364E+01       0.146E+01                  0         9   0.377E+01
         7      0.272E+01    *  0.135E+01      0.06059
         7      0.272E+01       0.184E+03                  0        10   0.115E+03
         7      0.272E+01    *  0.184E+03      0.00010
         7 Not accepted damping factor        0.00003     0        10
         7      0.272E+01       0.560E+00                  0         9   0.377E+01
         8      0.114E+01    *  0.344E+00      0.26952
         8      0.114E+01       0.198E+02                  0        10   0.175E+03
         9      0.114E+01    *  0.198E+02      0.00013
         9      0.114E+01       0.517E+02                  0        10   0.457E+03
         9      0.114E+01    *  0.517E+02      0.00010
         9 Not accepted damping factor        0.00003     0        10
         9      0.114E+01       0.129E+00                  0         9   0.377E+01
        10      0.236E+00    *  0.334E-01      1.00000
        10      0.236E+00       0.181E+00                  0        10   0.833E+02
        10      0.204E+00    *  0.182E+00      0.13225
        11      0.221E+00    *  0.175E+00      0.06226
        11      0.221E+00       0.771E+00                  0        10   0.432E+03
        12      0.220E+00    *  0.771E+00      0.00235
        12      0.220E+00       0.825E+01                  0        10   0.465E+04
        12      0.220E+00    *  0.825E+01      0.00010
        12 Not accepted damping factor        0.00002     0        10
        12      0.220E+00       0.193E-01                  0         9   0.377E+01
        13      0.435E-02    *  0.418E-03      1.00000
        13      0.435E-02       0.260E+01                  0        10   0.381E+04
        14      0.435E-02    *  0.260E+01      0.00010
        14      0.435E-02       0.108E+01                  0        10   0.158E+04
        15      0.435E-02    *  0.108E+01      0.00041
        15      0.435E-02       0.539E+00                  0        10   0.791E+03
        16      0.434E-02    *  0.538E+00      0.00165
        16      0.434E-02       0.269E+00                  0        10   0.396E+03
        17      0.431E-02    *  0.266E+00      0.00663
        17      0.431E-02       0.133E+00                  0        10   0.198E+03
        18      0.418E-02    *  0.128E+00      0.02680
        18      0.418E-02       0.639E-01                  0        10   0.988E+02
        19      0.367E-02    *  0.532E-01      0.11170
        19      0.367E-02       0.266E-01                  0        10   0.495E+02
        20      0.152E-02    *  0.488E-02      0.53756
        20      0.152E-02       0.247E-02                  0        10   0.249E+02
        21      0.947E-05    *  0.104E-03      1.00000
        21      0.947E-05       0.983E-04                  0        10   0.232E+02
        22      0.152E-07    *  0.158E-06      1.00000
        22      0.152E-07       0.159E-06                  0        10   0.233E+02
        23      0.396E-13    *  0.413E-12      1.00000



 Solution of nonlinear system of equations obtained within  23 iteration steps

 Achieved relative accuracy 0.413E-12

   Subcondition ( 1,  10)  0.233E+02
   Sensitivity ( 1,  10)  0.200E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      23  ***
   ***  Corrector steps   :      11  ***
   ***  Rejected rk-1 st. :       8  ***
   ***  Jacobian eval.    :      23  ***
   ***  Function eval.    :      35  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1685E-14

 problem: Discrete Boundary Function                      


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       302 (  0.4 percent)

 Integer Workspace declared as      1000 is used up to        62 (  6.2 percent)



 N =  10

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    10
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.888E-02       0.651E-01                  0        10   0.104E+02
         1      0.879E-02    *  0.644E-01      0.01000
         1      0.879E-02       0.644E-01                  0        10   0.104E+02
         2      0.758E-04    *  0.711E-03      1.00000
         2      0.758E-04       0.696E-03                  0        10   0.102E+02
         3      0.945E-08    *  0.883E-07      1.00000
         3      0.945E-08       0.884E-07                  0        10   0.102E+02
         4      0.156E-15    *  0.144E-14      1.00000



 Solution of nonlinear system of equations obtained within   4 iteration steps

 Achieved relative accuracy 0.144E-14

   Subcondition ( 1,  10)  0.102E+02
   Sensitivity ( 1,  10)  0.122E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       4  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       4  ***
   ***  Function eval.    :       5  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.2481E-16

 problem: Discrete Integral Function                      


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       302 (  0.4 percent)

 Integer Workspace declared as      1000 is used up to        62 (  6.2 percent)



 N =  10

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    10
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.796E-01       0.651E-01                  0        10   0.102E+01
         1      0.788E-01    *  0.644E-01      0.01000
         1      0.788E-01       0.644E-01                  0        10   0.102E+01
         2      0.886E-03    *  0.711E-03      1.00000
         2      0.886E-03       0.696E-03                  0        10   0.102E+01
         3      0.113E-06    *  0.883E-07      1.00000
         3      0.113E-06       0.884E-07                  0        10   0.102E+01
         4      0.186E-14    *  0.145E-14      1.00000



 Solution of nonlinear system of equations obtained within   4 iteration steps

 Achieved relative accuracy 0.145E-14

   Subcondition ( 1,  10)  0.102E+01
   Sensitivity ( 1,  10)  0.101E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       4  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       4  ***
   ***  Function eval.    :       5  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1335E-16

 problem: Trigonometric Function                          


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       302 (  0.4 percent)

 Integer Workspace declared as      1000 is used up to        62 (  6.2 percent)



 N =  10

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    10
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.266E-01       0.162E+00                  0        10   0.703E+01
         1      0.263E-01    *  0.161E+00      0.01000
         1      0.263E-01       0.141E+01                  0        10   0.597E+02
         2      0.263E-01    *  0.141E+01      0.00015
         2      0.263E-01       0.376E+01                  0        10   0.158E+03
         2      0.263E-01    *  0.376E+01      0.00010
         2 Not accepted damping factor        0.00004     0        10
         2      0.263E-01       0.348E-01                  0         9   0.218E+01
         3      0.192E-01    *  0.305E-01      0.45151
         3      0.192E-01       0.153E+01                  0        10   0.652E+02
         4      0.192E-01    *  0.153E+01      0.00020
         4      0.192E-01       0.108E+01                  0        10   0.464E+02
         5      0.191E-01    *  0.108E+01      0.00096
         5      0.191E-01       0.522E+00                  0        10   0.236E+02
         6      0.191E-01    *  0.519E+00      0.00384
         6      0.191E-01       0.239E+00                  0        10   0.119E+02
         7      0.187E-01    *  0.233E+00      0.01519
         7      0.187E-01       0.103E+00                  0        10   0.590E+01
         8      0.176E-01    *  0.944E-01      0.05413
         8      0.176E-01       0.671E-01                  0        10   0.271E+01
         9      0.154E-01    *  0.591E-01      0.10992
         9      0.154E-01       0.505E-01                  0        10   0.241E+01
        10      0.619E-02    *  0.138E-01      0.56850
        10      0.619E-02       0.950E-02                  0        10   0.170E+01
        11      0.610E-03    *  0.738E-03      1.00000
        11      0.610E-03       0.715E-03                  0        10   0.178E+01
        12      0.615E-05    *  0.135E-04      1.00000
        12      0.615E-05       0.142E-04                  0        10   0.183E+01
        13      0.348E-08    *  0.989E-08      1.00000
        13      0.348E-08       0.991E-08                  0        10   0.183E+01
        14      0.147E-14    *  0.451E-14      1.00000



 Solution of nonlinear system of equations obtained within  14 iteration steps

 Achieved relative accuracy 0.451E-14

   Subcondition ( 1,  10)  0.183E+01
   Sensitivity ( 1,  10)  0.136E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      14  ***
   ***  Corrector steps   :       1  ***
   ***  Rejected rk-1 st. :       1  ***
   ***  Jacobian eval.    :      14  ***
   ***  Function eval.    :      16  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1349E-14

 problem: Variably Dimensioned Function                   


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       302 (  0.4 percent)

 Integer Workspace declared as      1000 is used up to        62 (  6.2 percent)



 N =  10

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    10
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.708E+06       0.207E+00                  0        10   0.219E+07
         1      0.701E+06    *  0.205E+00      0.01000
         1      0.701E+06       0.206E+00                  0        10   0.218E+07
         2      0.208E+06    *  0.611E-01      1.00000
         2      0.208E+06       0.137E+00                  0        10   0.968E+06
         3      0.616E+05    *  0.407E-01      1.00000
         3      0.616E+05       0.917E-01                  0        10   0.430E+06
         4      0.183E+05    *  0.272E-01      1.00000
         4      0.183E+05       0.612E-01                  0        10   0.191E+06
         5      0.541E+04    *  0.181E-01      1.00000
         5      0.541E+04       0.409E-01                  0        10   0.848E+05
         6      0.161E+04    *  0.121E-01      1.00000
         6      0.161E+04       0.274E-01                  0        10   0.376E+05
         7      0.479E+03    *  0.815E-02      1.00000
         7      0.479E+03       0.184E-01                  0        10   0.166E+05
         8      0.143E+03    *  0.552E-02      1.00000
         8      0.143E+03       0.125E-01                  0        10   0.732E+04
         9      0.434E+02    *  0.380E-02      1.00000
         9      0.434E+02       0.874E-02                  0        10   0.318E+04
        10      0.134E+02    *  0.271E-02      1.00000
        10      0.134E+02       0.635E-02                  0        10   0.136E+04
        11      0.424E+01    *  0.200E-02      1.00000
        11      0.424E+01       0.471E-02                  0        10   0.577E+03
        12      0.120E+01    *  0.133E-02      0.99733
        12      0.120E+01       0.260E-02                  0        10   0.295E+03
        13      0.123E+00    *  0.266E-03      1.00000
        13      0.123E+00       0.317E-03                  0        10   0.247E+03
        14      0.189E-03    *  0.489E-06      1.00000
        14      0.189E-03       0.490E-06                  0        10   0.247E+03
        15      0.699E-12    *  0.181E-14      1.00000



 Solution of nonlinear system of equations obtained within  15 iteration steps

 Achieved relative accuracy 0.181E-14

   Subcondition ( 1,  10)  0.247E+03
   Sensitivity ( 1,  10)  0.316E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      15  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      15  ***
   ***  Function eval.    :      16  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.0000E+00

 problem: Broyden Tridiagonal Function                    


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       302 (  0.4 percent)

 Integer Workspace declared as      1000 is used up to        62 (  6.2 percent)



 N =  10

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    10
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.145E+01       0.303E+00                  0        10   0.121E+01
         1      0.143E+01    *  0.300E+00      0.01000
         1      0.143E+01       0.301E+00                  0        10   0.122E+01
         2      0.205E+00    *  0.431E-01      1.00000
         2      0.205E+00       0.608E-01                  0        10   0.133E+01
         3      0.915E-02    *  0.274E-02      1.00000
         3      0.915E-02       0.303E-02                  0        10   0.137E+01
         4      0.276E-04    *  0.927E-05      1.00000
         4      0.276E-04       0.933E-05                  0        10   0.138E+01
         5      0.304E-09    *  0.102E-09      1.00000
         5      0.304E-09       0.102E-09                  0        10   0.138E+01
         6      0.527E-15    *  0.889E-16      1.00000



 Solution of nonlinear system of equations obtained within   6 iteration steps

 Achieved relative accuracy 0.889E-16

   Subcondition ( 1,  10)  0.138E+01
   Sensitivity ( 1,  10)  0.108E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       6  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       6  ***
   ***  Function eval.    :       7  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.3988E-15

 problem: Broyden Banded Function                         


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       302 (  0.4 percent)

 Integer Workspace declared as      1000 is used up to        62 (  6.2 percent)



 N =  10

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    10
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.600E+01       0.281E+00                  0        10   0.102E+01
         1      0.594E+01    *  0.278E+00      0.01000
         1      0.594E+01       0.279E+00                  0        10   0.102E+01
         2      0.141E+01    *  0.668E-01      1.00000
         2      0.141E+01       0.125E+00                  0        10   0.102E+01
         3      0.232E+00    *  0.210E-01      1.00000
         3      0.232E+00       0.320E-01                  0        10   0.101E+01
         4      0.147E-01    *  0.214E-02      1.00000
         4      0.147E-01       0.249E-02                  0        10   0.101E+01
         5      0.938E-04    *  0.166E-04      1.00000
         5      0.938E-04       0.168E-04                  0        10   0.101E+01
         6      0.440E-08    *  0.808E-09      1.00000
         6      0.440E-08       0.808E-09                  0        10   0.101E+01
         7      0.283E-15    *  0.386E-16      1.00000



 Solution of nonlinear system of equations obtained within   7 iteration steps

 Achieved relative accuracy 0.386E-16

   Subcondition ( 1,  10)  0.101E+01
   Sensitivity ( 1,  10)  0.100E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       7  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       7  ***
   ***  Function eval.    :       8  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.2651E-15

 problem: Chemical Equilibrium 1                          


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to        94 (  0.1 percent)

 Integer Workspace declared as      1000 is used up to        54 (  5.4 percent)



 N =   2

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     2
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.354E+05       0.353E+05                  0         2   0.200E+01
         1      0.350E+05    *  0.350E+05      0.01000
         1      0.350E+05       0.116E+03                  0         2   0.383E+01
         2      0.310E+05    *  0.103E+03      0.07334
         2      0.310E+05       0.508E+01                  0         2   0.285E+01
         3      0.173E+05    *  0.300E+01      0.29948
         3      0.173E+05       0.142E+01                  0         2   0.180E+01
         4      0.102E+05    *  0.468E+00      0.97042
         4      0.102E+05       0.205E+00                  0         2   0.129E+01
         5      0.169E+03    *  0.356E-02      1.00000
         5      0.169E+03       0.391E-02                  0         2   0.230E+01
         6      0.000E+00    *  0.000E+00      1.00000
         6      0.000E+00       0.000E+00                  0         2   0.200E+01
         7      0.000E+00    *  0.000E+00      1.00000



 Solution of nonlinear system of equations obtained within   7 iteration steps

 Achieved relative accuracy 0.000E+00

   Subcondition ( 1,   2)  0.200E+01
   Sensitivity ( 1,   2)  0.141E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       7  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       7  ***
   ***  Function eval.    :       8  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.0000E+00

 problem: Chemical Equilibrium 2                          


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       182 (  0.3 percent)

 Integer Workspace declared as      1000 is used up to        58 (  5.8 percent)



 N =   6

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     6
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.225E+16       0.219E+02                  0         6   0.232E+01
         1      0.222E+16    *  0.217E+02      0.01000
         1      0.222E+16       0.171E+02                  0         6   0.247E+01
         2      0.495E+15    *  0.823E-01      1.00000
         2      0.495E+15       0.177E+00                  0         6   0.562E+02
         3      0.124E+15    *  0.443E-01      1.00000
         3      0.124E+15       0.885E-01                  0         6   0.110E+03
         4      0.309E+14    *  0.221E-01      1.00000
         4      0.309E+14       0.443E-01                  0         6   0.110E+03
         5      0.773E+13    *  0.111E-01      1.00000
         5      0.773E+13       0.221E-01                  0         6   0.110E+03
         6      0.193E+13    *  0.553E-02      1.00000
         6      0.193E+13       0.111E-01                  0         6   0.110E+03
         7      0.483E+12    *  0.277E-02      1.00000
         7      0.483E+12       0.553E-02                  0         6   0.110E+03
         8      0.121E+12    *  0.138E-02      1.00000
         8      0.121E+12       0.277E-02                  0         6   0.110E+03
         9      0.302E+11    *  0.692E-03      1.00000
         9      0.302E+11       0.138E-02                  0         6   0.110E+03
        10      0.755E+10    *  0.346E-03      1.00000
        10      0.755E+10       0.691E-03                  0         6   0.110E+03
        11      0.189E+10    *  0.172E-03      1.00000
        11      0.189E+10       0.344E-03                  0         6   0.110E+03
        12      0.472E+09    *  0.855E-04      1.00000
        12      0.472E+09       0.170E-03                  0         6   0.111E+03
        13      0.118E+09    *  0.414E-04      1.00000
        13      0.118E+09       0.809E-04                  0         6   0.113E+03
        14      0.293E+08    *  0.185E-04      1.00000
        14      0.293E+08       0.344E-04                  0         6   0.120E+03
        15      0.682E+07    *  0.698E-05      1.00000
        15      0.682E+07       0.122E-04                  0         6   0.134E+03
        16      0.969E+06    *  0.183E-05      1.00000
        16      0.969E+06       0.244E-05                  0         6   0.141E+03
        17      0.215E+05    *  0.592E-07      1.00000
        17      0.215E+05       0.614E-07                  0         6   0.141E+03
        18      0.667E+01    *  0.201E-10      1.00000



 Solution of nonlinear system of equations obtained within  18 iteration steps

 Achieved relative accuracy 0.201E-10

   Subcondition ( 1,   6)  0.141E+03
   Sensitivity ( 1,   6)  0.173E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      18  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      18  ***
   ***  Function eval.    :      19  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.3433E-02

 problem: Chemical Equilibrium 3, Variante 1              


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       302 (  0.4 percent)

 Integer Workspace declared as      1000 is used up to        62 (  6.2 percent)



 N =  10

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    10
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.137E+02       0.199E+01                  0        10   0.374E+01
         1      0.136E+02    *  0.197E+01      0.01000
         1      0.136E+02       0.182E+01                  0        10   0.374E+01
         2      0.999E+01    *  0.122E+01      0.19268
         2      0.999E+01       0.637E+00                  0        10   0.366E+01
         3      0.376E+01    *  0.307E+00      0.67397
         3      0.376E+01       0.633E+00                  0        10   0.335E+01
         4      0.208E+01    *  0.347E+00      0.58766
         4      0.208E+01       0.614E+00                  0        10   0.356E+01
         5      0.975E+00    *  0.275E+00      0.77705
         5      0.975E+00       0.578E+00                  0        10   0.434E+01
         6      0.478E+00    *  0.265E+00      0.74471
         6      0.478E+00       0.462E+00                  0        10   0.608E+01
         7      0.227E+00    *  0.201E+00      0.79137
         7      0.227E+00       0.330E+00                  0        10   0.748E+01
         8      0.112E+00    *  0.144E+00      0.78519
         8      0.112E+00       0.250E+00                  0        10   0.105E+02
         9      0.585E-01    *  0.109E+00      0.78438
         9      0.585E-01       0.198E+00                  0        10   0.126E+02
        10      0.330E-01    *  0.871E-01      0.78143
        10      0.330E-01       0.162E+00                  0        10   0.132E+02
        11      0.203E-01    *  0.736E-01      0.76378
        11      0.203E-01       0.141E+00                  0        10   0.143E+02
        12      0.137E-01    *  0.675E-01      0.72577
        12      0.137E-01       0.118E+00                  0        10   0.157E+02
        13      0.101E-01    *  0.572E-01      0.71612
        13      0.101E-01       0.102E+00                  0        10   0.176E+02
        14      0.793E-02    *  0.500E-01      0.70601
        14      0.793E-02       0.908E-01                  0        10   0.195E+02
        15      0.637E-02    *  0.465E-01      0.67940
        15      0.637E-02       0.836E-01                  0        10   0.211E+02
        16      0.513E-02    *  0.438E-01      0.66188
        16      0.513E-02       0.776E-01                  0        10   0.220E+02
        17      0.401E-02    *  0.411E-01      0.65283
        17      0.401E-02       0.714E-01                  0        10   0.218E+02
        18      0.294E-02    *  0.374E-01      0.65925
        18      0.294E-02       0.639E-01                  0        10   0.221E+02
        19      0.195E-02    *  0.316E-01      0.69723
        19      0.195E-02       0.540E-01                  0        10   0.240E+02
        20      0.110E-02    *  0.234E-01      0.78206
        20      0.110E-02       0.405E-01                  0        10   0.253E+02
        21      0.419E-03    *  0.120E-01      0.97427
        21      0.419E-03       0.210E-01                  0        10   0.262E+02
        22      0.775E-04    *  0.371E-02      1.00000
        22      0.775E-04       0.508E-02                  0        10   0.267E+02
        23      0.396E-05    *  0.259E-03      1.00000
        23      0.396E-05       0.277E-03                  0        10   0.268E+02
        24      0.118E-07    *  0.823E-06      1.00000
        24      0.118E-07       0.821E-06                  0        10   0.268E+02
        25      0.105E-12    *  0.731E-11      1.00000



 Solution of nonlinear system of equations obtained within  25 iteration steps

 Achieved relative accuracy 0.731E-11

   Subcondition ( 1,  10)  0.268E+02
   Sensitivity ( 1,  10)  0.159E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      25  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      25  ***
   ***  Function eval.    :      26  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1382E-14

 problem: Chemical Equilibrium 3, Variante 2              


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       302 (  0.4 percent)

 Integer Workspace declared as      1000 is used up to        62 (  6.2 percent)



 N =  10

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    10
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.117E+02       0.626E+01                  0        10   0.374E+01
         1      0.116E+02    *  0.620E+01      0.01000
         1      0.116E+02       0.569E+01                  0        10   0.360E+01
         2      0.194E+01    *  0.146E+01      0.85390
         2      0.194E+01       0.434E+00                  0        10   0.324E+01
         3      0.522E-01    *  0.399E-01      1.00000
         3      0.522E-01       0.700E-01                  0        10   0.351E+01
         4      0.515E-02    *  0.164E-01      1.00000
         4      0.515E-02       0.352E-01                  0        10   0.368E+01
         5      0.123E-02    *  0.997E-02      1.00000
         5      0.123E-02       0.220E-01                  0        10   0.397E+01
         6      0.443E-03    *  0.638E-02      1.00000
         6      0.443E-03       0.133E-01                  0        10   0.446E+01
         7      0.170E-03    *  0.311E-02      1.00000
         7      0.170E-03       0.526E-02                  0        10   0.507E+01
         8      0.248E-04    *  0.546E-03      1.00000
         8      0.248E-04       0.672E-03                  0        10   0.552E+01
         9      0.370E-06    *  0.873E-05      1.00000
         9      0.370E-06       0.896E-05                  0        10   0.568E+01
        10      0.626E-10    *  0.146E-08      1.00000
        10      0.626E-10       0.146E-08                  0        10   0.568E+01
        11      0.901E-15    *  0.544E-16      1.00000



 Solution of nonlinear system of equations obtained within  11 iteration steps

 Achieved relative accuracy 0.544E-16

   Subcondition ( 1,  10)  0.568E+01
   Sensitivity ( 1,  10)  0.142E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      11  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      11  ***
   ***  Function eval.    :      12  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.9222E-15

 problem: SST pollution, 0 dim., sst1=360                 


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       134 (  0.2 percent)

 Integer Workspace declared as      1000 is used up to        56 (  5.6 percent)



 N =   4

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     4
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.193E+12       0.253E+01                  0         4   0.660E+07
         1      0.191E+12    *  0.250E+01      0.01000
         1      0.191E+12       0.249E+01                  0         4   0.647E+07
         2      0.930E+11    *  0.183E+01      0.51221
         2      0.930E+11       0.223E+01                  0         4   0.437E+07
         3      0.728E+11    *  0.185E+01      0.21754
         3      0.728E+11       0.151E+01                  0         4   0.162E+07
         4      0.415E+11    *  0.110E+01      0.43036
         4      0.415E+11       0.142E+01                  0         4   0.107E+07
         5      0.269E+11    *  0.108E+01      0.35145
         5      0.269E+11       0.130E+01                  0         4   0.550E+06
         6      0.162E+11    *  0.960E+00      0.39914
         6      0.162E+11       0.140E+01                  0         4   0.336E+06
         7      0.112E+11    *  0.110E+01      0.30732
         7      0.112E+11       0.155E+01                  0         4   0.192E+06
         8      0.821E+10    *  0.126E+01      0.26656
         8      0.821E+10       0.167E+01                  0         4   0.126E+06
         9      0.615E+10    *  0.138E+01      0.25114
         9      0.615E+10       0.174E+01                  0         4   0.965E+05
        10      0.461E+10    *  0.143E+01      0.25014
        10      0.461E+10       0.175E+01                  0         4   0.750E+05
        11      0.344E+10    *  0.143E+01      0.25378
        11      0.344E+10       0.172E+01                  0         4   0.570E+05
        12      0.254E+10    *  0.140E+01      0.26155
        12      0.254E+10       0.168E+01                  0         4   0.434E+05
        13      0.185E+10    *  0.135E+01      0.27157
        13      0.185E+10       0.161E+01                  0         4   0.351E+05
        14      0.133E+10    *  0.129E+01      0.28432
        14      0.133E+10       0.153E+01                  0         4   0.352E+05
        15      0.928E+09    *  0.121E+01      0.30058
        15      0.928E+09       0.144E+01                  0         4   0.463E+05
        16      0.628E+09    *  0.111E+01      0.32442
        16      0.628E+09       0.131E+01                  0         4   0.702E+05
        17      0.399E+09    *  0.966E+00      0.36719
        17      0.399E+09       0.110E+01                  0         4   0.111E+06
        18      0.217E+09    *  0.734E+00      0.46390
        18      0.217E+09       0.775E+00                  0         4   0.183E+06
        19      0.604E+08    *  0.342E+00      0.75857
        19      0.604E+08       0.293E+00                  0         4   0.334E+06
        20      0.658E+07    *  0.226E-01      1.00000
        20      0.658E+07       0.155E-01                  0         4   0.421E+06
        21      0.666E+04    *  0.114E-04      1.00000
        21      0.666E+04       0.984E-05                  0         4   0.438E+06
        22      0.147E-01    *  0.402E-10      1.00000



 Solution of nonlinear system of equations obtained within  22 iteration steps

 Achieved relative accuracy 0.402E-10

   Subcondition ( 1,   4)  0.438E+06
   Sensitivity ( 1,   4)  0.173E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      22  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      22  ***
   ***  Function eval.    :      23  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.6888E-07

 problem: SST pollution, 0 dim., sst1=3250                


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       134 (  0.2 percent)

 Integer Workspace declared as      1000 is used up to        56 (  5.6 percent)



 N =   4

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     4
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.193E+12       0.252E+01                  0         4   0.660E+07
         1      0.191E+12    *  0.250E+01      0.01000
         1      0.191E+12       0.249E+01                  0         4   0.647E+07
         2      0.931E+11    *  0.183E+01      0.51183
         2      0.931E+11       0.223E+01                  0         4   0.437E+07
         3      0.730E+11    *  0.185E+01      0.21609
         3      0.730E+11       0.150E+01                  0         4   0.162E+07
         4      0.418E+11    *  0.110E+01      0.42725
         4      0.418E+11       0.141E+01                  0         4   0.107E+07
         5      0.273E+11    *  0.107E+01      0.34658
         5      0.273E+11       0.127E+01                  0         4   0.554E+06
         6      0.164E+11    *  0.938E+00      0.39844
         6      0.164E+11       0.132E+01                  0         4   0.341E+06
         7      0.111E+11    *  0.103E+01      0.32243
         7      0.111E+11       0.142E+01                  0         4   0.196E+06
         8      0.794E+10    *  0.113E+01      0.28703
         8      0.794E+10       0.152E+01                  0         4   0.126E+06
         9      0.576E+10    *  0.122E+01      0.27379
         9      0.576E+10       0.156E+01                  0         4   0.882E+05
        10      0.416E+10    *  0.126E+01      0.27758
        10      0.416E+10       0.156E+01                  0         4   0.665E+05
        11      0.298E+10    *  0.125E+01      0.28502
        11      0.298E+10       0.153E+01                  0         4   0.485E+05
        12      0.210E+10    *  0.121E+01      0.29630
        12      0.210E+10       0.148E+01                  0         4   0.350E+05
        13      0.145E+10    *  0.116E+01      0.30837
        13      0.145E+10       0.142E+01                  0         4   0.267E+05
        14      0.982E+09    *  0.110E+01      0.32317
        14      0.982E+09       0.134E+01                  0         4   0.250E+05
        15      0.644E+09    *  0.102E+01      0.34469
        15      0.644E+09       0.123E+01                  0         4   0.317E+05
        16      0.396E+09    *  0.895E+00      0.38673
        16      0.396E+09       0.105E+01                  0         4   0.475E+05
        17      0.202E+09    *  0.675E+00      0.49103
        17      0.202E+09       0.717E+00                  0         4   0.774E+05
        18      0.330E+08    *  0.256E+00      0.85476
        18      0.330E+08       0.202E+00                  0         4   0.148E+06
        19      0.154E+07    *  0.328E-02      1.00000
        19      0.154E+07       0.252E-02                  0         4   0.173E+06
        20      0.321E+03    *  0.347E-05      1.00000
        20      0.321E+03       0.320E-05                  0         4   0.166E+06
        21      0.545E-03    *  0.882E-11      1.00000



 Solution of nonlinear system of equations obtained within  21 iteration steps

 Achieved relative accuracy 0.882E-11

   Subcondition ( 1,   4)  0.166E+06
   Sensitivity ( 1,   4)  0.173E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      21  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      21  ***
   ***  Function eval.    :      22  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.6016E-07

 problem: Distillation Column, Hydrocarbon-6              


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to      1309 (  1.9 percent)

 Integer Workspace declared as      1000 is used up to        81 (  8.1 percent)



 N =  29

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    29
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.491E+06       0.168E+00                  0        29   0.106E+03
         1      0.486E+06    *  0.167E+00      0.01000
         1      0.486E+06       0.163E+00                  0        29   0.106E+03
         2      0.294E+06    *  0.964E-01      0.30806
         2      0.294E+06       0.796E-01                  0        29   0.102E+03
         3      0.626E+05    *  0.176E-01      1.00000
         3      0.626E+05       0.142E-01                  0        29   0.112E+03
         4      0.185E+04    *  0.530E-03      1.00000
         4      0.185E+04       0.566E-03                  0        29   0.139E+03
         5      0.761E+00    *  0.735E-06      1.00000
         5      0.761E+00       0.734E-06                  0        29   0.139E+03
         6      0.109E-05    *  0.888E-12      1.00000



 Solution of nonlinear system of equations obtained within   6 iteration steps

 Achieved relative accuracy 0.888E-12

   Subcondition ( 1,  29)  0.139E+03
   Sensitivity ( 1,  29)  0.235E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       6  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       6  ***
   ***  Function eval.    :       7  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1304E-08

 problem: Distillation Column, Hydrocarbon-20             


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to     11249 ( 16.1 percent)

 Integer Workspace declared as      1000 is used up to       151 ( 15.1 percent)



 N =  99

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    99
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.229E+06       0.234E+00                  0        99   0.254E+03
         1      0.226E+06    *  0.232E+00      0.01000
         1      0.226E+06       0.229E+00                  0        99   0.255E+03
         2      0.108E+06    *  0.112E+00      0.43656
         2      0.108E+06       0.114E+00                  0        99   0.304E+03
         3      0.251E+05    *  0.345E-01      1.00000
         3      0.251E+05       0.714E-01                  0        99   0.118E+04
         4      0.120E+05    *  0.385E-01      1.00000
         4      0.120E+05       0.784E-01                  0        99   0.384E+04
         5      0.623E+04    *  0.375E-01      0.83386
         5      0.623E+04       0.741E-01                  0        99   0.103E+05
         6      0.580E+04    *  0.334E-01      0.85322
         6      0.580E+04       0.576E-01                  0        99   0.209E+05
         7      0.373E+04    *  0.153E-01      1.00000
         7      0.373E+04       0.144E-01                  0        99   0.266E+05
         8      0.296E+03    *  0.654E-03      1.00000
         8      0.296E+03       0.673E-03                  0        99   0.206E+05
         9      0.563E+00    *  0.147E-05      1.00000
         9      0.563E+00       0.148E-05                  0        99   0.204E+05
        10      0.241E-05    *  0.540E-11      1.00000



 Solution of nonlinear system of equations obtained within  10 iteration steps

 Achieved relative accuracy 0.540E-11

   Subcondition ( 1,  99)  0.204E+05
   Sensitivity ( 1,  99)  0.237E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      10  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      10  ***
   ***  Function eval.    :      11  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.2593E-08

 problem: Distillation Column, Hydrocarbon-40             


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to     42449 ( 60.6 percent)

 Integer Workspace declared as      1000 is used up to       251 ( 25.1 percent)



 N = 199

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =   199
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.238E+06       0.298E+00                  0       199   0.358E+03
         1      0.235E+06    *  0.295E+00      0.01000
         1      0.235E+06       0.291E+00                  0       199   0.358E+03
         2      0.131E+06    *  0.148E+00      0.40936
         2      0.131E+06       0.152E+00                  0       199   0.419E+03
         3      0.224E+05    *  0.331E-01      1.00000
         3      0.224E+05       0.497E-01                  0       199   0.222E+04
         4      0.725E+04    *  0.268E-01      1.00000
         4      0.725E+04       0.550E-01                  0       199   0.640E+04
         5      0.404E+04    *  0.272E-01      0.82525
         5      0.404E+04       0.554E-01                  0       199   0.171E+05
         6      0.349E+04    *  0.280E-01      0.74958
         6      0.349E+04       0.547E-01                  0       199   0.390E+05
         7      0.340E+04    *  0.277E-01      0.75057
         7      0.340E+04       0.544E-01                  0       199   0.947E+05
         8      0.330E+04    *  0.277E-01      0.74216
         8      0.330E+04       0.542E-01                  0       199   0.220E+06
         9      0.325E+04    *  0.277E-01      0.73998
         9      0.325E+04       0.540E-01                  0       199   0.517E+06
        10      0.322E+04    *  0.276E-01      0.73903
        10      0.322E+04       0.537E-01                  0       199   0.122E+07
        11      0.322E+04    *  0.274E-01      0.74107
        11      0.322E+04       0.533E-01                  0       199   0.276E+07
        12      0.324E+04    *  0.270E-01      0.74672
        12      0.324E+04       0.525E-01                  0       199   0.615E+07
        13      0.320E+04    *  0.263E-01      0.75664
        13      0.320E+04       0.498E-01                  0       199   0.143E+08
        14      0.274E+04    *  0.219E-01      0.83072
        14      0.274E+04       0.310E-01                  0       199   0.264E+08
        15      0.156E+04    *  0.561E-02      1.00000
        15      0.156E+04       0.559E-02                  0       199   0.232E+08
        16      0.698E+02    *  0.236E-03      1.00000
        16      0.698E+02       0.240E-03                  0       199   0.211E+08
        17      0.102E+00    *  0.281E-06      1.00000
        17      0.102E+00       0.281E-06                  0       199   0.210E+08
        18      0.138E-06    *  0.200E-10      1.00000



 Solution of nonlinear system of equations obtained within  18 iteration steps

 Achieved relative accuracy 0.200E-10

   Subcondition ( 1, 199)  0.210E+08
   Sensitivity ( 1, 199)  0.237E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      18  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      18  ***
   ***  Function eval.    :      19  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.2482E-08

 problem: Distillation Column, Methanol-8                 


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to      1457 (  2.1 percent)

 Integer Workspace declared as      1000 is used up to        83 (  8.3 percent)



 N =  31

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    31
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.118E+07       0.385E-01                  0        31   0.746E+02
         1      0.117E+07    *  0.381E-01      0.01000
         1      0.117E+07       0.381E-01                  0        31   0.744E+02
         2      0.157E+06    *  0.563E-02      1.00000
         2      0.157E+06       0.457E-02                  0        31   0.554E+02
         3      0.380E+04    *  0.910E-04      1.00000
         3      0.380E+04       0.964E-04                  0        31   0.556E+02
         4      0.229E+01    *  0.634E-07      1.00000
         4      0.229E+01       0.637E-07                  0        31   0.555E+02
         5      0.834E-06    *  0.314E-13      1.00000



 Solution of nonlinear system of equations obtained within   5 iteration steps

 Achieved relative accuracy 0.314E-13

   Subcondition ( 1,  31)  0.555E+02
   Sensitivity ( 1,  31)  0.250E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       5  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       5  ***
   ***  Function eval.    :       6  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.4434E-08

 problem: Semiconductor boundary condition                


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       182 (  0.3 percent)

 Integer Workspace declared as      1000 is used up to        58 (  5.8 percent)



 N =   6

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     6
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.473E+07       0.612E+05                  0         6   0.137E+01
         0 FCN could not be evaluated          0.01000     0         6
         0 FCN could not be evaluated          0.00500     0         6
         0 FCN could not be evaluated          0.00250     0         6
         0 FCN could not be evaluated          0.00125     0         6
         0 FCN could not be evaluated          0.00063     0         6
         0 FCN could not be evaluated          0.00031     0         6
         0 FCN could not be evaluated          0.00016     0         6
         0 FCN could not be evaluated          0.00010     0         6

 Newton method fails to converge

   Subcondition ( 1,   6)  0.137E+01
   Sensitivity ( 1,   6)  0.112E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       1  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       1  ***
   ***  Function eval.    :       9  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.4440E-12

 problem: Exponentional/Sine function                     


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to        94 (  0.1 percent)

 Integer Workspace declared as      1000 is used up to        54 (  5.4 percent)



 N =   2

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     2
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.193E+01       0.443E+03                  0         2   0.164E+03
         0      0.301E+18    *  0.565E+19      0.01000
         1      0.193E+01    *  0.443E+03      0.00010
         1      0.193E+01       0.445E+02                  0         2   0.175E+02
         2      0.193E+01    *  0.444E+02      0.00111
         2      0.193E+01       0.200E+02                  0         2   0.924E+01
         3      0.194E+01    *  0.198E+02      0.00447
         3      0.194E+01       0.758E+01                  0         2   0.528E+01
         4      0.192E+01    *  0.732E+01      0.01893
         4      0.192E+01       0.223E+01                  0         2   0.348E+01
         5      0.177E+01    *  0.185E+01      0.09090
         5      0.177E+01       0.545E+00                  0         2   0.266E+01
         6      0.103E+01    *  0.125E+00      0.50993
         6      0.103E+01       0.865E-01                  0         2   0.200E+01
         7      0.174E+00    *  0.140E-01      1.00000
         7      0.174E+00       0.219E-01                  0         2   0.178E+01
         8      0.854E-02    *  0.121E-02      1.00000
         8      0.854E-02       0.142E-02                  0         2   0.165E+01
         9      0.258E-04    *  0.460E-05      1.00000
         9      0.258E-04       0.466E-05                  0         2   0.162E+01
        10      0.250E-09    *  0.460E-10      1.00000



 Solution of nonlinear system of equations obtained within  10 iteration steps

 Achieved relative accuracy 0.460E-10

   Subcondition ( 1,   2)  0.162E+01
   Sensitivity ( 1,   2)  0.141E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      10  ***
   ***  Corrector steps   :       1  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      10  ***
   ***  Function eval.    :      12  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.5341E-14

 problem: Combustion of Propane - Full formulation        


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       337 (  0.5 percent)

 Integer Workspace declared as      1000 is used up to        63 (  6.3 percent)



 N =  11

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    11
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.120E+23       0.220E+02                  0        11   0.670E+15
         0 FCN could not be evaluated          0.01000     0        11
         0 FCN could not be evaluated          0.00500     0        11
         0 FCN could not be evaluated          0.00250     0        11
         0 FCN could not be evaluated          0.00125     0        11
         0 FCN could not be evaluated          0.00063     0        11
         0 FCN could not be evaluated          0.00031     0        11
         0 FCN could not be evaluated          0.00016     0        11
         0 FCN could not be evaluated          0.00010     0        11

 Newton method fails to converge

   Subcondition ( 1,  11)  0.670E+15
   Sensitivity ( 1,  11)  0.245E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       1  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       1  ***
   ***  Function eval.    :       9  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.6839E-12

 problem: Combustion of Propane - Reduced formulation     


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       157 (  0.2 percent)

 Integer Workspace declared as      1000 is used up to        57 (  5.7 percent)



 N =   5

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     5
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.804E+00       0.153E+02                  0         5   0.319E+04
         1      0.795E+00    *  0.151E+02      0.01000
         1      0.795E+00       0.108E+02                  0         5   0.413E+04
         2      0.757E+00    *  0.103E+02      0.04779
         2      0.757E+00       0.114E+02                  0         5   0.987E+04
         3      0.403E+00    *  0.785E+01      0.43508
         3      0.403E+00       0.169E+01                  0         5   0.260E+05
         4      0.109E+00    *  0.150E+00      0.82612
         4      0.109E+00       0.138E-01                  0         5   0.272E+05
         5      0.295E-01    *  0.731E-02      1.00000
         5      0.295E-01       0.344E-02                  0         5   0.172E+05
         6      0.106E-02    *  0.160E-03      1.00000
         6      0.106E-02       0.173E-03                  0         5   0.176E+05
         7      0.155E-05    *  0.287E-06      1.00000
         7      0.155E-05       0.287E-06                  0         5   0.176E+05
         8      0.326E-11    *  0.648E-12      1.00000



 Solution of nonlinear system of equations obtained within   8 iteration steps

 Achieved relative accuracy 0.648E-12

   Subcondition ( 1,   5)  0.176E+05
   Sensitivity ( 1,   5)  0.173E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       8  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       8  ***
   ***  Function eval.    :       9  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.3565E-16

 problem: Intersection of an Ellipsoid with a Hyperboloid 


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       209 (  0.3 percent)

 Integer Workspace declared as      1000 is used up to        59 (  5.9 percent)



 N =   7

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     7
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.471E+01       0.144E+01                  0         7   0.128E+02
         1      0.466E+01    *  0.143E+01      0.01000
         1      0.466E+01       0.139E+01                  0         7   0.126E+02
         2      0.347E+00    *  0.689E-01      0.63516
         2      0.347E+00       0.109E+00                  0         7   0.227E+02
         3      0.169E+00    *  0.729E-01      1.00000
         3      0.169E+00       0.396E-01                  0         7   0.102E+02
         4      0.226E-01    *  0.113E-02      1.00000
         4      0.226E-01       0.115E-02                  0         7   0.132E+02
         5      0.241E-04    *  0.125E-05      1.00000
         5      0.241E-04       0.128E-05                  0         7   0.125E+02
         6      0.288E-10    *  0.152E-11      1.00000



 Solution of nonlinear system of equations obtained within   6 iteration steps

 Achieved relative accuracy 0.152E-11

   Subcondition ( 1,   7)  0.125E+02
   Sensitivity ( 1,   7)  0.174E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       6  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       6  ***
   ***  Function eval.    :       7  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1056E-14

 problem: Esterification Reaction                         


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       269 (  0.4 percent)

 Integer Workspace declared as      1000 is used up to        61 (  6.1 percent)



 N =   9

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     9
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.458E+07       0.250E+01                  0         9   0.124E+04
         1      0.453E+07    *  0.248E+01      0.01000
         1      0.453E+07       0.256E+01                  0         9   0.128E+04
         2      0.113E+07    *  0.913E+00      1.00000
         2      0.113E+07       0.158E+01                  0         9   0.353E+03
         3      0.560E+06    *  0.974E+00      0.59395
         3      0.560E+06       0.358E+00                  0         9   0.790E+02
         4      0.140E+06    *  0.976E-01      1.00000
         4      0.140E+06       0.262E+00                  0         9   0.696E+02
         5      0.350E+05    *  0.747E-01      1.00000
         5      0.350E+05       0.250E+00                  0         9   0.490E+02
         6      0.875E+04    *  0.747E-01      1.00000
         6      0.875E+04       0.270E+00                  0         9   0.418E+02
         7      0.219E+04    *  0.906E-01      1.00000
         7      0.219E+04       0.372E+00                  0         9   0.341E+02
         8      0.844E+03    *  0.193E+00      0.75926
         8      0.844E+03       0.111E+01                  0         9   0.746E+02
         9      0.721E+03    *  0.978E+00      0.15329
         9      0.721E+03       0.452E+01                  0         9   0.230E+03
        10      0.701E+03    *  0.440E+01      0.02846
        10      0.701E+03       0.222E+02                  0         9   0.870E+03
        11      0.699E+03    *  0.221E+02      0.00292
        11      0.699E+03       0.169E+03                  0         9   0.574E+04
        12      0.699E+03    *  0.169E+03      0.00010
        12      0.699E+03       0.198E+04                  0         9   0.635E+05
        12      0.699E+03    *  0.200E+04      0.00010
        12 Not accepted damping factor        0.00000     0         9
        12      0.699E+03       0.323E+00                  0         8   0.240E+02
        13      0.200E+03    *  0.705E-01      1.00000
        13      0.200E+03       0.318E+02                  0         9   0.798E+03
        14      0.200E+03    *  0.318E+02      0.00010
        14      0.200E+03       0.332E+02                  0         9   0.776E+03
        15      0.199E+03    *  0.332E+02      0.00401
        15      0.199E+03       0.337E+03                  0         9   0.859E+04
        16      0.199E+03    *  0.337E+03      0.00010
        16      0.199E+03       0.125E+03                  0         9   0.339E+04
        16      0.199E+03    *  0.125E+03      0.00018
        16      0.199E+03    *  0.125E+03      0.00010
        16 Not accepted damping factor        0.00003     0         9
        16      0.199E+03       0.366E+00                  0         8   0.524E+02
        17      0.126E+03    *  0.170E-01      1.00000
        17      0.126E+03       0.791E+03                  0         9   0.234E+05
        17      0.127E+03    *  0.792E+03      0.00010
        17 Not accepted damping factor        0.00000     0         9
        17      0.126E+03       0.259E+00                  0         8   0.527E+02
        18      0.121E+03    *  0.178E+00      0.30842
        18      0.121E+03       0.111E+04                  0         9   0.270E+05
        18      0.122E+03    *  0.111E+04      0.00010
        18 Not accepted damping factor        0.00000     0         9
        18      0.121E+03       0.172E+00                  0         8   0.132E+03
        19      0.110E+03    *  0.202E-01      1.00000
        19      0.110E+03       0.292E+04                  0         9   0.638E+05
        19      0.113E+03    *  0.294E+04      0.00010
        19 Not accepted damping factor        0.00000     0         9
        19      0.110E+03       0.126E-01                  0         8   0.133E+03
        20      0.108E+03    *  0.708E-03      1.00000
        20      0.108E+03       0.598E+03                  0         9   0.171E+05
        21      0.108E+03    *  0.598E+03      0.00010
        21      0.108E+03       0.686E+02                  0         9   0.189E+04
        21      0.108E+03    *  0.686E+02      0.00081
        22      0.108E+03    *  0.686E+02      0.00038
        22      0.108E+03       0.129E+03                  0         9   0.347E+04
        22      0.108E+03    *  0.129E+03      0.00015
        22      0.108E+03    *  0.129E+03      0.00010
        22 Not accepted damping factor        0.00004     0         9
        22      0.108E+03       0.598E-02                  0         8   0.895E+02
        23      0.108E+03    *  0.114E-03      1.00000
        23      0.108E+03       0.139E+03                  0         9   0.380E+04
        23      0.108E+03    *  0.139E+03      0.00010
        23 Not accepted damping factor        0.00003     0         9
        23      0.108E+03       0.184E-02                  0         8   0.916E+02
        24      0.108E+03    *  0.150E-02      0.18726
        24      0.108E+03       0.139E+03                  0         9   0.379E+04
        24      0.108E+03    *  0.139E+03      0.00010
        24 Not accepted damping factor        0.00003     0         9
        24      0.108E+03       0.156E-02                  0         8   0.916E+02
        25      0.108E+03    *  0.239E-04      1.00000
        25      0.108E+03       0.137E+03                  0         9   0.374E+04
        25      0.108E+03    *  0.137E+03      0.00010
        25 Not accepted damping factor        0.00003     0         9
        25      0.108E+03       0.118E-03                  0         8   0.916E+02
        26      0.108E+03    *  0.987E-06      1.00000
        26      0.108E+03       0.137E+03                  0         9   0.375E+04
        26      0.108E+03    *  0.137E+03      0.00010
        26 Not accepted damping factor        0.00003     0         9
        26      0.108E+03       0.508E-04                  0         8   0.916E+02
        27      0.108E+03    *  0.486E-04      0.04443
        27      0.108E+03       0.137E+03                  0         9   0.375E+04
        27      0.108E+03    *  0.137E+03      0.00010
        27 Not accepted damping factor        0.00003     0         9
        27      0.108E+03       0.512E-04                  0         8   0.916E+02
        28      0.108E+03    *  0.140E-04      0.73338
        28      0.108E+03       0.137E+03                  0         9   0.375E+04
        28      0.108E+03    *  0.137E+03      0.00010
        28 Not accepted damping factor        0.00003     0         9
        28      0.108E+03       0.183E-04                  0         8   0.916E+02
        29      0.108E+03    *  0.188E-06      1.00000
        29      0.108E+03       0.137E+03                  0         9   0.375E+04
        29      0.108E+03    *  0.137E+03      0.00010
        29 Not accepted damping factor        0.00003     0         9
        29      0.108E+03       0.397E-05                  0         8   0.916E+02
        30      0.108E+03    *  0.308E-05      0.22516
        30      0.108E+03       0.137E+03                  0         9   0.375E+04
        30      0.108E+03    *  0.137E+03      0.00010
        30 Not accepted damping factor        0.00003     0         9
        30      0.108E+03       0.348E-05                  0         8   0.916E+02
        31      0.108E+03    *  0.295E-07      1.00000
        31      0.108E+03       0.137E+03                  0         9   0.375E+04
        31      0.108E+03    *  0.137E+03      0.00010
        31 Not accepted damping factor        0.00003     0         9
        31      0.108E+03       0.849E-06                  0         8   0.916E+02
        32      0.108E+03    *  0.727E-06      0.14524
        32      0.108E+03       0.137E+03                  0         9   0.375E+04
        32      0.108E+03    *  0.137E+03      0.00010
        32 Not accepted damping factor        0.00003     0         9
        32      0.108E+03       0.796E-06                  0         8   0.916E+02
        33      0.108E+03    *  0.680E-08      1.00000
        33      0.108E+03       0.137E+03                  0         9   0.375E+04
        33      0.108E+03    *  0.137E+03      0.00010
        33 Not accepted damping factor        0.00003     0         9
        33      0.108E+03       0.195E-06                  0         8   0.916E+02
        34      0.108E+03    *  0.167E-06      0.14496
        34      0.108E+03       0.137E+03                  0         9   0.375E+04
        34      0.108E+03    *  0.137E+03      0.00010
        34 Not accepted damping factor        0.00003     0         9
        34      0.108E+03       0.183E-06                  0         8   0.916E+02
        35      0.108E+03    *  0.156E-08      1.00000
        35      0.108E+03       0.137E+03                  0         9   0.375E+04
        35      0.108E+03    *  0.137E+03      0.00010
        35 Not accepted damping factor        0.00003     0         9
        35      0.108E+03       0.450E-07                  0         8   0.916E+02
        36      0.108E+03    *  0.385E-07      0.14439
        36      0.108E+03       0.137E+03                  0         9   0.375E+04
        36      0.108E+03    *  0.137E+03      0.00010
        36 Not accepted damping factor        0.00003     0         9
        36      0.108E+03       0.422E-07                  0         8   0.916E+02
        37      0.108E+03    *  0.360E-09      1.00000
        37      0.108E+03       0.137E+03                  0         9   0.375E+04
        37      0.108E+03    *  0.137E+03      0.00010
        37 Not accepted damping factor        0.00003     0         9
        37      0.108E+03       0.104E-07                  0         8   0.916E+02
        38      0.108E+03    *  0.887E-08      0.14463
        38      0.108E+03       0.137E+03                  0         9   0.375E+04
        38      0.108E+03    *  0.137E+03      0.00010
        38 Not accepted damping factor        0.00003     0         9
        38      0.108E+03       0.972E-08                  0         8   0.916E+02

 Iteration terminates at stationary point


   Subcondition ( 1,   8)  0.916E+02
   Sensitivity ( 1,   8)  0.240E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      39  ***
   ***  Corrector steps   :      25  ***
   ***  Rejected rk-1 st. :      22  ***
   ***  Jacobian eval.    :      39  ***
   ***  Function eval.    :      65  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1076E+03

 problem: Flow in a Driven Cavity                         


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to     11462 ( 16.4 percent)

 Integer Workspace declared as      1000 is used up to       152 ( 15.2 percent)



 N = 100

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =   100
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.656E-01       0.255E-01                  0       100   0.714E+02
         1      0.650E-01    *  0.253E-01      0.01000
         1      0.650E-01       0.253E-01                  0       100   0.706E+02
         2      0.153E-01    *  0.249E-02      1.00000
         2      0.153E-01       0.553E-02                  0       100   0.569E+02
         3      0.189E-02    *  0.952E-03      1.00000
         3      0.189E-02       0.834E-03                  0       100   0.676E+02
         4      0.892E-04    *  0.486E-04      1.00000
         4      0.892E-04       0.471E-04                  0       100   0.446E+02
         5      0.192E-06    *  0.859E-07      1.00000
         5      0.192E-06       0.857E-07                  0       100   0.446E+02
         6      0.703E-12    *  0.323E-12      1.00000



 Solution of nonlinear system of equations obtained within   6 iteration steps

 Achieved relative accuracy 0.323E-12

   Subcondition ( 1, 100)  0.446E+02
   Sensitivity ( 1, 100)  0.144E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       6  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       6  ***
   ***  Function eval.    :       7  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.3109E-16

 problem: Flow in a Channel                               


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to      7582 ( 10.8 percent)

 Integer Workspace declared as      1000 is used up to       132 ( 13.2 percent)



 N =  80

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    80
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.454E+03       0.664E+03                  0        80   0.508E+06
         1      0.450E+03    *  0.658E+03      0.01000
         1      0.450E+03       0.815E+02                  0        80   0.465E+06
         2      0.393E+03    *  0.724E+02      0.12333
         2      0.393E+03       0.174E+02                  0        80   0.393E+06
         3      0.385E+03    *  0.171E+02      0.01987
         3      0.385E+03       0.133E+02                  0        80   0.370E+06
         4      0.337E+03    *  0.123E+02      0.12365
         4      0.337E+03       0.487E+01                  0        80   0.330E+06
         5      0.291E+03    *  0.432E+01      0.13571
         5      0.291E+03       0.229E+01                  0        80   0.294E+06
         6      0.173E+03    *  0.159E+01      0.39694
         6      0.173E+03       0.972E+00                  0        80   0.208E+06
         7      0.465E+02    *  0.374E+00      0.73563
         7      0.465E+02       0.334E+00                  0        80   0.140E+06
         8      0.131E+01    *  0.111E-01      1.00000
         8      0.131E+01       0.938E-02                  0        80   0.117E+06
         9      0.221E-03    *  0.456E-05      1.00000
         9      0.221E-03       0.367E-05                  0        80   0.117E+06
        10      0.552E-11    *  0.172E-11      1.00000



 Solution of nonlinear system of equations obtained within  10 iteration steps

 Achieved relative accuracy 0.172E-11

   Subcondition ( 1,  80)  0.117E+06
   Sensitivity ( 1,  80)  0.245E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      10  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      10  ***
   ***  Function eval.    :      11  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.9221E-13

 problem: Gupta Problem                                   


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       374 (  0.5 percent)

 Integer Workspace declared as      1000 is used up to        64 (  6.4 percent)



 N =  12

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    12
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.724E+03       0.217E+01                  0        12   0.271E+02
         1      0.717E+03    *  0.215E+01      0.01000
         1      0.717E+03       0.210E+01                  0        12   0.277E+02
         1 FCN could not be evaluated          0.98393     0        12
         1 FCN could not be evaluated          0.49196     0        12
         1 FCN could not be evaluated          0.24598     0        12
         2      0.629E+03    *  0.185E+01      0.12299
         2      0.629E+03       0.139E+01                  0        12   0.115E+03
         2 FCN could not be evaluated          1.00000     0        12
         3      0.315E+03    *  0.878E+00      0.50000
         3      0.315E+03       0.563E+00                  0        12   0.192E+04
         4      0.315E+02    *  0.121E+00      0.90000
         4      0.315E+02       0.789E-01                  0        12   0.471E+03
         5      0.671E+00    *  0.282E-01      1.00000
         5      0.671E+00       0.298E-01                  0        12   0.117E+04
         6      0.374E-01    *  0.870E-05      1.00000
         6      0.374E-01       0.871E-05                  0        12   0.938E+03
         7      0.180E-03    *  0.176E-07      1.00000
         7      0.180E-03       0.177E-07                  0        12   0.975E+03
         8      0.420E-08    *  0.414E-12      1.00000



 Solution of nonlinear system of equations obtained within   8 iteration steps

 Achieved relative accuracy 0.414E-12

   Subcondition ( 1,  12)  0.975E+03
   Sensitivity ( 1,  12)  0.173E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       8  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       8  ***
   ***  Function eval.    :      13  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1956E-12

 problem: Human Heart Dipole                              


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       238 (  0.3 percent)

 Integer Workspace declared as      1000 is used up to        60 (  6.0 percent)



 N =   8

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     8
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.154E+00       0.115E+01                  0         8   0.288E+01
         1      0.153E+00    *  0.114E+01      0.01000
         1      0.153E+00       0.114E+01                  0         8   0.288E+01
         2      0.100E+00    *  0.701E+00      0.26769
         2      0.100E+00       0.119E+01                  0         8   0.108E+02
         3      0.223E+00    *  0.116E+01      0.22464
         3      0.223E+00       0.441E+01                  0         8   0.111E+03
         4      0.219E+00    *  0.433E+01      0.01352
         4      0.219E+00       0.293E+01                  0         8   0.700E+02
         5      0.220E+00    *  0.269E+01      0.05996
         5      0.220E+00       0.121E+01                  0         8   0.379E+02
         6      0.372E+00    *  0.813E+00      0.21389
         6      0.372E+00       0.315E+00                  0         8   0.257E+02
         7      0.316E+00    *  0.136E+00      0.95796
         7      0.316E+00       0.567E+00                  0         8   0.312E+02
         8      0.292E+00    *  0.481E+00      0.16301
         8      0.292E+00       0.589E+00                  0         8   0.395E+02
         9      0.380E+00    *  0.365E+00      0.40921
         9      0.380E+00       0.392E+00                  0         8   0.354E+02
        10      0.500E+00    *  0.373E+00      1.00000
        10      0.500E+00       0.226E+00                  0         8   0.166E+02
        11      0.305E+00    *  0.124E+00      1.00000
        11      0.305E+00       0.524E+00                  0         8   0.236E+02
        12      0.232E+00    *  0.421E+00      0.25096
        12      0.232E+00       0.659E+00                  0         8   0.215E+02
        13      0.145E+00    *  0.548E+00      0.33324
        13      0.145E+00       0.170E+01                  0         8   0.254E+02
        14      0.145E+00    *  0.162E+01      0.06613
        14      0.145E+00       0.351E+01                  0         8   0.316E+02
        15      0.147E+00    *  0.344E+01      0.02314
        15      0.147E+00       0.425E+01                  0         8   0.379E+02
        16      0.132E+00    *  0.413E+01      0.02339
        16      0.132E+00       0.315E+01                  0         8   0.465E+02
        17      0.107E+00    *  0.297E+01      0.04669
        17      0.107E+00       0.212E+01                  0         8   0.839E+02
        17      0.284E+00    *  0.226E+01      0.10170
        18      0.111E+00    *  0.207E+01      0.02859
        18      0.111E+00       0.190E+01                  0         8   0.125E+03
        18      0.284E+00    *  0.212E+01      0.11131
        19      0.113E+00    *  0.186E+01      0.02676
        19      0.113E+00       0.182E+01                  0         8   0.185E+03
        19      0.271E+00    *  0.206E+01      0.11674
        20      0.114E+00    *  0.178E+01      0.02485
        20      0.114E+00       0.178E+01                  0         8   0.265E+03
        20      0.206E+00    *  0.190E+01      0.10416
        21      0.114E+00    *  0.174E+01      0.02551
        21      0.114E+00       0.184E+01                  0         8   0.382E+03
        21      0.153E+00    *  0.185E+01      0.08470
        22      0.115E+00    *  0.179E+01      0.02876
        22      0.115E+00       0.206E+01                  0         8   0.566E+03
        23      0.128E+00    *  0.201E+01      0.06739
        23      0.128E+00       0.366E+01                  0         8   0.128E+04
        24      0.125E+00    *  0.355E+01      0.02944
        24      0.125E+00       0.350E+01                  0         8   0.137E+04
        24      0.393E+01    *  0.528E+02      0.55810
        25      0.123E+00    *  0.346E+01      0.01004
        25      0.123E+00       0.346E+01                  0         8   0.140E+04
        25      0.122E+02    *  0.189E+03      0.80163
        26      0.123E+00    *  0.344E+01      0.00586
        26      0.123E+00       0.343E+01                  0         8   0.141E+04
        26      0.148E+02    *  0.234E+03      0.86099
        27      0.122E+00    *  0.341E+01      0.00543
        27      0.122E+00       0.341E+01                  0         8   0.142E+04
        27      0.147E+02    *  0.233E+03      0.87292
        28      0.121E+00    *  0.339E+01      0.00556
        28      0.121E+00       0.338E+01                  0         8   0.144E+04
        28      0.140E+02    *  0.222E+03      0.87461
        29      0.121E+00    *  0.336E+01      0.00581
        29      0.121E+00       0.336E+01                  0         8   0.145E+04
        29      0.131E+02    *  0.207E+03      0.87090
        30      0.120E+00    *  0.333E+01      0.00615
        30      0.120E+00       0.333E+01                  0         8   0.146E+04
        30      0.119E+02    *  0.187E+03      0.86217
        31      0.119E+00    *  0.330E+01      0.00659
        31      0.119E+00       0.329E+01                  0         8   0.148E+04
        31      0.105E+02    *  0.165E+03      0.84830
        32      0.118E+00    *  0.327E+01      0.00716
        32      0.118E+00       0.325E+01                  0         8   0.150E+04
        32      0.910E+01    *  0.141E+03      0.82912
        33      0.117E+00    *  0.323E+01      0.00790
        33      0.117E+00       0.321E+01                  0         8   0.153E+04
        33      0.762E+01    *  0.116E+03      0.80457
        34      0.116E+00    *  0.318E+01      0.00887
        34      0.116E+00       0.314E+01                  0         8   0.156E+04
        34      0.611E+01    *  0.910E+02      0.77216
        35      0.115E+00    *  0.311E+01      0.01020
        35      0.115E+00       0.283E+01                  0         8   0.148E+04
        35      0.411E+01    *  0.543E+02      0.70530
        36      0.114E+00    *  0.280E+01      0.01279
        36      0.114E+00       0.252E+01                  0         8   0.142E+04
        36      0.262E+01    *  0.300E+02      0.63689
        37      0.112E+00    *  0.248E+01      0.01657
        37      0.112E+00       0.220E+01                  0         8   0.135E+04
        37      0.161E+01    *  0.154E+02      0.57220
        38      0.109E+00    *  0.215E+01      0.02204
        38      0.109E+00       0.185E+01                  0         8   0.128E+04
        38      0.954E+00    *  0.726E+01      0.51510
        39      0.105E+00    *  0.179E+01      0.03013
        39      0.105E+00       0.150E+01                  0         8   0.120E+04
        39      0.561E+00    *  0.315E+01      0.47144
        40      0.100E+00    *  0.143E+01      0.04225
        40      0.100E+00       0.115E+01                  0         8   0.110E+04
        40      0.343E+00    *  0.130E+01      0.45074
        41      0.925E-01    *  0.107E+01      0.06075
        41      0.925E-01       0.831E+00                  0         8   0.989E+03
        42      0.239E+00    *  0.590E+00      0.46992
        42      0.239E+00       0.142E+00                  0         8   0.692E+03
        43      0.101E+00    *  0.611E-01      1.00000
        43      0.101E+00       0.188E+00                  0         8   0.196E+04
        44      0.158E-01    *  0.491E-01      1.00000
        44      0.158E-01       0.757E-01                  0         8   0.185E+04
        45      0.207E-02    *  0.274E-02      1.00000
        45      0.207E-02       0.324E-02                  0         8   0.174E+04
        46      0.190E-04    *  0.141E-04      1.00000
        46      0.190E-04       0.146E-04                  0         8   0.174E+04
        47      0.961E-09    *  0.209E-08      1.00000
        47      0.961E-09       0.210E-08                  0         8   0.174E+04
        48      0.226E-16    *  0.188E-17      1.00000



 Solution of nonlinear system of equations obtained within  48 iteration steps

 Achieved relative accuracy 0.188E-17

   Subcondition ( 1,   8)  0.174E+04
   Sensitivity ( 1,   8)  0.200E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      48  ***
   ***  Corrector steps   :      22  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      48  ***
   ***  Function eval.    :      71  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.7666E-19

 problem: Incompressible Elastic Rods                     


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to      7238 ( 10.3 percent)

 Integer Workspace declared as      1000 is used up to       130 ( 13.0 percent)



 N =  78

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    78
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.103E+00       0.346E+00                  0        77   0.767E+03
         1      0.102E+00    *  0.343E+00      0.01000
         1      0.102E+00       0.544E+07                  0        78   0.661E+09
         2      0.617E+02    *  0.535E+07      0.00010
         2      0.617E+02       0.546E+03                  0        78   0.749E+07
         3      0.608E+02    *  0.546E+03      0.00120
         3      0.608E+02       0.208E+05                  0        78   0.202E+09
         3      0.608E+02    *  0.208E+05      0.00010
         3 Not accepted damping factor        0.00002     0        78
         3      0.608E+02       0.504E+02                  0        77   0.282E+07
         3      0.427E+03    *  0.853E+02      0.03805
         4      0.605E+02    *  0.503E+02      0.00086
         4      0.605E+02       0.802E+03                  0        78   0.809E+07
         5      0.605E+02    *  0.802E+03      0.00010
         5      0.605E+02       0.711E+03                  0        78   0.839E+07
         6      0.604E+02    *  0.710E+03      0.00064
         6      0.604E+02       0.274E+03                  0        78   0.678E+07
         7      0.586E+02    *  0.273E+03      0.00276
         7      0.586E+02       0.114E+04                  0        78   0.112E+08
         8      0.586E+02    *  0.114E+04      0.00010
         8      0.586E+02       0.130E+04                  0        78   0.124E+08
         9      0.586E+02    *  0.130E+04      0.00038
         9      0.586E+02       0.188E+04                  0        78   0.101E+08
        10      0.586E+02    *  0.188E+04      0.00010
        10      0.586E+02       0.320E+04                  0        78   0.161E+08
        11      0.585E+02    *  0.320E+04      0.00014
        11      0.585E+02       0.402E+04                  0        78   0.147E+08
        12      0.585E+02    *  0.402E+04      0.00010
        12      0.585E+02       0.413E+04                  0        78   0.578E+07
        13      0.585E+02    *  0.413E+04      0.00010
        13      0.585E+02       0.873E+04                  0        78   0.120E+08
        14      0.585E+02    *  0.873E+04      0.00010
        14      0.585E+02       0.217E+05                  0        78   0.721E+07
        15      0.585E+02    *  0.217E+05      0.00010
        15      0.585E+02       0.151E+04                  0        78   0.476E+06
        16      0.585E+02    *  0.151E+04      0.00062
        16      0.585E+02       0.368E+03                  0        78   0.358E+06
        17      0.583E+02    *  0.366E+03      0.00250
        17      0.583E+02       0.977E+02                  0        78   0.240E+06
        18      0.578E+02    *  0.962E+02      0.01003
        18      0.578E+02       0.254E+02                  0        78   0.136E+06
        19      0.552E+02    *  0.239E+02      0.03997
        19      0.552E+02       0.713E+01                  0        78   0.733E+05
        20      0.427E+02    *  0.595E+01      0.14280
        20      0.427E+02       0.565E+01                  0        78   0.148E+06
        21      0.389E+02    *  0.532E+01      0.06561
        21      0.389E+02       0.694E+01                  0        78   0.195E+06
        22      0.257E+02    *  0.574E+01      0.24029
        22      0.257E+02       0.643E+02                  0        78   0.741E+07
        23      0.257E+02    *  0.642E+02      0.00034
        23      0.257E+02       0.278E+02                  0        78   0.344E+07
        24      0.257E+02    *  0.278E+02      0.00167
        24      0.257E+02       0.101E+02                  0        78   0.133E+07
        25      0.255E+02    *  0.100E+02      0.00601
        25      0.255E+02       0.719E+01                  0        78   0.450E+06
        26      0.252E+02    *  0.710E+01      0.01213
        26      0.252E+02       0.732E+01                  0        78   0.230E+06
        27      0.243E+02    *  0.684E+01      0.07263
        27      0.243E+02       0.349E+01                  0        78   0.603E+05
        28      0.518E+02    *  0.274E+01      0.31158
        28      0.518E+02       0.707E+01                  0        78   0.271E+05
        28      0.505E+02    *  0.724E+01      0.02492
        29      0.515E+02    *  0.705E+01      0.00624
        29      0.515E+02       0.838E+03                  0        78   0.570E+07
        29      0.515E+02    *  0.839E+03      0.00010
        29 Not accepted damping factor        0.00001     0        78
        29      0.515E+02       0.169E+01                  0        77   0.322E+05
        30      0.244E+03    *  0.102E+01      1.00000
        30      0.244E+03       0.193E+01                  0        78   0.247E+05
        30      0.225E+03    *  0.197E+01      0.07957
        31      0.237E+03    *  0.190E+01      0.02983
        31      0.237E+03       0.169E+02                  0        78   0.297E+06
        32      0.237E+03    *  0.169E+02      0.00035
        32      0.237E+03       0.246E+02                  0        78   0.464E+06
        33      0.237E+03    *  0.246E+02      0.00011
        33      0.237E+03       0.838E+03                  0        78   0.158E+08
        33      0.237E+03    *  0.882E+03      0.00010
        33 Not accepted damping factor        0.00000     0        78
        33      0.237E+03       0.445E+00                  0        77   0.112E+05
        33      0.104E+03    *  0.491E+01      0.55175
        34      0.234E+03    *  0.438E+00      0.01366
        34      0.234E+03       0.256E+02                  0        78   0.527E+06
        35      0.234E+03    *  0.256E+02      0.00010
        35      0.234E+03       0.211E+03                  0        78   0.433E+07
        35      0.234E+03    *  0.212E+03      0.00010
        35 Not accepted damping factor        0.00000     0        78
        35      0.234E+03       0.446E+00                  0        77   0.791E+04
        35      0.654E+02    *  0.173E+01      0.71328
        36      0.218E+03    *  0.420E+00      0.06692
        36      0.218E+03       0.118E+02                  0        78   0.253E+06
        37      0.218E+03    *  0.118E+02      0.00010
        37      0.218E+03       0.161E+02                  0        78   0.335E+06
        38      0.218E+03    *  0.161E+02      0.00020
        38      0.218E+03       0.152E+05                  0        78   0.316E+09
        38      0.216E+03    *  0.359E+06      0.00010
        38 Not accepted damping factor        0.00000     0        78
        38      0.218E+03       0.455E+00                  0        77   0.493E+04
        39      0.111E+02    *  0.180E+00      1.00000
        39      0.111E+02       0.666E+00                  0        78   0.176E+05
        40      0.900E+01    *  0.621E+00      0.19228
        40      0.900E+01       0.196E+01                  0        78   0.751E+05
        41      0.708E+01    *  0.133E+01      0.21475
        41      0.708E+01       0.865E+00                  0        78   0.261E+05
        42      0.217E+01    *  0.167E+00      0.70551
        42      0.217E+01       0.155E+00                  0        78   0.186E+05
        43      0.341E-01    *  0.559E-01      1.00000
        43      0.341E-01       0.110E+00                  0        78   0.234E+05
        44      0.991E-02    *  0.116E-01      1.00000
        44      0.991E-02       0.157E-01                  0        78   0.333E+05
        45      0.109E-03    *  0.306E-04      1.00000
        45      0.109E-03       0.420E-04                  0        78   0.531E+05
        46      0.109E-08    *  0.278E-08      1.00000
        46      0.109E-08       0.303E-08                  0        78   0.580E+05
        47      0.300E-12    *  0.976E-13      1.00000



 Solution of nonlinear system of equations obtained within  47 iteration steps

 Achieved relative accuracy 0.976E-13

   Subcondition ( 1,  78)  0.580E+05
   Sensitivity ( 1,  78)  0.242E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      47  ***
   ***  Corrector steps   :      10  ***
   ***  Rejected rk-1 st. :       5  ***
   ***  Jacobian eval.    :      47  ***
   ***  Function eval.    :      58  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1666E-12

 problem: Dew Point Temperature                           


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       542 (  0.8 percent)

 Integer Workspace declared as      1000 is used up to        68 (  6.8 percent)



 N =  16

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    16
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.570E-01       0.306E+00                  0        16   0.133E+03
         1      0.564E-01    *  0.303E+00      0.01000
         1      0.564E-01       0.300E+00                  0        16   0.132E+03
         2      0.107E+00    *  0.783E-01      0.77483
         2      0.107E+00       0.424E+00                  0        16   0.269E+03
         3      0.102E+00    *  0.409E+00      0.14366
         3      0.102E+00       0.444E+01                  0        16   0.267E+04
         4      0.102E+00    *  0.443E+01      0.00259
         4      0.102E+00       0.578E+01                  0        16   0.350E+04
         5      0.102E+00    *  0.578E+01      0.00410
         5      0.102E+00       0.762E+01                  0        16   0.481E+04
         6      0.102E+00    *  0.762E+01      0.00144
         6      0.102E+00       0.674E+02                  0        16   0.418E+05
         6      0.102E+00    *  0.674E+02      0.00010
         6 Not accepted damping factor        0.00002     0        16
         6      0.102E+00       0.748E-01                  0        15   0.225E+02
         7      0.931E-02    *  0.185E-01      1.00000
         7      0.931E-02       0.579E-01                  0        16   0.160E+03
         8      0.572E-02    *  0.310E-01      0.44318
         8      0.572E-02       0.265E-01                  0        16   0.154E+03
         9      0.155E-02    *  0.312E-02      1.00000
         9      0.155E-02       0.457E-02                  0        16   0.151E+03
        10      0.167E-04    *  0.921E-04      1.00000
        10      0.167E-04       0.904E-04                  0        16   0.150E+03
        11      0.119E-07    *  0.464E-07      1.00000
        11      0.119E-07       0.462E-07                  0        16   0.149E+03
        12      0.202E-14    *  0.938E-14      1.00000



 Solution of nonlinear system of equations obtained within  12 iteration steps

 Achieved relative accuracy 0.938E-14

   Subcondition ( 1,  16)  0.149E+03
   Sensitivity ( 1,  16)  0.245E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      12  ***
   ***  Corrector steps   :       1  ***
   ***  Rejected rk-1 st. :       1  ***
   ***  Jacobian eval.    :      12  ***
   ***  Function eval.    :      14  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1953E-15

 problem: Partial Oxidation of Methane with Oxygen        


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       209 (  0.3 percent)

 Integer Workspace declared as      1000 is used up to        59 (  5.9 percent)



 N =   7

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     7
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.352E+05       0.833E+00                  0         7   0.399E+01
         1      0.348E+05    *  0.820E+00      0.01000
         1      0.348E+05       0.829E+00                  0         7   0.399E+01
         1      0.188E+05    *  0.952E+00      0.42088
         2      0.300E+05    *  0.615E+00      0.13265
         2      0.300E+05       0.123E+01                  0         7   0.628E+01
         3      0.292E+05    *  0.106E+01      0.02578
         3      0.292E+05       0.354E+01                  0         7   0.184E+02
         4      0.291E+05    *  0.335E+01      0.00373
         4      0.291E+05       0.472E+02                  0         7   0.248E+03
         5      0.291E+05    *  0.464E+02      0.00010
         5      0.291E+05       0.147E+02                  0         7   0.773E+02
         5      0.291E+05    *  0.149E+02      0.00024
         5      0.291E+05    *  0.148E+02      0.00010
         5 Not accepted damping factor        0.00000     0         7
         5      0.291E+05       0.346E+00                  0         6   0.511E+01
         6      0.136E+05    *  0.136E+00      0.86014
         6      0.136E+05       0.460E+00                  0         7   0.471E+01
         7      0.102E+05    *  0.369E+00      0.24647
         7      0.102E+05       0.350E+00                  0         7   0.468E+01
         8      0.346E+03    *  0.998E-01      1.00000
         8      0.346E+03       0.135E+00                  0         7   0.938E+01
         8      0.564E+02    *  0.161E+00      1.00000
         9      0.212E+03    *  0.130E+00      0.42057
         9      0.212E+03       0.166E+01                  0         7   0.285E+02
        10      0.211E+03    *  0.166E+01      0.00064
        10      0.211E+03       0.151E+01                  0         7   0.255E+02
        11      0.211E+03    *  0.148E+01      0.01043
        11      0.211E+03       0.866E+00                  0         7   0.145E+02
        12      0.214E+03    *  0.817E+00      0.04277
        12      0.214E+03       0.696E+00                  0         7   0.115E+02
        13      0.737E+03    *  0.398E+00      0.29521
        13      0.737E+03       0.309E+00                  0         7   0.481E+01
        14      0.125E+04    *  0.934E-01      1.00000
        14      0.125E+04       0.731E+00                  0         7   0.398E+01
        15      0.118E+04    *  0.625E+00      0.10476
        15      0.118E+04       0.593E+00                  0         7   0.434E+01
        16      0.201E+04    *  0.353E+00      0.75016
        16      0.201E+04       0.727E+00                  0         7   0.802E+01
        16      0.164E+04    *  0.848E+00      0.20510
        16      0.191E+04    *  0.755E+00      0.05310
        16      0.198E+04    *  0.734E+00      0.01401
        16      0.200E+04    *  0.729E+00      0.00370
        16      0.201E+04    *  0.727E+00      0.00098
        16      0.201E+04    *  0.727E+00      0.00026
        16      0.201E+04    *  0.727E+00      0.00010
        16 Not accepted damping factor        0.00003     0         7
        16      0.201E+04       0.272E+00                  0         6   0.516E+01
        17      0.884E+04    *  0.106E+00      1.00000
        17      0.884E+04       0.392E+00                  0         7   0.692E+01
        18      0.766E+04    *  0.379E+00      0.13329
        18      0.766E+04       0.337E+00                  0         7   0.703E+01
        18      0.762E+02    *  0.241E+01      1.00000
        19      0.713E+04    *  0.332E+00      0.06998
        19      0.713E+04       0.313E+00                  0         7   0.713E+01
        19      0.228E+04    *  0.370E+00      0.67543
        20      0.529E+04    *  0.298E+00      0.25714
        20      0.529E+04       0.227E+00                  0         7   0.748E+01
        20      0.141E+04    *  0.869E+00      0.72929
        21      0.490E+04    *  0.222E+00      0.07388
        21      0.490E+04       0.207E+00                  0         7   0.757E+01
        21      0.186E+04    *  0.242E+00      0.61781
        22      0.376E+04    *  0.196E+00      0.23190
        22      0.376E+04       0.151E+00                  0         7   0.781E+01
        22      0.119E+04    *  0.255E+00      0.68064
        23      0.314E+04    *  0.144E+00      0.16390
        23      0.314E+04       0.122E+00                  0         7   0.796E+01
        23      0.102E+04    *  0.185E+00      0.67273
        24      0.257E+04    *  0.115E+00      0.18187
        24      0.257E+04       0.956E-01                  0         7   0.810E+01
        24      0.742E+03    *  0.183E+00      0.70969
        25      0.219E+04    *  0.904E-01      0.14998
        25      0.219E+04       0.789E-01                  0         7   0.819E+01
        25      0.553E+03    *  0.212E+00      0.74527
        26      0.194E+04    *  0.754E-01      0.11184
        26      0.194E+04       0.687E-01                  0         7   0.826E+01
        26      0.428E+03    *  0.337E+00      0.77809
        27      0.182E+04    *  0.668E-01      0.06393
        27      0.182E+04       0.637E-01                  0         7   0.829E+01
        27      0.359E+03    *  0.101E+01      0.80086
        28      0.178E+04    *  0.631E-01      0.02048
        28      0.178E+04       0.622E-01                  0         7   0.830E+01
        28      0.335E+03    *  0.321E+01      0.81181
        29      0.177E+04    *  0.620E-01      0.00637
        29      0.177E+04       0.617E-01                  0         7   0.831E+01
        29      0.324E+03    *  0.130E+01      0.81536
        30      0.174E+04    *  0.613E-01      0.01567
        30      0.174E+04       0.607E-01                  0         7   0.831E+01
        30      0.314E+03    *  0.984E+00      0.81849
        31      0.170E+04    *  0.601E-01      0.02048
        31      0.170E+04       0.592E-01                  0         7   0.832E+01
        31      0.298E+03    *  0.663E+00      0.82353
        32      0.165E+04    *  0.584E-01      0.02996
        32      0.165E+04       0.572E-01                  0         7   0.834E+01
        32      0.278E+03    *  0.453E+00      0.83060
        33      0.158E+04    *  0.561E-01      0.04291
        33      0.158E+04       0.545E-01                  0         7   0.836E+01
        33      0.250E+03    *  0.304E+00      0.84083
        34      0.148E+04    *  0.529E-01      0.06201
        34      0.148E+04       0.507E-01                  0         7   0.838E+01
        34      0.212E+03    *  0.201E+00      0.85557
        35      0.135E+04    *  0.486E-01      0.08973
        35      0.135E+04       0.456E-01                  0         7   0.842E+01
        35      0.165E+03    *  0.131E+00      0.87678
        36      0.118E+04    *  0.430E-01      0.12976
        36      0.118E+04       0.391E-01                  0         7   0.847E+01
        36      0.108E+03    *  0.839E-01      0.90692
        37      0.958E+03    *  0.359E-01      0.18549
        37      0.958E+03       0.311E-01                  0         7   0.853E+01
        37      0.488E+02    *  0.535E-01      0.94806
        38      0.713E+03    *  0.277E-01      0.25579
        38      0.713E+03       0.225E-01                  0         7   0.861E+01
        38      0.247E+01    *  0.344E-01      0.99898
        39      0.480E+03    *  0.195E-01      0.32623
        39      0.480E+03       0.147E-01                  0         7   0.869E+01
        39      0.609E+01    *  0.409E-01      1.00000
        40      0.394E+03    *  0.134E-01      0.17979
        40      0.394E+03       0.119E-01                  0         7   0.872E+01
        40      0.926E+01    *  0.433E-01      1.00000
        41      0.340E+03    *  0.111E-01      0.13767
        41      0.340E+03       0.102E-01                  0         7   0.874E+01
        41      0.126E+02    *  0.448E-01      1.00000
        42      0.301E+03    *  0.959E-02      0.11386
        42      0.301E+03       0.900E-02                  0         7   0.875E+01
        42      0.161E+02    *  0.459E-01      1.00000
        43      0.271E+03    *  0.852E-02      0.09798
        43      0.271E+03       0.808E-02                  0         7   0.876E+01
        43      0.199E+02    *  0.468E-01      1.00000
        44      0.248E+03    *  0.770E-02      0.08643
        44      0.248E+03       0.736E-02                  0         7   0.877E+01
        44      0.239E+02    *  0.475E-01      1.00000
        45      0.229E+03    *  0.705E-02      0.07758
        45      0.229E+03       0.677E-02                  0         7   0.878E+01
        45      0.281E+02    *  0.480E-01      1.00000
        46      0.213E+03    *  0.651E-02      0.07053
        46      0.213E+03       0.628E-02                  0         7   0.879E+01
        46      0.326E+02    *  0.485E-01      1.00000
        47      0.199E+03    *  0.606E-02      0.06477
        47      0.199E+03       0.586E-02                  0         7   0.879E+01
        47      0.373E+02    *  0.489E-01      1.00000
        48      0.187E+03    *  0.567E-02      0.05995
        48      0.187E+03       0.550E-02                  0         7   0.880E+01
        48      0.422E+02    *  0.493E-01      1.00000
        49      0.176E+03    *  0.533E-02      0.05586
        49      0.176E+03       0.519E-02                  0         7   0.880E+01
        49      0.474E+02    *  0.496E-01      1.00000
        50      0.167E+03    *  0.504E-02      0.05232
        50      0.167E+03       0.491E-02                  0         7   0.880E+01
        50      0.528E+02    *  0.498E-01      1.00000
        51      0.159E+03    *  0.477E-02      0.04924
        51      0.159E+03       0.466E-02                  0         7   0.881E+01
        51      0.584E+02    *  0.501E-01      1.00000
        52      0.152E+03    *  0.454E-02      0.04653
        52      0.152E+03       0.444E-02                  0         7   0.881E+01
        52      0.642E+02    *  0.503E-01      1.00000
        53      0.145E+03    *  0.433E-02      0.04412
        53      0.145E+03       0.424E-02                  0         7   0.881E+01
        53      0.703E+02    *  0.505E-01      1.00000
        54      0.139E+03    *  0.414E-02      0.04196
        54      0.139E+03       0.406E-02                  0         7   0.882E+01
        54      0.766E+02    *  0.507E-01      1.00000
        55      0.133E+03    *  0.397E-02      0.04002
        55      0.133E+03       0.389E-02                  0         7   0.882E+01
        55      0.831E+02    *  0.509E-01      1.00000
        56      0.128E+03    *  0.381E-02      0.03826
        56      0.128E+03       0.374E-02                  0         7   0.882E+01
        56      0.899E+02    *  0.510E-01      1.00000
        57      0.124E+03    *  0.366E-02      0.03665
        57      0.124E+03       0.360E-02                  0         7   0.882E+01
        57      0.968E+02    *  0.512E-01      1.00000
        58      0.119E+03    *  0.353E-02      0.03518
        58      0.119E+03       0.347E-02                  0         7   0.883E+01
        58      0.104E+03    *  0.513E-01      1.00000
        59      0.115E+03    *  0.341E-02      0.03383
        59      0.115E+03       0.335E-02                  0         7   0.883E+01
        59      0.111E+03    *  0.514E-01      1.00000
        60      0.111E+03    *  0.329E-02      0.03259
        60      0.111E+03       0.324E-02                  0         7   0.883E+01
        60      0.119E+03    *  0.515E-01      1.00000
        61      0.108E+03    *  0.318E-02      0.03144
        61      0.108E+03       0.314E-02                  0         7   0.883E+01
        61      0.127E+03    *  0.517E-01      1.00000
        62      0.105E+03    *  0.308E-02      0.03037
        62      0.105E+03       0.304E-02                  0         7   0.883E+01
        62      0.135E+03    *  0.518E-01      1.00000
        63      0.102E+03    *  0.299E-02      0.02938
        63      0.102E+03       0.295E-02                  0         7   0.883E+01
        63      0.143E+03    *  0.518E-01      1.00000
        64      0.987E+02    *  0.290E-02      0.02845
        64      0.987E+02       0.286E-02                  0         7   0.883E+01
        64      0.152E+03    *  0.519E-01      1.00000
        65      0.959E+02    *  0.282E-02      0.02758
        65      0.959E+02       0.278E-02                  0         7   0.884E+01
        65      0.160E+03    *  0.520E-01      1.00000
        66      0.934E+02    *  0.274E-02      0.02676
        66      0.934E+02       0.271E-02                  0         7   0.884E+01
        66      0.169E+03    *  0.521E-01      1.00000
        67      0.910E+02    *  0.267E-02      0.02600
        67      0.910E+02       0.264E-02                  0         7   0.884E+01
        67      0.179E+03    *  0.522E-01      1.00000
        68      0.887E+02    *  0.260E-02      0.02527
        68      0.887E+02       0.257E-02                  0         7   0.884E+01
        68      0.188E+03    *  0.523E-01      1.00000
        69      0.865E+02    *  0.253E-02      0.02459
        69      0.865E+02       0.251E-02                  0         7   0.884E+01
        69      0.197E+03    *  0.523E-01      1.00000
        70      0.844E+02    *  0.247E-02      0.02395
        70      0.844E+02       0.244E-02                  0         7   0.884E+01
        70      0.207E+03    *  0.524E-01      1.00000
        71      0.825E+02    *  0.241E-02      0.02334
        71      0.825E+02       0.239E-02                  0         7   0.884E+01
        71      0.217E+03    *  0.524E-01      1.00000
        72      0.806E+02    *  0.236E-02      0.02276
        72      0.806E+02       0.233E-02                  0         7   0.884E+01
        72      0.227E+03    *  0.525E-01      1.00000
        73      0.788E+02    *  0.230E-02      0.02221
        73      0.788E+02       0.228E-02                  0         7   0.884E+01
        73      0.238E+03    *  0.526E-01      1.00000
        74      0.771E+02    *  0.225E-02      0.02169
        74      0.771E+02       0.223E-02                  0         7   0.884E+01
        74      0.249E+03    *  0.526E-01      1.00000
        75      0.755E+02    *  0.220E-02      0.02119
        75      0.755E+02       0.218E-02                  0         7   0.884E+01
        75      0.259E+03    *  0.527E-01      1.00000
        76      0.739E+02    *  0.216E-02      0.02071
        76      0.739E+02       0.214E-02                  0         7   0.885E+01
        76      0.270E+03    *  0.527E-01      1.00000
        77      0.724E+02    *  0.211E-02      0.02026
        77      0.724E+02       0.209E-02                  0         7   0.885E+01
        77      0.282E+03    *  0.528E-01      1.00000
        78      0.710E+02    *  0.207E-02      0.01983
        78      0.710E+02       0.205E-02                  0         7   0.885E+01
        78      0.293E+03    *  0.528E-01      1.00000
        79      0.696E+02    *  0.203E-02      0.01941
        79      0.696E+02       0.201E-02                  0         7   0.885E+01
        79      0.305E+03    *  0.528E-01      1.00000
        80      0.683E+02    *  0.199E-02      0.01902
        80      0.683E+02       0.197E-02                  0         7   0.885E+01
        80      0.317E+03    *  0.529E-01      1.00000
        81      0.670E+02    *  0.195E-02      0.01864
        81      0.670E+02       0.193E-02                  0         7   0.885E+01
        81      0.329E+03    *  0.529E-01      1.00000
        82      0.658E+02    *  0.191E-02      0.01827
        82      0.658E+02       0.190E-02                  0         7   0.885E+01
        82      0.341E+03    *  0.530E-01      1.00000
        83      0.647E+02    *  0.188E-02      0.01792
        83      0.647E+02       0.186E-02                  0         7   0.885E+01
        83      0.354E+03    *  0.530E-01      1.00000
        84      0.635E+02    *  0.185E-02      0.01758
        84      0.635E+02       0.183E-02                  0         7   0.885E+01
        84      0.367E+03    *  0.530E-01      1.00000
        85      0.624E+02    *  0.181E-02      0.01726
        85      0.624E+02       0.180E-02                  0         7   0.885E+01
        85      0.380E+03    *  0.531E-01      1.00000
        86      0.614E+02    *  0.178E-02      0.01695
        86      0.614E+02       0.177E-02                  0         7   0.885E+01
        86      0.393E+03    *  0.531E-01      1.00000
        87      0.604E+02    *  0.175E-02      0.01665
        87      0.604E+02       0.174E-02                  0         7   0.885E+01
        87      0.406E+03    *  0.531E-01      1.00000
        88      0.594E+02    *  0.172E-02      0.01636
        88      0.594E+02       0.171E-02                  0         7   0.885E+01
        88      0.420E+03    *  0.532E-01      1.00000
        89      0.585E+02    *  0.169E-02      0.01608
        89      0.585E+02       0.168E-02                  0         7   0.885E+01
        89      0.433E+03    *  0.532E-01      1.00000
        90      0.575E+02    *  0.167E-02      0.01581
        90      0.575E+02       0.166E-02                  0         7   0.885E+01
        90      0.447E+03    *  0.532E-01      1.00000
        91      0.567E+02    *  0.164E-02      0.01555
        91      0.567E+02       0.163E-02                  0         7   0.885E+01
        91      0.462E+03    *  0.533E-01      1.00000
        92      0.558E+02    *  0.162E-02      0.01530
        92      0.558E+02       0.160E-02                  0         7   0.885E+01
        92      0.476E+03    *  0.533E-01      1.00000
        93      0.550E+02    *  0.159E-02      0.01505
        93      0.550E+02       0.158E-02                  0         7   0.885E+01
        93      0.491E+03    *  0.533E-01      1.00000
        94      0.542E+02    *  0.157E-02      0.01482
        94      0.542E+02       0.156E-02                  0         7   0.885E+01
        94      0.506E+03    *  0.533E-01      1.00000
        95      0.534E+02    *  0.154E-02      0.01459
        95      0.534E+02       0.153E-02                  0         7   0.886E+01
        95      0.521E+03    *  0.534E-01      1.00000
        96      0.526E+02    *  0.152E-02      0.01437
        96      0.526E+02       0.151E-02                  0         7   0.886E+01
        96      0.536E+03    *  0.534E-01      1.00000
        97      0.519E+02    *  0.150E-02      0.01415
        97      0.519E+02       0.149E-02                  0         7   0.886E+01
        97      0.551E+03    *  0.534E-01      1.00000
        98      0.512E+02    *  0.148E-02      0.01394
        98      0.512E+02       0.147E-02                  0         7   0.886E+01
        98      0.567E+03    *  0.534E-01      1.00000
        99      0.505E+02    *  0.146E-02      0.01374
        99      0.505E+02       0.145E-02                  0         7   0.886E+01
        99      0.583E+03    *  0.535E-01      1.00000
       100      0.498E+02    *  0.144E-02      0.01354

 Iteration terminates after NITMAX =100  Iteration steps

   Subcondition ( 1,   7)  0.886E+01
   Sensitivity ( 1,   7)  0.224E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :     101  ***
   ***  Corrector steps   :      93  ***
   ***  Rejected rk-1 st. :       2  ***
   ***  Jacobian eval.    :     100  ***
   ***  Function eval.    :     194  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.4984E+02

 problem: Swirling Flow between Disks                     


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to     11038 ( 15.8 percent)

 Integer Workspace declared as      1000 is used up to       150 ( 15.0 percent)



 N =  98

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    98
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.619E+00       0.920E+04                  0        98   0.750E+07
         1      0.613E+00    *  0.911E+04      0.01000
         1      0.613E+00       0.105E+03                  0        98   0.233E+07
         2      0.307E+00    *  0.594E+02      0.48859
         2      0.307E+00       0.421E+01                  0        98   0.166E+07
         3      0.244E+00    *  0.343E+01      0.20470
         3      0.244E+00       0.154E+01                  0        98   0.131E+07
         4      0.153E+00    *  0.113E+01      0.37302
         4      0.153E+00       0.220E+01                  0        98   0.946E+06
         5      0.117E+00    *  0.178E+01      0.22977
         5      0.117E+00       0.188E+01                  0        98   0.289E+06
         6      0.257E-01    *  0.766E+00      0.90294
         6      0.257E-01       0.524E+00                  0        98   0.120E+06
         7      0.590E-02    *  0.181E+00      1.00000
         7      0.590E-02       0.242E+00                  0        98   0.522E+05
         8      0.324E-03    *  0.357E-01      1.00000
         8      0.324E-03       0.357E-01                  0        98   0.796E+05
         9      0.163E-05    *  0.112E-03      1.00000
         9      0.163E-05       0.142E-03                  0        98   0.830E+05
        10      0.214E-10    *  0.585E-08      1.00000
        10      0.214E-10       0.561E-08                  0        98   0.830E+05
        11      0.253E-13    *  0.844E-14      1.00000



 Solution of nonlinear system of equations obtained within  11 iteration steps

 Achieved relative accuracy 0.844E-14

   Subcondition ( 1,  98)  0.830E+05
   Sensitivity ( 1,  98)  0.316E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      11  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      11  ***
   ***  Function eval.    :      12  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.2872E-13

 problem: Solid Fuel Ignition                             


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to     11462 ( 16.4 percent)

 Integer Workspace declared as      1000 is used up to       152 ( 15.2 percent)



 N = 100

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =   100
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.116E+00       0.798E-01                  0       100   0.794E+01
         1      0.115E+00    *  0.790E-01      0.01000
         1      0.115E+00       0.790E-01                  0       100   0.847E+01
         2      0.262E-03    *  0.803E-03      1.00000
         2      0.262E-03       0.788E-03                  0       100   0.785E+01
         3      0.213E-07    *  0.213E-06      1.00000
         3      0.213E-07       0.213E-06                  0       100   0.785E+01
         4      0.207E-14    *  0.198E-13      1.00000



 Solution of nonlinear system of equations obtained within   4 iteration steps

 Achieved relative accuracy 0.198E-13

   Subcondition ( 1, 100)  0.785E+01
   Sensitivity ( 1, 100)  0.112E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       4  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       4  ***
   ***  Function eval.    :       5  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.8720E-16

 problem: Sulphur Dioxide to Sulphur Trioxide             


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       113 (  0.2 percent)

 Integer Workspace declared as      1000 is used up to        55 (  5.5 percent)



 N =   3

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     3
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.414E-01       0.194E+01                  0         3   0.344E+02
         1      0.458E-01    *  0.192E+01      0.01000
         1      0.458E-01       0.180E+01                  0         3   0.334E+02
         2      0.307E+01    *  0.134E+01      0.13949
         2      0.307E+01       0.448E+00                  0         3   0.228E+02
         2      0.116E+02    *  0.131E+02      0.98023
         3      0.302E+01    *  0.440E+00      0.01639
         3      0.302E+01       0.583E+00                  0         3   0.316E+02
         4      0.344E+01    *  0.226E+00      0.35043
         4      0.344E+01       0.581E-01                  0         3   0.243E+02
         5      0.964E+00    *  0.388E-01      1.00000
         5      0.964E+00       0.369E-01                  0         3   0.263E+02
         6      0.276E-03    *  0.125E-02      1.00000
         6      0.276E-03       0.143E-02                  0         3   0.283E+02
         7      0.569E-04    *  0.502E-05      1.00000
         7      0.569E-04       0.489E-05                  0         3   0.274E+02
         8      0.293E-09    *  0.412E-10      1.00000



 Solution of nonlinear system of equations obtained within   8 iteration steps

 Achieved relative accuracy 0.412E-10

   Subcondition ( 1,   3)  0.274E+02
   Sensitivity ( 1,   3)  0.141E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       8  ***
   ***  Corrector steps   :       1  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       8  ***
   ***  Function eval.    :      10  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.8402E-11

 problem: Compressibility factor from the RK equation     


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to        77 (  0.1 percent)

 Integer Workspace declared as      1000 is used up to        53 (  5.3 percent)



 N =   1

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     1
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.178E+00       0.226E+01                  0         1   0.100E+01
         1      0.176E+00    *  0.223E+01      0.01000
         1      0.176E+00       0.146E+01                  0         1   0.100E+01
         2      0.165E+00    *  0.137E+01      0.04485
         2      0.165E+00       0.738E+00                  0         1   0.100E+01
         3      0.122E+00    *  0.546E+00      0.19301
         3      0.122E+00       0.341E+00                  0         1   0.100E+01
         4      0.968E-02    *  0.271E-01      1.00000
         4      0.968E-02       0.525E-01                  0         1   0.100E+01
         5      0.225E-02    *  0.122E-01      1.00000
         5      0.225E-02       0.235E-01                  0         1   0.100E+01
         6      0.524E-03    *  0.547E-02      1.00000
         6      0.524E-03       0.103E-01                  0         1   0.100E+01
         7      0.107E-03    *  0.211E-02      1.00000
         7      0.107E-03       0.359E-02                  0         1   0.100E+01
         8      0.132E-04    *  0.443E-03      1.00000
         8      0.132E-04       0.589E-03                  0         1   0.100E+01
         9      0.359E-06    *  0.160E-04      1.00000
         9      0.359E-06       0.169E-04                  0         1   0.100E+01
        10      0.296E-09    *  0.139E-07      1.00000
        10      0.296E-09       0.140E-07                  0         1   0.100E+01
        11      0.202E-15    *  0.952E-14      1.00000



 Solution of nonlinear system of equations obtained within  11 iteration steps

 Achieved relative accuracy 0.952E-14

   Subcondition ( 1,   1)  0.100E+01
   Sensitivity ( 1,   1)  0.100E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      11  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      11  ***
   ***  Function eval.    :      12  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.0000E+00

 problem: Chem. Reactor Equil. Conversion, Shacham 1983   


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to        94 (  0.1 percent)

 Integer Workspace declared as      1000 is used up to        54 (  5.4 percent)



 N =   2

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     2
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.440E+04       0.679E-01                  0         2   0.476E+01
         1      0.435E+04    *  0.672E-01      0.01000
         1      0.435E+04       0.672E-01                  0         2   0.477E+01
         2      0.384E+01    *  0.407E-02      1.00000
         2      0.384E+01       0.143E-01                  0         2   0.503E+01
         3      0.394E+00    *  0.492E-02      1.00000
         3      0.394E+00       0.124E-01                  0         2   0.475E+01
         4      0.216E+00    *  0.551E-02      0.75448
         4      0.216E+00       0.106E-01                  0         2   0.449E+01
         5      0.140E+00    *  0.335E-02      0.97235
         5      0.140E+00       0.722E-02                  0         2   0.404E+01
         6      0.645E-01    *  0.176E-02      1.00000
         6      0.645E-01       0.321E-02                  0         2   0.349E+01
         7      0.126E-01    *  0.421E-03      1.00000
         7      0.126E-01       0.564E-03                  0         2   0.313E+01
         8      0.389E-03    *  0.147E-04      1.00000
         8      0.389E-03       0.155E-04                  0         2   0.305E+01
         9      0.295E-06    *  0.115E-07      1.00000
         9      0.295E-06       0.115E-07                  0         2   0.305E+01
        10      0.103E-10    *  0.616E-14      1.00000



 Solution of nonlinear system of equations obtained within  10 iteration steps

 Achieved relative accuracy 0.616E-14

   Subcondition ( 1,   2)  0.305E+01
   Sensitivity ( 1,   2)  0.141E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      10  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      10  ***
   ***  Function eval.    :      11  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1029E-10

 problem: Chemical Reactor Steady State, Shacham 1983     


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to        94 (  0.1 percent)

 Integer Workspace declared as      1000 is used up to        54 (  5.4 percent)



 N =   2

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     2
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.841E+02       0.209E-01                  0         2   0.294E+01
         1      0.833E+02    *  0.207E-01      0.01000
         1      0.833E+02       0.205E-01                  0         2   0.296E+01
         2      0.376E+02    *  0.834E-02      0.50365
         2      0.376E+02       0.118E-01                  0         2   0.402E+01
         3      0.359E+01    *  0.944E-03      0.90776
         3      0.359E+01       0.950E-03                  0         2   0.307E+01
         4      0.367E-02    *  0.131E-05      1.00000
         4      0.367E-02       0.125E-05                  0         2   0.302E+01
         5      0.217E-07    *  0.789E-11      1.00000



 Solution of nonlinear system of equations obtained within   5 iteration steps

 Achieved relative accuracy 0.789E-11

   Subcondition ( 1,   2)  0.302E+01
   Sensitivity ( 1,   2)  0.141E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       5  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       5  ***
   ***  Function eval.    :       6  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.6343E-12

 problem: Cutlips steady state for reaction rate equations


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       182 (  0.3 percent)

 Integer Workspace declared as      1000 is used up to        58 (  5.8 percent)



 N =   6

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     6
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.150E+02       0.366E+02                  0         6   0.139E+05
         1      0.149E+02    *  0.362E+02      0.01000
         1      0.149E+02       0.362E+02                  0         6   0.138E+05
         2      0.263E-01    *  0.307E+01      1.00000
         2      0.263E-01       0.865E+00                  0         6   0.147E+04
         3      0.524E-03    *  0.303E-02      1.00000
         3      0.524E-03       0.209E-02                  0         6   0.113E+02
         4      0.193E-07    *  0.856E-07      1.00000
         4      0.193E-07       0.791E-07                  0         6   0.599E+01
         5      0.268E-16    *  0.791E-16      1.00000



 Solution of nonlinear system of equations obtained within   5 iteration steps

 Achieved relative accuracy 0.791E-16

   Subcondition ( 1,   6)  0.599E+01
   Sensitivity ( 1,   6)  0.201E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       5  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       5  ***
   ***  Function eval.    :       6  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.2647E-16

 problem: Hodgkin Huxley 1                                


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       134 (  0.2 percent)

 Integer Workspace declared as      1000 is used up to        56 (  5.6 percent)



 N =   4

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by numerical differentiation (without feedback strategy)
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     4
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.225E+01       0.109E+02                  0         4   0.206E+03
         1      0.222E+01    *  0.108E+02      0.01000
         1      0.222E+01       0.106E+02                  0         4   0.202E+03
         1      0.683E+02    *  0.153E+03      0.82089
         2      0.217E+01    *  0.104E+02      0.02310
         2      0.217E+01       0.100E+02                  0         4   0.194E+03
         2      0.270E+02    *  0.539E+02      0.65780
         3      0.207E+01    *  0.961E+01      0.03773
         3      0.207E+01       0.744E+01                  0         4   0.151E+03
         3      0.955E+01    *  0.124E+02      0.52135
         4      0.190E+01    *  0.692E+01      0.06332
         4      0.190E+01       0.318E+01                  0         4   0.762E+02
         5      0.406E+01    *  0.183E+01      0.46352
         5      0.406E+01       0.410E+00                  0         4   0.196E+02
         6      0.435E+01    *  0.150E+00      1.00000
         6      0.435E+01       0.242E+00                  0         4   0.165E+02
         7      0.995E+00    *  0.548E-01      1.00000
         7      0.995E+00       0.212E+00                  0         4   0.438E+02
         8      0.111E+00    *  0.227E-01      1.00000
         8      0.111E+00       0.838E-01                  0         4   0.135E+03
         9      0.171E-02    *  0.134E-02      1.00000
         9      0.171E-02       0.138E-02                  0         4   0.136E+03
        10      0.425E-06    *  0.338E-06      1.00000
        10      0.425E-06       0.338E-06                  0         4   0.136E+03
        11      0.503E-14    *  0.943E-13      1.00000



 Solution of nonlinear system of equations obtained within  11 iteration steps

 Achieved relative accuracy 0.943E-13

   Subcondition ( 1,   4)  0.136E+03
   Sensitivity ( 1,   4)  0.141E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      11  ***
   ***  Corrector steps   :       3  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      11  ***
   ***  Function eval.    :      15  ***
   ***  ...  for Jacobian :      44  ***
   *************************************

 norm of residuum =   0.4596E-15

 problem: Inf Reflux                                      


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to        77 (  0.1 percent)

 Integer Workspace declared as      1000 is used up to        53 (  5.3 percent)



 N =   1

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     1
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.190E-01       0.340E+02                  0         1   0.100E+01
         1      0.138E-01    *  0.246E+02      0.01000
         1      0.138E-01       0.571E-01                  0         1   0.100E+01
         2      0.212E-02    *  0.880E-02      1.00000
         2      0.212E-02       0.124E-01                  0         1   0.100E+01
         3      0.758E-04    *  0.443E-03      1.00000
         3      0.758E-04       0.476E-03                  0         1   0.100E+01
         4      0.107E-06    *  0.670E-06      1.00000
         4      0.107E-06       0.672E-06                  0         1   0.100E+01
         5      0.212E-12    *  0.133E-11      1.00000



 Solution of nonlinear system of equations obtained within   5 iteration steps

 Achieved relative accuracy 0.133E-11

   Subcondition ( 1,   1)  0.100E+01
   Sensitivity ( 1,   1)  0.100E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       5  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       5  ***
   ***  Function eval.    :       6  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.5551E-16

 problem: Order 10 to 11                                  


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to        77 (  0.1 percent)

 Integer Workspace declared as      1000 is used up to        53 (  5.3 percent)



 N =   1

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     1
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.231E+11       0.658E-02                  0         1   0.100E+01
         1      0.228E+11    *  0.652E-02      0.01000
         1      0.228E+11       0.650E-02                  0         1   0.100E+01
         2      0.342E+10    *  0.973E-03      1.00000
         2      0.342E+10       0.743E-03                  0         1   0.100E+01
         3      0.529E+08    *  0.115E-04      1.00000
         3      0.529E+08       0.119E-04                  0         1   0.100E+01
         4      0.132E+05    *  0.297E-08      1.00000
         4      0.132E+05       0.297E-08                  0         1   0.100E+01
         5      0.166E-02    *  0.374E-15      1.00000



 Solution of nonlinear system of equations obtained within   5 iteration steps

 Achieved relative accuracy 0.374E-15

   Subcondition ( 1,   1)  0.100E+01
   Sensitivity ( 1,   1)  0.100E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       5  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       5  ***
   ***  Function eval.    :       6  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.7935E-03

 problem: Chemical equil. from a Partial Methane Oxidation


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       209 (  0.3 percent)

 Integer Workspace declared as      1000 is used up to        59 (  5.9 percent)



 N =   7

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     7
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.767E+02       0.145E-01                  0         7   0.897E+01
         1      0.759E+02    *  0.144E-01      0.01000
         1      0.759E+02       0.232E-01                  0         7   0.233E+02
         2      0.756E+02    *  0.231E-01      0.00362
         2      0.756E+02       0.180E-01                  0         7   0.272E+02
         3      0.753E+02    *  0.178E-01      0.00448
         3      0.753E+02       0.211E-01                  0         7   0.145E+02
         3      0.751E+02    *  0.212E-01      0.00199
         3      0.753E+02    *  0.211E-01      0.00034
         3      0.753E+02    *  0.211E-01      0.00010
         3 Not accepted damping factor        0.00002     0         7
         3      0.753E+02       0.124E-01                  0         6   0.730E+01
         4      0.752E+02    *  0.123E-01      0.00756
         4      0.752E+02       0.624E-01                  0         7   0.448E+02
         4      0.751E+02    *  0.625E-01      0.00031
         4      0.751E+02    *  0.624E-01      0.00010
         4 Not accepted damping factor        0.00001     0         7
         4      0.752E+02       0.704E-02                  0         6   0.771E+01
         5      0.796E+02    *  0.686E-02      0.02681
         5      0.796E+02       0.612E+00                  0         7   0.939E+03
         5      0.796E+02    *  0.619E+00      0.00010
         5 Not accepted damping factor        0.00000     0         7
         5      0.796E+02       0.748E-02                  0         6   0.463E+01
         6      0.780E+02    *  0.739E-02      0.02568
         6      0.780E+02       0.825E+00                  0         7   0.764E+03
         6      0.780E+02    *  0.832E+00      0.00010
         6 Not accepted damping factor        0.00000     0         7
         6      0.780E+02       0.873E-02                  0         6   0.521E+01
         7      0.798E+02    *  0.862E-02      0.01910
         7      0.798E+02       0.138E-01                  0         7   0.815E+01
         8      0.793E+02    *  0.137E-01      0.00641
         8      0.793E+02       0.610E-01                  0         7   0.380E+02
         8      0.793E+02    *  0.611E-01      0.00029
         8      0.793E+02    *  0.611E-01      0.00010
         8 Not accepted damping factor        0.00001     0         7
         8      0.793E+02       0.953E-02                  0         6   0.754E+01
         9      0.805E+02    *  0.948E-02      0.00708
         9      0.805E+02       0.360E-01                  0         7   0.556E+02
        10      0.805E+02    *  0.358E-01      0.00048
        10      0.805E+02       0.110E+00                  0         7   0.702E+02
        10      0.805E+02    *  0.110E+00      0.00010
        10 Not accepted damping factor        0.00000     0         7
        10      0.805E+02       0.133E-01                  0         6   0.154E+02
        11      0.805E+02    *  0.133E-01      0.00133
        11      0.805E+02       0.157E-01                  0         7   0.100E+02
        12      0.804E+02    *  0.157E-01      0.00136
        12      0.804E+02       0.159E+00                  0         7   0.647E+02
        12      0.804E+02    *  0.159E+00      0.00010
        12 Not accepted damping factor        0.00001     0         7
        12      0.804E+02       0.764E-02                  0         6   0.222E+02
        13      0.812E+02    *  0.758E-02      0.00474
        13      0.812E+02       0.213E-01                  0         7   0.139E+02
        14      0.811E+02    *  0.212E-01      0.00074
        14      0.811E+02       0.177E+00                  0         7   0.714E+02
        15      0.811E+02    *  0.177E+00      0.00010
        15      0.811E+02       0.174E-01                  0         7   0.130E+02
        16      0.810E+02    *  0.174E-01      0.00102
        16      0.810E+02       0.321E+00                  0         7   0.140E+03
        16      0.810E+02    *  0.321E+00      0.00010
        16 Not accepted damping factor        0.00000     0         7
        16      0.810E+02       0.102E-01                  0         6   0.141E+02
        17      0.814E+02    *  0.102E-01      0.00182
        17      0.814E+02       0.327E+00                  0         7   0.311E+03
        17      0.814E+02    *  0.328E+00      0.00010
        17 Not accepted damping factor        0.00000     0         7
        17      0.814E+02       0.104E-01                  0         6   0.655E+01
        18      0.813E+02    *  0.104E-01      0.00349
        18      0.813E+02       0.131E+00                  0         7   0.609E+02
        18      0.813E+02    *  0.131E+00      0.00010
        18 Not accepted damping factor        0.00001     0         7
        18      0.813E+02       0.458E-01                  0         6   0.361E+02
        19      0.814E+02    *  0.458E-01      0.00018
        19      0.814E+02       0.306E-01                  0         7   0.231E+02
        20      0.813E+02    *  0.306E-01      0.00019
        20      0.813E+02       0.476E-01                  0         7   0.395E+02
        20      0.813E+02    *  0.476E-01      0.00010
        20 Not accepted damping factor        0.00003     0         7
        20      0.813E+02       0.141E-01                  0         6   0.102E+02
        21      0.813E+02    *  0.141E-01      0.00046
        21      0.813E+02       0.140E+00                  0         7   0.131E+03
        22      0.813E+02    *  0.139E+00      0.00010
        22      0.813E+02       0.551E-01                  0         7   0.681E+02
        22      0.813E+02    *  0.552E-01      0.00020
        22      0.813E+02    *  0.551E-01      0.00010
        22 Not accepted damping factor        0.00001     0         7
        22      0.813E+02       0.157E-01                  0         6   0.175E+02
        23      0.813E+02    *  0.157E-01      0.00137
        23      0.813E+02       0.194E-01                  0         7   0.158E+02
        24      0.813E+02    *  0.194E-01      0.00110
        24      0.813E+02       0.382E+00                  0         7   0.168E+03
        24      0.813E+02    *  0.382E+00      0.00010
        24 Not accepted damping factor        0.00000     0         7
        24      0.813E+02       0.966E-02                  0         6   0.123E+02
        25      0.818E+02    *  0.964E-02      0.00227
        25      0.818E+02       0.667E-01                  0         7   0.617E+02
        25      0.818E+02    *  0.667E-01      0.00010
        25 Not accepted damping factor        0.00001     0         7
        25      0.818E+02       0.959E-02                  0         6   0.641E+01
        26      0.816E+02    *  0.956E-02      0.00431
        26      0.816E+02       0.307E+00                  0         7   0.109E+03
        26      0.816E+02    *  0.307E+00      0.00010
        26 Not accepted damping factor        0.00000     0         7
        26      0.816E+02       0.394E-01                  0         6   0.414E+02
        26      0.816E+02    *  0.394E-01      0.00026
        26 Not accepted damping factor        0.00003     0         6
        26      0.816E+02       0.620E-02                  0         5   0.766E+01
        27      0.846E+02    *  0.615E-02      0.01441
        27      0.846E+02       0.197E-01                  0         7   0.223E+02
        28      0.845E+02    *  0.196E-01      0.00142
        28      0.845E+02       0.237E-01                  0         7   0.164E+02
        28      0.844E+02    *  0.238E-01      0.00100
        28      0.844E+02    *  0.237E-01      0.00015
        28      0.845E+02    *  0.237E-01      0.00010
        28 Not accepted damping factor        0.00001     0         7
        28      0.845E+02       0.131E-01                  0         6   0.111E+02
        29      0.844E+02    *  0.130E-01      0.00419
        29      0.844E+02       0.134E-01                  0         7   0.797E+01
        30      0.840E+02    *  0.134E-01      0.00400
        30      0.840E+02       0.307E-01                  0         7   0.207E+02
        30      0.839E+02    *  0.308E-01      0.00103
        30      0.840E+02    *  0.307E-01      0.00019
        30      0.840E+02    *  0.307E-01      0.00010
        30 Not accepted damping factor        0.00002     0         7
        30      0.840E+02       0.110E-01                  0         6   0.694E+01
        31      0.846E+02    *  0.109E-01      0.00716
        31      0.846E+02       0.115E-01                  0         7   0.185E+02
        32      0.841E+02    *  0.114E-01      0.00662
        32      0.841E+02       0.276E-01                  0         7   0.184E+02
        32      0.840E+02    *  0.277E-01      0.00136
        32      0.841E+02    *  0.276E-01      0.00019
        32      0.841E+02    *  0.276E-01      0.00010
        32 Not accepted damping factor        0.00001     0         7
        32      0.841E+02       0.117E-01                  0         6   0.849E+01
        33      0.843E+02    *  0.115E-01      0.01328
        33      0.843E+02       0.121E-01                  0         7   0.827E+01
        34      0.830E+02    *  0.120E-01      0.01567
        34      0.830E+02       0.155E-01                  0         7   0.204E+02
        35      0.820E+02    *  0.154E-01      0.01194
        35      0.820E+02       0.322E+00                  0         7   0.158E+03
        35      0.820E+02    *  0.322E+00      0.00010
        35 Not accepted damping factor        0.00001     0         7
        35      0.820E+02       0.826E-02                  0         6   0.604E+01
        36      0.866E+02    *  0.804E-02      0.03447
        36      0.866E+02       0.147E-01                  0         7   0.109E+02
        37      0.855E+02    *  0.145E-01      0.01241
        37      0.855E+02       0.238E-01                  0         7   0.123E+02
        37      0.852E+02    *  0.239E-01      0.00399
        37      0.854E+02    *  0.238E-01      0.00090
        37      0.855E+02    *  0.238E-01      0.00020
        37      0.855E+02    *  0.238E-01      0.00010
        37 Not accepted damping factor        0.00002     0         7
        37      0.855E+02       0.782E-02                  0         6   0.867E+01
        38      0.907E+02    *  0.757E-02      0.03701
        38      0.907E+02       0.687E-01                  0         7   0.432E+02
        39      0.907E+02    *  0.687E-01      0.00043
        39      0.907E+02       0.432E-01                  0         7   0.273E+02
        39      0.907E+02    *  0.433E-01      0.00061
        39      0.907E+02    *  0.432E-01      0.00010
        39 Not accepted damping factor        0.00001     0         7
        39      0.907E+02       0.100E-01                  0         6   0.101E+02
        40      0.905E+02    *  0.100E-01      0.00273
        40      0.905E+02       0.179E+00                  0         7   0.574E+03
        40      0.905E+02    *  0.180E+00      0.00010
        40 Not accepted damping factor        0.00000     0         7
        40      0.905E+02       0.307E-01                  0         6   0.269E+02
        40      0.905E+02    *  0.307E-01      0.00028
        40 Not accepted damping factor        0.00004     0         6
        40      0.905E+02       0.387E-02                  0         5   0.326E+01
        41      0.932E+02    *  0.384E-02      0.01788
        41      0.932E+02       0.218E-01                  0         7   0.169E+02
        42      0.931E+02    *  0.218E-01      0.00053
        42      0.931E+02       0.286E-01                  0         7   0.206E+02
        42      0.931E+02    *  0.286E-01      0.00025
        42      0.931E+02    *  0.286E-01      0.00010
        42 Not accepted damping factor        0.00001     0         7
        42      0.931E+02       0.106E-01                  0         6   0.102E+02
        43      0.931E+02    *  0.106E-01      0.00135
        43      0.931E+02       0.162E-01                  0         7   0.102E+02
        43      0.931E+02    *  0.163E-01      0.00054
        43      0.931E+02    *  0.162E-01      0.00012
        43      0.931E+02    *  0.162E-01      0.00010
        43 Not accepted damping factor        0.00002     0         7
        43      0.931E+02       0.102E-01                  0         6   0.769E+01
        44      0.930E+02    *  0.102E-01      0.00185
        44      0.930E+02       0.605E-01                  0         7   0.457E+02
        44      0.930E+02    *  0.606E-01      0.00010
        44 Not accepted damping factor        0.00001     0         7
        44      0.930E+02       0.726E-02                  0         6   0.678E+01
        45      0.935E+02    *  0.723E-02      0.00337
        45      0.935E+02       0.475E+00                  0         7   0.944E+03
        45      0.935E+02    *  0.480E+00      0.00010
        45 Not accepted damping factor        0.00000     0         7
        45      0.935E+02       0.646E-02                  0         6   0.468E+01
        46      0.931E+02    *  0.645E-02      0.00452
        46      0.931E+02       0.407E+00                  0         7   0.435E+03
        46      0.931E+02    *  0.409E+00      0.00010
        46 Not accepted damping factor        0.00000     0         7
        46      0.931E+02       0.768E-02                  0         6   0.554E+01
        47      0.933E+02    *  0.766E-02      0.00279
        47      0.933E+02       0.122E-01                  0         7   0.851E+01
        48      0.932E+02    *  0.122E-01      0.00087
        48      0.932E+02       0.657E-01                  0         7   0.439E+02
        48      0.932E+02    *  0.657E-01      0.00010
        48 Not accepted damping factor        0.00001     0         7
        48      0.932E+02       0.797E-02                  0         6   0.724E+01
        49      0.934E+02    *  0.796E-02      0.00098
        49      0.934E+02       0.431E-01                  0         7   0.579E+02
        50      0.934E+02    *  0.431E-01      0.00010
        50      0.934E+02       0.120E+00                  0         7   0.849E+02
        50      0.933E+02    *  0.120E+00      0.00010
        50 Not accepted damping factor        0.00000     0         7
        50      0.934E+02       0.134E-01                  0         6   0.155E+02
        51      0.933E+02    *  0.134E-01      0.00032
        51      0.933E+02       0.145E-01                  0         7   0.997E+01
        52      0.933E+02    *  0.145E-01      0.00038
        52      0.933E+02       0.135E+00                  0         7   0.553E+02
        52      0.933E+02    *  0.135E+00      0.00010
        52 Not accepted damping factor        0.00001     0         7
        52      0.933E+02       0.621E-02                  0         6   0.224E+02
        53      0.935E+02    *  0.620E-02      0.00143
        53      0.935E+02       0.195E-01                  0         7   0.135E+02
        54      0.935E+02    *  0.194E-01      0.00017
        54      0.935E+02       0.171E+00                  0         7   0.745E+02
        55      0.935E+02    *  0.171E+00      0.00010
        55      0.935E+02       0.163E-01                  0         7   0.130E+02
        56      0.934E+02    *  0.162E-01      0.00106
        56      0.934E+02       0.307E+00                  0         7   0.142E+03
        56      0.934E+02    *  0.308E+00      0.00010
        56 Not accepted damping factor        0.00000     0         7
        56      0.934E+02       0.873E-02                  0         6   0.144E+02
        57      0.937E+02    *  0.871E-02      0.00202
        57      0.937E+02       0.329E+00                  0         7   0.325E+03
        57      0.937E+02    *  0.330E+00      0.00010
        57 Not accepted damping factor        0.00000     0         7
        57      0.937E+02       0.924E-02                  0         6   0.656E+01
        58      0.936E+02    *  0.921E-02      0.00349
        58      0.936E+02       0.130E+00                  0         7   0.621E+02
        58      0.936E+02    *  0.130E+00      0.00010
        58 Not accepted damping factor        0.00001     0         7
        58      0.936E+02       0.447E-01                  0         6   0.390E+02
        59      0.936E+02    *  0.447E-01      0.00015
        59      0.936E+02       0.268E-01                  0         7   0.215E+02
        60      0.936E+02    *  0.268E-01      0.00018
        60      0.936E+02       0.444E-01                  0         7   0.375E+02
        60      0.936E+02    *  0.444E-01      0.00010
        60 Not accepted damping factor        0.00002     0         7
        60      0.936E+02       0.124E-01                  0         6   0.957E+01
        61      0.936E+02    *  0.124E-01      0.00045
        61      0.936E+02       0.834E-01                  0         7   0.821E+02
        62      0.936E+02    *  0.833E-01      0.00010
        62      0.936E+02       0.455E-01                  0         7   0.733E+02
        62      0.936E+02    *  0.455E-01      0.00014
        62      0.936E+02    *  0.455E-01      0.00010
        62 Not accepted damping factor        0.00001     0         7
        62      0.936E+02       0.168E-01                  0         6   0.195E+02
        63      0.936E+02    *  0.168E-01      0.00085
        63      0.936E+02       0.198E-01                  0         7   0.173E+02
        64      0.935E+02    *  0.198E-01      0.00064
        64      0.935E+02       0.352E+00                  0         7   0.168E+03
        64      0.935E+02    *  0.352E+00      0.00010
        64 Not accepted damping factor        0.00000     0         7
        64      0.935E+02       0.843E-02                  0         6   0.116E+02
        65      0.938E+02    *  0.842E-02      0.00149
        65      0.938E+02       0.547E-01                  0         7   0.399E+02
        66      0.938E+02    *  0.547E-01      0.00010
        66      0.938E+02       0.124E+00                  0         7   0.590E+02
        66      0.938E+02    *  0.124E+00      0.00010
        66 Not accepted damping factor        0.00001     0         7
        66      0.938E+02       0.145E-01                  0         6   0.217E+02
        66      0.939E+02    *  0.145E-01      0.00036
        66 Not accepted damping factor        0.00007     0         6
        66      0.938E+02       0.551E-02                  0         5   0.120E+02
        67      0.943E+02    *  0.550E-02      0.00258
        67      0.943E+02       0.110E-01                  0         7   0.134E+02
        68      0.943E+02    *  0.109E-01      0.00069
        68      0.943E+02       0.217E-01                  0         7   0.143E+02
        68      0.942E+02    *  0.217E-01      0.00021
        68      0.943E+02    *  0.217E-01      0.00010
        68 Not accepted damping factor        0.00002     0         7
        68      0.943E+02       0.115E-01                  0         6   0.102E+02
        69      0.943E+02    *  0.115E-01      0.00209
        69      0.943E+02       0.140E-01                  0         7   0.861E+01
        70      0.941E+02    *  0.140E-01      0.00140
        70      0.941E+02       0.129E+00                  0         7   0.611E+02
        71      0.941E+02    *  0.129E+00      0.00010
        71      0.941E+02       0.103E+00                  0         7   0.181E+03
        72      0.941E+02    *  0.102E+00      0.00013
        72      0.941E+02       0.252E+00                  0         7   0.248E+03
        73      0.941E+02    *  0.251E+00      0.00010
        73      0.941E+02       0.166E-01                  0         7   0.124E+02
        74      0.940E+02    *  0.166E-01      0.00149
        74      0.940E+02       0.198E+01                  0         7   0.876E+03
        75      0.939E+02    *  0.198E+01      0.00010
        75      0.939E+02       0.280E-01                  0         7   0.272E+02
        76      0.933E+02    *  0.275E-01      0.00711
        76      0.933E+02       0.950E-01                  0         7   0.624E+02
        77      0.932E+02    *  0.947E-01      0.00069
        77      0.932E+02       0.234E-01                  0         7   0.166E+02
        78      0.930E+02    *  0.232E-01      0.00277
        78      0.930E+02       0.938E-01                  0         7   0.422E+02
        78      0.929E+02    *  0.938E-01      0.00018
        78      0.929E+02    *  0.938E-01      0.00010
        78 Not accepted damping factor        0.00001     0         7
        78      0.930E+02       0.745E-02                  0         6   0.186E+02
        79      0.944E+02    *  0.735E-02      0.00990
        79      0.944E+02       0.101E-01                  0         7   0.761E+01
        80      0.939E+02    *  0.100E-01      0.00599
        80      0.939E+02       0.517E-01                  0         7   0.346E+02
        80      0.939E+02    *  0.517E-01      0.00023
        80      0.939E+02    *  0.517E-01      0.00010
        80 Not accepted damping factor        0.00001     0         7
        80      0.939E+02       0.106E-01                  0         6   0.888E+01
        81      0.947E+02    *  0.105E-01      0.00545
        81      0.947E+02       0.304E-01                  0         7   0.412E+02
        82      0.947E+02    *  0.304E-01      0.00054
        82      0.947E+02       0.140E-01                  0         7   0.204E+02
        83      0.945E+02    *  0.140E-01      0.00167
        83      0.945E+02       0.330E-01                  0         7   0.408E+02
        83      0.945E+02    *  0.330E-01      0.00035
        83      0.945E+02    *  0.330E-01      0.00010
        83 Not accepted damping factor        0.00001     0         7
        83      0.945E+02       0.116E-01                  0         6   0.664E+01
        84      0.942E+02    *  0.116E-01      0.00249
        84      0.942E+02       0.597E-01                  0         7   0.698E+02
        84      0.942E+02    *  0.597E-01      0.00010
        84 Not accepted damping factor        0.00001     0         7
        84      0.942E+02       0.701E-02                  0         6   0.108E+02
        85      0.945E+02    *  0.699E-02      0.00664
        85      0.945E+02       0.255E-01                  0         7   0.275E+02
        86      0.944E+02    *  0.255E-01      0.00052
        86      0.944E+02       0.286E-01                  0         7   0.244E+02
        86      0.944E+02    *  0.286E-01      0.00035
        86      0.944E+02    *  0.286E-01      0.00010
        86 Not accepted damping factor        0.00001     0         7
        86      0.944E+02       0.929E-02                  0         6   0.117E+02
        87      0.944E+02    *  0.928E-02      0.00175
        87      0.944E+02       0.891E-02                  0         7   0.786E+01
        88      0.942E+02    *  0.890E-02      0.00221
        88      0.942E+02       0.929E+02                  0         7   0.497E+05
        88      0.942E+02    *  0.130E+03      0.00010
        88 Not accepted damping factor        0.00000     0         7
        88      0.942E+02       0.549E-02                  0         6   0.613E+01
        89      0.949E+02    *  0.547E-02      0.00506
        89      0.949E+02       0.109E-01                  0         7   0.967E+01
        90      0.948E+02    *  0.109E-01      0.00150
        90      0.948E+02       0.633E-01                  0         7   0.401E+02
        90      0.948E+02    *  0.633E-01      0.00010
        90 Not accepted damping factor        0.00001     0         7
        90      0.948E+02       0.117E-01                  0         6   0.109E+02
        91      0.950E+02    *  0.117E-01      0.00105
        91      0.950E+02       0.233E-01                  0         7   0.208E+02
        92      0.950E+02    *  0.233E-01      0.00022
        92      0.950E+02       0.491E-01                  0         7   0.318E+02
        92      0.949E+02    *  0.491E-01      0.00010
        92 Not accepted damping factor        0.00001     0         7
        92      0.950E+02       0.130E-01                  0         6   0.171E+02
        93      0.950E+02    *  0.130E-01      0.00032
        93      0.950E+02       0.463E-01                  0         7   0.411E+02
        94      0.950E+02    *  0.463E-01      0.00010
        94      0.950E+02       0.621E-01                  0         7   0.486E+02
        94      0.950E+02    *  0.621E-01      0.00010
        94 Not accepted damping factor        0.00003     0         7
        94      0.950E+02       0.230E-01                  0         6   0.191E+02
        94      0.950E+02    *  0.230E-01      0.00017
        94 Not accepted damping factor        0.00002     0         6
        94      0.950E+02       0.518E-02                  0         5   0.571E+01
        95      0.953E+02    *  0.518E-02      0.00157
        95      0.953E+02       0.122E-01                  0         7   0.172E+02
        96      0.953E+02    *  0.122E-01      0.00027
        96      0.953E+02       0.166E-01                  0         7   0.142E+02
        96      0.952E+02    *  0.166E-01      0.00020
        96      0.953E+02    *  0.166E-01      0.00010
        96 Not accepted damping factor        0.00002     0         7
        96      0.953E+02       0.119E-01                  0         6   0.103E+02
        97      0.952E+02    *  0.119E-01      0.00072
        97      0.952E+02       0.140E-01                  0         7   0.809E+01
        97      0.952E+02    *  0.140E-01      0.00046
        97      0.952E+02    *  0.140E-01      0.00016
        97      0.952E+02    *  0.140E-01      0.00010
        97 Not accepted damping factor        0.00003     0         7
        97      0.952E+02       0.121E-01                  0         6   0.669E+01
        97      0.952E+02    *  0.121E-01      0.00054
        97      0.952E+02    *  0.121E-01      0.00017
        97 Not accepted damping factor        0.00005     0         6
        97      0.952E+02       0.803E-02                  0         5   0.548E+01
        98      0.953E+02    *  0.803E-02      0.00107
        98      0.953E+02       0.276E-01                  0         7   0.279E+02
        98      0.953E+02    *  0.276E-01      0.00010
        98 Not accepted damping factor        0.00001     0         7
        98      0.953E+02       0.127E-01                  0         6   0.133E+02
        98      0.953E+02    *  0.127E-01      0.00032
        98      0.953E+02    *  0.127E-01      0.00010
        98 Not accepted damping factor        0.00003     0         6
        98      0.953E+02       0.391E-02                  0         5   0.283E+01
        99      0.953E+02    *  0.390E-02      0.00447
        99      0.953E+02       0.559E-01                  0         7   0.441E+02
       100      0.953E+02    *  0.559E-01      0.00010

 Iteration terminates after NITMAX =100  Iteration steps

   Subcondition ( 1,   7)  0.441E+02
   Sensitivity ( 1,   7)  0.220E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :     101  ***
   ***  Corrector steps   :      85  ***
   ***  Rejected rk-1 st. :      56  ***
   ***  Jacobian eval.    :     100  ***
   ***  Function eval.    :     186  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.9532E+02

 problem: Typical Rect. Singly Reinforced Concrete Beam   


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       413 (  0.6 percent)

 Integer Workspace declared as      1000 is used up to        65 (  6.5 percent)



 N =  13

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by numerical differentiation (without feedback strategy)
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =    13
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.332E+06       0.622E+06                  0        13   0.294E+05
         0      0.287E+08    *  0.187E+08      0.01000
         1      0.330E+06    *  0.622E+06      0.00010
         1      0.330E+06       0.634E+04                  0        13   0.362E+03
         2      0.321E+06    *  0.634E+04      0.00040
         2      0.321E+06       0.803E+03                  0        13   0.762E+02
         3      0.297E+06    *  0.801E+03      0.00159
         3      0.297E+06       0.170E+03                  0        13   0.291E+02
         4      0.254E+06    *  0.169E+03      0.00620
         4      0.254E+06       0.457E+02                  0        13   0.144E+02
         5      0.211E+06    *  0.442E+02      0.02450
         5      0.211E+06       0.157E+02                  0        13   0.886E+01
         6      0.182E+06    *  0.138E+02      0.11527
         6      0.182E+06       0.675E+01                  0        13   0.840E+01
         7      0.180E+06    *  0.501E+01      0.34314
         7      0.180E+06       0.343E+01                  0        13   0.766E+01
         8      0.117E+06    *  0.285E+01      0.21304
         8      0.117E+06       0.217E+01                  0        13   0.795E+01
         9      0.759E+05    *  0.169E+01      0.28477
         9      0.759E+05       0.161E+01                  0        13   0.870E+01
        10      0.512E+05    *  0.117E+01      0.34429
        10      0.512E+05       0.113E+01                  0        13   0.932E+01
        11      0.325E+05    *  0.724E+00      0.45813
        11      0.325E+05       0.744E+00                  0        13   0.988E+01
        12      0.104E+05    *  0.407E+00      0.59517
        12      0.104E+05       0.495E+00                  0        13   0.865E+01
        13      0.390E+05    *  0.274E+00      0.61427
        13      0.390E+05       0.389E+00                  0        13   0.803E+01
        14      0.673E+05    *  0.228E+00      0.56275
        14      0.673E+05       0.335E+00                  0        13   0.773E+01
        15      0.945E+05    *  0.188E+00      0.60398
        15      0.945E+05       0.289E+00                  0        13   0.796E+01
        16      0.118E+06    *  0.148E+00      0.66937
        16      0.118E+06       0.232E+00                  0        13   0.883E+01
        17      0.135E+06    *  0.956E-01      0.81121
        17      0.135E+06       0.155E+00                  0        13   0.107E+02
        18      0.107E+06    *  0.403E-01      1.00000
        18      0.107E+06       0.659E-01                  0        13   0.134E+02
        19      0.270E+05    *  0.101E-01      1.00000
        19      0.270E+05       0.134E-01                  0        13   0.151E+02
        20      0.143E+04    *  0.559E-03      1.00000
        20      0.143E+04       0.601E-03                  0        13   0.159E+02
        21      0.313E+01    *  0.121E-05      1.00000
        21      0.313E+01       0.122E-05                  0        13   0.160E+02
        22      0.133E-04    *  0.497E-11      1.00000



 Solution of nonlinear system of equations obtained within  22 iteration steps

 Achieved relative accuracy 0.497E-11

   Subcondition ( 1,  13)  0.160E+02
   Sensitivity ( 1,  13)  0.150E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      22  ***
   ***  Corrector steps   :       1  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      22  ***
   ***  Function eval.    :      24  ***
   ***  ...  for Jacobian :     286  ***
   *************************************

 norm of residuum =   0.2583E-09

 problem: Wallis Function of 1685                         


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to        77 (  0.1 percent)

 Integer Workspace declared as      1000 is used up to        53 (  5.3 percent)



 N =   1

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by a user subroutine
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     1
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.609E+01       0.354E+03                  0         1   0.100E+01
         0      0.692E+02    *  0.402E+04      0.01000
         1      0.608E+01    *  0.354E+03      0.00010
         1      0.608E+01       0.312E+02                  0         1   0.100E+01
         2      0.607E+01    *  0.311E+02      0.00124
         2      0.607E+01       0.152E+02                  0         1   0.100E+01
         3      0.603E+01    *  0.151E+02      0.00499
         3      0.603E+01       0.731E+01                  0         1   0.100E+01
         4      0.584E+01    *  0.709E+01      0.02014
         4      0.584E+01       0.320E+01                  0         1   0.100E+01
         5      0.507E+01    *  0.278E+01      0.08341
         5      0.507E+01       0.105E+01                  0         1   0.100E+01
         6      0.188E+01    *  0.389E+00      0.38760
         6      0.188E+01       0.128E+00                  0         1   0.100E+01
         7      0.264E+00    *  0.180E-01      1.00000
         7      0.264E+00       0.115E-01                  0         1   0.100E+01
         8      0.337E-02    *  0.146E-03      1.00000
         8      0.337E-02       0.143E-03                  0         1   0.100E+01
         9      0.572E-06    *  0.243E-07      1.00000
         9      0.572E-06       0.245E-07                  0         1   0.100E+01
        10      0.195E-13    *  0.836E-15      1.00000



 Solution of nonlinear system of equations obtained within  10 iteration steps

 Achieved relative accuracy 0.836E-15

   Subcondition ( 1,   1)  0.100E+01
   Sensitivity ( 1,   1)  0.100E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      10  ***
   ***  Corrector steps   :       1  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :      10  ***
   ***  Function eval.    :      12  ***
   ***  ...  for Jacobian :       0  ***
   *************************************

 norm of residuum =   0.1776E-14

 problem: Inters. of Ellipsoid with a Hyperboloid - red.  


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       113 (  0.2 percent)

 Integer Workspace declared as      1000 is used up to        55 (  5.5 percent)



 N =   3

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by numerical differentiation (without feedback strategy)
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     3
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.719E+01       0.154E+01                  0         3   0.253E+01
         1      0.712E+01    *  0.152E+01      0.01000
         1      0.712E+01       0.148E+01                  0         3   0.248E+01
         2      0.645E+00    *  0.193E+00      0.57413
         2      0.645E+00       0.110E+00                  0         3   0.133E+01
         3      0.114E+00    *  0.651E-02      1.00000
         3      0.114E+00       0.423E-02                  0         3   0.152E+01
         4      0.398E-03    *  0.150E-04      1.00000
         4      0.398E-03       0.146E-04                  0         3   0.153E+01
         5      0.514E-08    *  0.187E-09      1.00000
         5      0.514E-08       0.188E-09                  0         3   0.153E+01
         6      0.657E-15    *  0.223E-15      1.00000



 Solution of nonlinear system of equations obtained within   6 iteration steps

 Achieved relative accuracy 0.223E-15

   Subcondition ( 1,   3)  0.153E+01
   Sensitivity ( 1,   3)  0.140E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       6  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       6  ***
   ***  Function eval.    :       7  ***
   ***  ...  for Jacobian :      18  ***
   *************************************

 norm of residuum =   0.5813E-15

 problem: Esterification Reaction - reduced system        


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to        94 (  0.1 percent)

 Integer Workspace declared as      1000 is used up to        54 (  5.4 percent)



 N =   2

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by numerical differentiation (without feedback strategy)
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     2
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.971E+07       0.412E+01                  0         2   0.596E+03
         1      0.962E+07    *  0.408E+01      0.01000
         1      0.962E+07       0.428E+01                  0         2   0.619E+03
         2      0.549E+07    *  0.285E+01      0.50172
         2      0.549E+07       0.421E+01                  0         2   0.498E+03
         3      0.383E+07    *  0.253E+01      0.33745
         3      0.383E+07       0.201E+01                  0         2   0.131E+03
         4      0.959E+06    *  0.111E+01      1.00000
         4      0.959E+06       0.263E+00                  0         2   0.322E+02
         5      0.240E+06    *  0.990E-01      1.00000
         5      0.240E+06       0.178E+00                  0         2   0.150E+02
         6      0.729E+05    *  0.720E-01      0.89746
         6      0.729E+05       0.196E+00                  0         2   0.960E+01
         7      0.244E+05    *  0.862E-01      0.84282
         7      0.244E+05       0.230E+00                  0         2   0.532E+01
         8      0.950E+04    *  0.114E+00      0.75224
         8      0.950E+04       0.300E+00                  0         2   0.109E+01
         9      0.450E+04    *  0.174E+00      0.62349
         9      0.450E+04       0.467E+00                  0         2   0.760E+01
        10      0.274E+04    *  0.328E+00      0.43968
        10      0.274E+04       0.927E+00                  0         2   0.206E+02
        11      0.213E+04    *  0.780E+00      0.23632
        11      0.213E+04       0.248E+01                  0         2   0.537E+02
        12      0.196E+04    *  0.234E+01      0.08318
        12      0.196E+04       0.893E+01                  0         2   0.164E+03
        13      0.193E+04    *  0.884E+01      0.01652
        13      0.193E+04       0.440E+02                  0         2   0.675E+03
        14      0.192E+04    *  0.440E+02      0.00146
        14      0.192E+04       0.358E+03                  0         2   0.495E+04
        14      0.192E+04    *  0.358E+03      0.00010
        14 Not accepted damping factor        0.00004     0         2
        14      0.192E+04       0.124E-01                  0         1   0.100E+01
        15      0.272E+04    *  0.224E-03      1.00000
        15      0.272E+04       0.656E+02                  0         2   0.866E+03
        16      0.272E+04    *  0.656E+02      0.00010
        16      0.272E+04       0.730E+02                  0         2   0.950E+03
        17      0.272E+04    *  0.730E+02      0.00082
        17      0.272E+04       0.216E+04                  0         2   0.274E+05
        17      0.272E+04    *  0.217E+04      0.00010
        17 Not accepted damping factor        0.00000     0         2
        17      0.272E+04       0.664E-02                  0         1   0.100E+01
        18      0.317E+04    *  0.237E-04      1.00000
        18      0.317E+04       0.153E+03                  0         2   0.188E+04
        19      0.317E+04    *  0.153E+03      0.00010
        19      0.317E+04       0.290E+03                  0         2   0.351E+04
        20      0.317E+04    *  0.290E+03      0.00010
        20      0.317E+04       0.493E+03                  0         2   0.588E+04
        20      0.317E+04    *  0.493E+03      0.00010
        20 Not accepted damping factor        0.00002     0         2
        20      0.317E+04       0.350E-02                  0         1   0.100E+01
        21      0.340E+04    *  0.196E-05      1.00000
        21      0.340E+04       0.113E+04                  0         2   0.133E+05
        21      0.340E+04    *  0.113E+04      0.00010
        21 Not accepted damping factor        0.00000     0         2
        21      0.340E+04       0.180E-02                  0         1   0.100E+01
        22      0.340E+04    *  0.180E-02      0.00217
        22      0.340E+04       0.113E+04                  0         2   0.132E+05
        22      0.340E+04    *  0.113E+04      0.00010
        22 Not accepted damping factor        0.00000     0         2
        22      0.340E+04       0.180E-02                  0         1   0.100E+01
        23      0.352E+04    *  0.391E-07      1.00000
        23      0.352E+04       0.412E+03                  0         2   0.481E+04
        23      0.352E+04    *  0.412E+03      0.00010
        23 Not accepted damping factor        0.00003     0         2
        23      0.352E+04       0.916E-03                  0         1   0.100E+01
        24      0.352E+04    *  0.916E-03      0.00010
        24      0.352E+04       0.412E+03                  0         2   0.480E+04
        24      0.352E+04    *  0.412E+03      0.00010
        24 Not accepted damping factor        0.00003     0         2
        24      0.352E+04       0.917E-03                  0         1   0.100E+01
        25      0.359E+04    *  0.823E-07      1.00000
        25      0.359E+04       0.310E+03                  0         2   0.360E+04
        26      0.359E+04    *  0.310E+03      0.00010
        26      0.359E+04       0.422E+03                  0         2   0.485E+04
        26      0.359E+04    *  0.423E+03      0.00010
        26 Not accepted damping factor        0.00003     0         2
        26      0.359E+04       0.414E-03                  0         1   0.100E+01
        27      0.362E+04    *  0.229E-07      1.00000
        27      0.362E+04       0.500E+03                  0         2   0.567E+04
        27      0.361E+04    *  0.500E+03      0.00010
        27 Not accepted damping factor        0.00002     0         2
        27      0.362E+04       0.214E-03                  0         1   0.100E+01
        28      0.362E+04    *  0.214E-03      0.00021
        28      0.362E+04       0.500E+03                  0         2   0.567E+04
        28      0.361E+04    *  0.500E+03      0.00010
        28 Not accepted damping factor        0.00002     0         2
        28      0.362E+04       0.214E-03                  0         1   0.100E+01
        29      0.363E+04    *  0.686E-08      1.00000
        29      0.363E+04       0.546E+03                  0         2   0.618E+04
        29      0.363E+04    *  0.546E+03      0.00010
        29 Not accepted damping factor        0.00002     0         2
        29      0.363E+04       0.109E-03                  0         1   0.100E+01
        30      0.363E+04    *  0.109E-03      0.00012
        30      0.363E+04       0.546E+03                  0         2   0.618E+04
        30      0.363E+04    *  0.546E+03      0.00010
        30 Not accepted damping factor        0.00002     0         2
        30      0.363E+04       0.109E-03                  0         1   0.100E+01
        31      0.364E+04    *  0.188E-08      1.00000
        31      0.364E+04       0.572E+03                  0         2   0.648E+04
        31      0.364E+04    *  0.572E+03      0.00010
        31 Not accepted damping factor        0.00002     0         2
        31      0.364E+04       0.551E-04                  0         1   0.100E+01
        32      0.364E+04    *  0.551E-04      0.00010
        32      0.364E+04       0.572E+03                  0         2   0.648E+04
        32      0.364E+04    *  0.572E+03      0.00010
        32 Not accepted damping factor        0.00002     0         2
        32      0.364E+04       0.551E-04                  0         1   0.100E+01
        33      0.364E+04    *  0.487E-09      1.00000
        33      0.364E+04       0.587E+03                  0         2   0.665E+04
        33      0.364E+04    *  0.587E+03      0.00010
        33 Not accepted damping factor        0.00002     0         2
        33      0.364E+04       0.280E-04                  0         1   0.100E+01
        34      0.364E+04    *  0.280E-04      0.00010
        34      0.364E+04       0.587E+03                  0         2   0.665E+04
        34      0.364E+04    *  0.587E+03      0.00010
        34 Not accepted damping factor        0.00002     0         2
        34      0.364E+04       0.280E-04                  0         1   0.100E+01
        35      0.364E+04    *  0.136E-09      1.00000
        35      0.364E+04       0.595E+03                  0         2   0.673E+04
        35      0.364E+04    *  0.595E+03      0.00010
        35 Not accepted damping factor        0.00002     0         2
        35      0.364E+04       0.142E-04                  0         1   0.100E+01
        36      0.364E+04    *  0.142E-04      0.00010
        36      0.364E+04       0.595E+03                  0         2   0.673E+04
        36      0.364E+04    *  0.595E+03      0.00010
        36 Not accepted damping factor        0.00002     0         2
        36      0.364E+04       0.142E-04                  0         1   0.100E+01

 Iteration terminates at stationary point


   Subcondition ( 1,   1)  0.100E+01
   Sensitivity ( 1,   1)  0.141E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :      37  ***
   ***  Corrector steps   :      18  ***
   ***  Rejected rk-1 st. :      18  ***
   ***  Jacobian eval.    :      37  ***
   ***  Function eval.    :      56  ***
   ***  ...  for Jacobian :      74  ***
   *************************************

 norm of residuum =   0.3644E+04

 problem: Gupta problem - reduced system                  


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to        94 (  0.1 percent)

 Integer Workspace declared as      1000 is used up to        54 (  5.4 percent)



 N =   2

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by numerical differentiation (without feedback strategy)
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     2
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.177E+04       0.417E+00                  0         2   0.149E+01
         1      0.176E+04    *  0.413E+00      0.01000
         1      0.176E+04       0.374E+00                  0         2   0.144E+01
         2      0.155E+04    *  0.269E+00      0.11515
         2      0.155E+04       0.765E-01                  0         2   0.107E+01
         2 FCN could not be evaluated          0.70747     0         2
         3      0.992E+03    *  0.385E-01      0.35374
         3      0.992E+03       0.384E-01                  0         2   0.100E+01
         4      0.155E+02    *  0.334E-02      1.00000
         4      0.155E+02       0.370E-02                  0         2   0.127E+01
         4 FCN could not be evaluated          1.00000     0         2
         5      0.115E+01    *  0.274E-03      0.50000
         5      0.115E+01       0.594E-04                  0         2   0.118E+01
         6      0.409E-01    *  0.211E-05      1.00000
         6      0.409E-01       0.227E-05                  0         2   0.114E+01
         7      0.554E-04    *  0.308E-08      1.00000
         7      0.554E-04       0.309E-08                  0         2   0.113E+01
         8      0.132E-09    *  0.736E-14      1.00000



 Solution of nonlinear system of equations obtained within   8 iteration steps

 Achieved relative accuracy 0.736E-14

   Subcondition ( 1,   2)  0.113E+01
   Sensitivity ( 1,   2)  0.115E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       8  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       8  ***
   ***  Function eval.    :      11  ***
   ***  ...  for Jacobian :      16  ***
   *************************************

 norm of residuum =   0.2247E-13

 problem: Dew Point Temperature - reduced system          


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       113 (  0.2 percent)

 Integer Workspace declared as      1000 is used up to        55 (  5.5 percent)



 N =   3

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by numerical differentiation (without feedback strategy)
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     3
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.132E+00       0.846E-01                  0         3   0.442E+01
         1      0.130E+00    *  0.838E-01      0.01000
         1      0.130E+00       0.821E-01                  0         3   0.451E+01
         2      0.658E-01    *  0.371E-01      0.42172
         2      0.658E-01       0.299E-01                  0         3   0.585E+01
         3      0.316E-02    *  0.152E-01      1.00000
         3      0.316E-02       0.232E-01                  0         3   0.499E+01
         4      0.594E-03    *  0.198E-02      1.00000
         4      0.594E-03       0.243E-02                  0         3   0.460E+01
         5      0.969E-05    *  0.319E-04      1.00000
         5      0.969E-05       0.328E-04                  0         3   0.454E+01
         6      0.215E-08    *  0.682E-08      1.00000
         6      0.215E-08       0.683E-08                  0         3   0.454E+01
         7      0.795E-15    *  0.694E-15      1.00000



 Solution of nonlinear system of equations obtained within   7 iteration steps

 Achieved relative accuracy 0.694E-15

   Subcondition ( 1,   3)  0.454E+01
   Sensitivity ( 1,   3)  0.142E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       7  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       7  ***
   ***  Function eval.    :       8  ***
   ***  ...  for Jacobian :      21  ***
   *************************************

 norm of residuum =   0.6410E-16

 problem: Part. Oxid. of Methane with Oxygen - reduced    


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to       113 (  0.2 percent)

 Integer Workspace declared as      1000 is used up to        55 (  5.5 percent)



 N =   3

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by numerical differentiation (without feedback strategy)
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     3
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.109E+05       0.317E+00                  0         3   0.279E+01
         1      0.107E+05    *  0.314E+00      0.01000
         1      0.107E+05       0.297E+00                  0         3   0.277E+01
         2      0.348E+04    *  0.983E-01      0.19876
         2      0.348E+04       0.240E-01                  0         3   0.228E+02
         3      0.309E+04    *  0.630E-02      1.00000
         3      0.309E+04       0.159E-02                  0         3   0.924E+02
         4      0.537E+03    *  0.372E-03      1.00000
         4      0.537E+03       0.544E-03                  0         3   0.620E+02
         5      0.237E+02    *  0.210E-04      1.00000
         5      0.237E+02       0.229E-04                  0         3   0.565E+02
         6      0.494E-01    *  0.465E-07      1.00000
         6      0.494E-01       0.467E-07                  0         3   0.562E+02
         7      0.170E-06    *  0.162E-12      1.00000



 Solution of nonlinear system of equations obtained within   7 iteration steps

 Achieved relative accuracy 0.162E-12

   Subcondition ( 1,   3)  0.562E+02
   Sensitivity ( 1,   3)  0.173E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       7  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       7  ***
   ***  Function eval.    :       8  ***
   ***  ...  for Jacobian :      21  ***
   *************************************

 norm of residuum =   0.1050E-11

 problem: Sulphur Dioxide to Sulphur Trioxide - reduced   


    N L E Q 2  *****  V e r s i o n  2 . 3 ***

 Newton-Method for the solution of nonlinear systems


 Real    Workspace declared as     70000 is used up to        77 (  0.1 percent)

 Integer Workspace declared as      1000 is used up to        53 (  5.3 percent)



 N =   1

 Prescribed relative precision  0.10E-09

 The Jacobian is supplied by numerical differentiation (without feedback strategy)
 Automatic row scaling of the Jacobian is allowed

 Rank-1 updates are inhibited
 Problem is specified as being highly nonlinear
 Bounded damping strategy is off
 Maximum permitted number of iteration steps :    100


 Internal parameters:

 Starting value for damping factor FCSTART =  0.10E-01
 Minimum allowed damping factor FCMIN =  0.10E-03
 Rank-1 updates decision parameter SIGMA =  0.10E+06
 Initial Jacobian pseudo-rank IRANK =     1
 Maximum permitted subcondition COND =  0.10E+18



  ******************************************************************

        It       Normf           Normx         Damp.Fct.   New      Rank        Cond
         0      0.717E-01       0.378E+00                  0         1   0.100E+01
         1      0.710E-01    *  0.374E+00      0.01000
         1      0.710E-01       0.357E+00                  0         1   0.100E+01
         2      0.398E-01    *  0.201E+00      0.23750
         2      0.398E-01       0.634E-01                  0         1   0.100E+01
         3      0.265E-01    *  0.422E-01      1.00000
         3      0.265E-01       0.163E-01                  0         1   0.100E+01
         4      0.305E-02    *  0.187E-02      1.00000
         4      0.305E-02       0.240E-02                  0         1   0.100E+01
         5      0.552E-04    *  0.434E-04      1.00000
         5      0.552E-04       0.450E-04                  0         1   0.100E+01
         6      0.189E-07    *  0.154E-07      1.00000
         6      0.189E-07       0.154E-07                  0         1   0.100E+01
         7      0.376E-14    *  0.307E-14      1.00000



 Solution of nonlinear system of equations obtained within   7 iteration steps

 Achieved relative accuracy 0.307E-14

   Subcondition ( 1,   1)  0.100E+01
   Sensitivity ( 1,   1)  0.100E+01


   ******  Statistics * NLEQ2    *******
   ***  Newton iterations :       7  ***
   ***  Corrector steps   :       0  ***
   ***  Rejected rk-1 st. :       0  ***
   ***  Jacobian eval.    :       7  ***
   ***  Function eval.    :       8  ***
   ***  ...  for Jacobian :       7  ***
   *************************************

 norm of residuum =   0.5551E-16

