pcprb3.out.txt

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

  
 PCPRB3
   PITCON sample program.
   Two point boundary value problem.
  
   Number of equations is  21
   Number of variables is  22
  
   Seek limit points in lambda.
  
   This problem will be run six times.
  
  
 This is run number  1
  
 This run uses the full linear solver DENSLV.
 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:       56
 Factorizations:  56
 Solves:          78
 Functions:       81
  
 This is run number  2
  
 This run uses the full linear solver DENSLV.
 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.923281    
 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:  57
 Solves:          79
 Functions:       1393
  
 This is run number  3
  
 This run uses the full linear solver DENSLV.
 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.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