pcprb3.out.txt

求解非线性方程组的一组源代码,FORTRAN90.用于解决N个未知数,N-1个方程.

  
 Jacobians:       0
 Factorizations:  55
 Solves:          77
 Functions:       2500
  
 This is run number  4
  
 This run uses the banded linear solver BANSLV.
 This run assumes that the user supplies the
 jacobian matrix via a subroutine.
  
 Step  Type of point     Lambda
  
   0  Start point      0.00000    
  
 PITCON 7.0
   University of Pittsburgh continuation code
  
   Modified:        12 November 1999
   Linear algebra:  LAPACK
   Precision:       double
  
   1  Corrected        0.00000    
 TRYSTP - DEBUG - about to return.
   2  Continuation    0.264435E-01
 TRYSTP - DEBUG - about to return.
   3  Continuation    0.104600    
 TRYSTP - DEBUG - about to return.
   4  Continuation    0.219655    
 TRYSTP - DEBUG - about to return.
   5  Continuation    0.327553    
 TRYSTP - DEBUG - about to return.
   6  Continuation    0.425474    
 TRYSTP - DEBUG - about to return.
   7  Continuation    0.513404    
 TRYSTP - DEBUG - about to return.
   8  Continuation    0.591468    
 TRYSTP - DEBUG - about to return.
   9  Continuation    0.659916    
 TRYSTP - DEBUG - about to return.
  10  Continuation    0.719109    
 TRYSTP - DEBUG - about to return.
  11  Continuation    0.769501    
 TRYSTP - DEBUG - about to return.
  12  Continuation    0.811615    
  13  Target point    0.800000    
  
 Value of target point  1
    0.00000      0.602686E-01  0.118413      0.174305      0.227816    
   0.278815      0.327169      0.372748      0.415422      0.455064    
   0.491551      0.524765      0.554597      0.580942      0.603707    
   0.622811      0.638182      0.649762      0.657507      0.661387    
   0.661387      0.800000    
  
 TRYSTP - DEBUG - about to return.
  14  Continuation    0.846018    
 TRYSTP - DEBUG - about to return.
  15  Continuation    0.873302    
 TRYSTP - DEBUG - about to return.
  16  Continuation    0.894067    
 TRYSTP - DEBUG - about to return.
  17  Continuation    0.908904    
 TRYSTP - DEBUG - about to return.
  18  Continuation    0.918385    
 TRYSTP - DEBUG - about to return.
  19  Continuation    0.923057    
 TRYSTP - DEBUG - about to return.
  20  Continuation    0.923436    
  21  Limit point     0.923283    
 TRYSTP - DEBUG - about to return.
  22  Continuation    0.920003    
 TRYSTP - DEBUG - about to return.
  23  Continuation    0.913206    
 TRYSTP - DEBUG - about to return.
  24  Continuation    0.904857    
 TRYSTP - DEBUG - about to return.
  25  Continuation    0.891133    
 TRYSTP - DEBUG - about to return.
  26  Continuation    0.877626    
 TRYSTP - DEBUG - about to return.
  27  Continuation    0.860121    
 TRYSTP - DEBUG - about to return.
  28  Continuation    0.842782    
 TRYSTP - DEBUG - about to return.
  29  Continuation    0.822585    
 TRYSTP - DEBUG - about to return.
  30  Continuation    0.802513    
 TRYSTP - DEBUG - about to return.
  31  Continuation    0.780519    
  32  Target point    0.800000    
  
 Value of target point  2
    0.00000      0.152806      0.303281      0.451047      0.595672    
   0.736669      0.873486       1.00551       1.13207       1.25242    
    1.36578       1.47130       1.56811       1.65533       1.73208    
    1.79754       1.85093       1.89160       1.91902       1.93283    
    1.93283      0.800000    
  
  
 Jacobians:       55
 Factorizations:  55
 Solves:          130
 Functions:       80
  
 This is run number  5
  
 This run uses the banded linear solver BANSLV.
 This run assumes that PITCON will approximate
 the jacobian using forward differences.
  
 Step  Type of point     Lambda
  
   0  Start point      0.00000    
  
 PITCON 7.0
   University of Pittsburgh continuation code
  
   Modified:        12 November 1999
   Linear algebra:  LAPACK
   Precision:       double
  
   1  Corrected        0.00000    
 TRYSTP - DEBUG - about to return.
   2  Continuation    0.264435E-01
 TRYSTP - DEBUG - about to return.
   3  Continuation    0.104600    
 TRYSTP - DEBUG - about to return.
   4  Continuation    0.219655    
 TRYSTP - DEBUG - about to return.
   5  Continuation    0.327553    
 TRYSTP - DEBUG - about to return.
   6  Continuation    0.425474    
 TRYSTP - DEBUG - about to return.
   7  Continuation    0.513405    
 TRYSTP - DEBUG - about to return.
   8  Continuation    0.591469    
 TRYSTP - DEBUG - about to return.
   9  Continuation    0.659917    
 TRYSTP - DEBUG - about to return.
  10  Continuation    0.719110    
 TRYSTP - DEBUG - about to return.
  11  Continuation    0.769503    
 TRYSTP - DEBUG - about to return.
  12  Continuation    0.811616    
  13  Target point    0.800000    
  
 Value of target point  1
    0.00000      0.602686E-01  0.118413      0.174305      0.227816    
   0.278815      0.327169      0.372749      0.415422      0.455064    
   0.491551      0.524765      0.554597      0.580942      0.603707    
   0.622811      0.638182      0.649762      0.657508      0.661388    
   0.661388      0.800000    
  
 TRYSTP - DEBUG - about to return.
  14  Continuation    0.846019    
 TRYSTP - DEBUG - about to return.
  15  Continuation    0.873303    
 TRYSTP - DEBUG - about to return.
  16  Continuation    0.894068    
 TRYSTP - DEBUG - about to return.
  17  Continuation    0.908904    
 TRYSTP - DEBUG - about to return.
  18  Continuation    0.918385    
 TRYSTP - DEBUG - about to return.
  19  Continuation    0.923057    
 TRYSTP - DEBUG - about to return.
  20  Continuation    0.923436    
  21  Limit point     0.923282    
 TRYSTP - DEBUG - about to return.
  22  Continuation    0.920003    
 TRYSTP - DEBUG - about to return.
  23  Continuation    0.913205    
 TRYSTP - DEBUG - about to return.
  24  Continuation    0.904848    
 TRYSTP - DEBUG - about to return.
  25  Continuation    0.891132    
 TRYSTP - DEBUG - about to return.
  26  Continuation    0.877616    
 TRYSTP - DEBUG - about to return.
  27  Continuation    0.860119    
 TRYSTP - DEBUG - about to return.
  28  Continuation    0.842772    
 TRYSTP - DEBUG - about to return.
  29  Continuation    0.822583    
 TRYSTP - DEBUG - about to return.
  30  Continuation    0.802502    
 TRYSTP - DEBUG - about to return.
  31  Continuation    0.780517    
  32  Target point    0.800000    
  
 Value of target point  2
    0.00000      0.152804      0.303277      0.451041      0.595665    
   0.736659      0.873474       1.00550       1.13206       1.25241    
    1.36576       1.47128       1.56809       1.65530       1.73206    
    1.79751       1.85090       1.89157       1.91899       1.93280    
    1.93280      0.800000    
  
  
 Jacobians:       0
 Factorizations:  56
 Solves:          132
 Functions:       361
  
 This is run number  6
  
 This run uses the banded linear solver BANSLV.
 This run assumes that PITCON will approximate
 the jacobian using central differences.
  
 Step  Type of point     Lambda
  
   0  Start point      0.00000    
  
 PITCON 7.0
   University of Pittsburgh continuation code
  
   Modified:        12 November 1999
   Linear algebra:  LAPACK
   Precision:       double
  
   1  Corrected        0.00000    
 TRYSTP - DEBUG - about to return.
   2  Continuation    0.264435E-01
 TRYSTP - DEBUG - about to return.
   3  Continuation    0.104600    
 TRYSTP - DEBUG - about to return.
   4  Continuation    0.219655    
 TRYSTP - DEBUG - about to return.
   5  Continuation    0.327553    
 TRYSTP - DEBUG - about to return.
   6  Continuation    0.425474    
 TRYSTP - DEBUG - about to return.
   7  Continuation    0.513404    
 TRYSTP - DEBUG - about to return.
   8  Continuation    0.591468    
 TRYSTP - DEBUG - about to return.
   9  Continuation    0.659916    
 TRYSTP - DEBUG - about to return.
  10  Continuation    0.719109    
 TRYSTP - DEBUG - about to return.
  11  Continuation    0.769501    
 TRYSTP - DEBUG - about to return.
  12  Continuation    0.811615    
  13  Target point    0.800000    
  
 Value of target point  1
    0.00000      0.602686E-01  0.118413      0.174305      0.227816    
   0.278815      0.327169      0.372748      0.415422      0.455064    
   0.491551      0.524765      0.554597      0.580942      0.603707    
   0.622811      0.638182      0.649762      0.657507      0.661387    
   0.661387      0.800000    
  
 TRYSTP - DEBUG - about to return.
  14  Continuation    0.846018    
 TRYSTP - DEBUG - about to return.
  15  Continuation    0.873302    
 TRYSTP - DEBUG - about to return.
  16  Continuation    0.894067    
 TRYSTP - DEBUG - about to return.
  17  Continuation    0.908904    
 TRYSTP - DEBUG - about to return.
  18  Continuation    0.918385    
 TRYSTP - DEBUG - about to return.
  19  Continuation    0.923057    
 TRYSTP - DEBUG - about to return.
  20  Continuation    0.923436    
  21  Limit point     0.923282    
 TRYSTP - DEBUG - about to return.
  22  Continuation    0.920003    
 TRYSTP - DEBUG - about to return.
  23  Continuation    0.913206    
 TRYSTP - DEBUG - about to return.
  24  Continuation    0.904857    
 TRYSTP - DEBUG - about to return.
  25  Continuation    0.891133    
 TRYSTP - DEBUG - about to return.
  26  Continuation    0.877626    
 TRYSTP - DEBUG - about to return.
  27  Continuation    0.860121    
 TRYSTP - DEBUG - about to return.
  28  Continuation    0.842782    
 TRYSTP - DEBUG - about to return.
  29  Continuation    0.822585    
 TRYSTP - DEBUG - about to return.
  30  Continuation    0.802513    
 TRYSTP - DEBUG - about to return.
  31  Continuation    0.780519    
  32  Target point    0.800000    
  
 Value of target point  2
    0.00000      0.152806      0.303281      0.451047      0.595672    
   0.736669      0.873486       1.00551       1.13207       1.25242    
    1.36578       1.47130       1.56811       1.65533       1.73208    
    1.79754       1.85093       1.89160       1.91902       1.93283    
    1.93283      0.800000    
  
  
 Jacobians:       0
 Factorizations:  57
 Solves:          134
 Functions:       538