pcprb4.out.txt

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

  
 PCPRB4:
   PITCON test problem
   Freudenstein-Roth function
  
   Number of equations is  2
   Number of variables is  3
  
 This is run number  1
 PITCON is free to choose parameterization.
 No limit points are sought.
  
 Step  Type of point     X(1)         X(2)         X(3)
  
   0  Start point      15.0000      -2.00000       0.00000    
  
 PITCON 7.0
   University of Pittsburgh continuation code
  
   Modified:        12 November 1999
   Linear algebra:  LAPACK
   Precision:       double
  
   1  Corrected        15.0000      -2.00000       0.00000    
 TRYSTP - DEBUG - about to return.
   2  Continuation     14.7105      -1.94205      0.653814E-01
 TRYSTP - DEBUG - about to return.
   3  Continuation     14.2846      -1.72915      0.268744    
 TRYSTP - DEBUG - about to return.
   4  Continuation     16.9061      -1.20942      0.546844    
 TRYSTP - DEBUG - about to return.
   5  Continuation     24.9177     -0.599075      0.555139    
 TRYSTP - DEBUG - about to return.
   6  Continuation     34.8979     -0.405252E-01  0.353317    
 TRYSTP - DEBUG - about to return.
   7  Continuation     44.8807      0.487801      0.594475E-01
 TRYSTP - DEBUG - about to return.
   8  Continuation     54.8605       1.09475     -0.304452    
 TRYSTP - DEBUG - about to return.
   9  Continuation     61.3809       1.81337     -0.624493    
 TRYSTP - DEBUG - about to return.
  10  Continuation    -48.5825       4.67841       2.87942    
  11  Target point     5.00000       4.00000       1.00000    
  
 We have reached the point we wanted.
 The code may stop now.
  
 This is run number  2
 The user fixes the parameterization.
 No limit points are sought.
  
 Step  Type of point     X(1)         X(2)         X(3)
  
   0  Start point      15.0000      -2.00000       0.00000    
  
 PITCON 7.0
   University of Pittsburgh continuation code
  
   Modified:        12 November 1999
   Linear algebra:  LAPACK
   Precision:       double
  
   1  Corrected        15.0000      -2.00000       0.00000    
 TRYSTP - DEBUG - about to return.
   2  Continuation     14.7407      -1.94890      0.578825E-01
 TRYSTP - DEBUG - about to return.
   3  Continuation     14.2963      -1.77723      0.227735    
 TRYSTP - DEBUG - about to return.
   4  Continuation     15.8338      -1.33743      0.504570    
 TRYSTP - DEBUG - about to return.
   5  Continuation     18.6039      -1.04741      0.578437    
 TRYSTP - DEBUG - about to return.
   6  Continuation     29.5771     -0.327668      0.476341    
 TRYSTP - DEBUG - about to return.
   7  Continuation     40.0071      0.227762      0.211791    
 TRYSTP - DEBUG - about to return.
   8  Continuation     49.5624      0.752615     -0.101854    
 TRYSTP - DEBUG - about to return.
   9  Continuation     57.9488       1.34333     -0.437080    
 TRYSTP - DEBUG - about to return.
  10  Continuation     60.8565       2.25880     -0.686030    
 TRYSTP - DEBUG - about to return.
  11  Continuation     55.6686       2.70414     -0.589661    
 TRYSTP - DEBUG - about to return.
  12  Continuation     40.7916       3.27084     -0.169093    
 TRYSTP - DEBUG - about to return.
  13  Continuation     29.3855       3.55244      0.189747    
 TRYSTP - DEBUG - about to return.
  14  Continuation     18.4706       3.77123      0.546943    
 TRYSTP - DEBUG - about to return.
  15  Continuation     7.78375       3.95562      0.905475    
 TRYSTP - DEBUG - about to return.
  16  Continuation    -2.76866       4.11736       1.26589    
  17  Target point     5.00000       4.00000       1.00000    
  
 We have reached the point we wanted.
 The code may stop now.
  
 This is run number  3
 PITCON is free to choose parameterization.
 Seek limit points in index  1
  
 Step  Type of point     X(1)         X(2)         X(3)
  
   0  Start point      15.0000      -2.00000       0.00000    
  
 PITCON 7.0
   University of Pittsburgh continuation code
  
   Modified:        12 November 1999
   Linear algebra:  LAPACK
   Precision:       double
  
   1  Corrected        15.0000      -2.00000       0.00000    
 TRYSTP - DEBUG - about to return.
   2  Continuation     14.7105      -1.94205      0.653814E-01
 TRYSTP - DEBUG - about to return.
   3  Continuation     14.2846      -1.72915      0.268744    
   4  Limit point      14.2831      -1.74138      0.258578    
 TRYSTP - DEBUG - about to return.
   5  Continuation     16.9061      -1.20942      0.546844    
 TRYSTP - DEBUG - about to return.
   6  Continuation     24.9177     -0.599075      0.555139    
 TRYSTP - DEBUG - about to return.
   7  Continuation     34.8979     -0.405252E-01  0.353317    
 TRYSTP - DEBUG - about to return.
   8  Continuation     44.8807      0.487801      0.594475E-01
 TRYSTP - DEBUG - about to return.
   9  Continuation     54.8605       1.09475     -0.304452    
 TRYSTP - DEBUG - about to return.
  10  Continuation     61.3809       1.81337     -0.624493    
 TRYSTP - DEBUG - about to return.
  11  Continuation    -48.5825       4.67841       2.87942    
  12  Target point     5.00000       4.00000       1.00000    
  
 We have reached the point we wanted.
 The code may stop now.
  
 This is run number  4
 The user fixes the parameterization.
 Seek limit points in index  1
  
 Step  Type of point     X(1)         X(2)         X(3)
  
   0  Start point      15.0000      -2.00000       0.00000    
  
 PITCON 7.0
   University of Pittsburgh continuation code
  
   Modified:        12 November 1999
   Linear algebra:  LAPACK
   Precision:       double
  
   1  Corrected        15.0000      -2.00000       0.00000    
 TRYSTP - DEBUG - about to return.
   2  Continuation     14.7407      -1.94890      0.578825E-01
 TRYSTP - DEBUG - about to return.
   3  Continuation     14.2963      -1.77723      0.227735    
 TRYSTP - DEBUG - about to return.
   4  Continuation     15.8338      -1.33743      0.504570    
   5  Limit point      14.2831      -1.74138      0.258578    
 TRYSTP - DEBUG - about to return.
   6  Continuation     18.6039      -1.04741      0.578437    
 TRYSTP - DEBUG - about to return.
   7  Continuation     29.5771     -0.327668      0.476341    
 TRYSTP - DEBUG - about to return.
   8  Continuation     40.0071      0.227762      0.211791    
 TRYSTP - DEBUG - about to return.
   9  Continuation     49.5624      0.752615     -0.101854    
 TRYSTP - DEBUG - about to return.
  10  Continuation     57.9488       1.34333     -0.437080    
 TRYSTP - DEBUG - about to return.
  11  Continuation     60.8565       2.25880     -0.686030    
  12  Limit point      61.6694       1.98380     -0.663880    
 TRYSTP - DEBUG - about to return.
  13  Continuation     55.6686       2.70414     -0.589661    
 TRYSTP - DEBUG - about to return.
  14  Continuation     40.7916       3.27084     -0.169093    
 TRYSTP - DEBUG - about to return.
  15  Continuation     29.3855       3.55244      0.189747    
 TRYSTP - DEBUG - about to return.
  16  Continuation     18.4706       3.77123      0.546943    
 TRYSTP - DEBUG - about to return.
  17  Continuation     7.78375       3.95562      0.905475    
 TRYSTP - DEBUG - about to return.
  18  Continuation    -2.76866       4.11736       1.26589    
  19  Target point     5.00000       4.00000       1.00000    
  
 We have reached the point we wanted.
 The code may stop now.
  
 This is run number  5
 PITCON is free to choose parameterization.
 Seek limit points in index  3
  
 Step  Type of point     X(1)         X(2)         X(3)
  
   0  Start point      15.0000      -2.00000       0.00000    
  
 PITCON 7.0
   University of Pittsburgh continuation code
  
   Modified:        12 November 1999
   Linear algebra:  LAPACK
   Precision:       double
  
   1  Corrected        15.0000      -2.00000       0.00000    
 TRYSTP - DEBUG - about to return.
   2  Continuation     14.7105      -1.94205      0.653814E-01
 TRYSTP - DEBUG - about to return.
   3  Continuation     14.2846      -1.72915      0.268744    
 TRYSTP - DEBUG - about to return.
   4  Continuation     16.9061      -1.20942      0.546844    
 TRYSTP - DEBUG - about to return.
   5  Continuation     24.9177     -0.599075      0.555139    
   6  Limit point      20.4860     -0.896795      0.587587    
 TRYSTP - DEBUG - about to return.
   7  Continuation     34.8979     -0.405252E-01  0.353317    
 TRYSTP - DEBUG - about to return.
   8  Continuation     44.8807      0.487801      0.594475E-01
 TRYSTP - DEBUG - about to return.
   9  Continuation     54.8605       1.09475     -0.304452    
 TRYSTP - DEBUG - about to return.
  10  Continuation     61.3809       1.81337     -0.624493    
 TRYSTP - DEBUG - about to return.
  11  Continuation    -48.5825       4.67841       2.87942    
  12  Target point     5.00000       4.00000       1.00000    
  
 We have reached the point we wanted.
 The code may stop now.
  
 This is run number  6
 The user fixes the parameterization.
 Seek limit points in index  3
  
 Step  Type of point     X(1)         X(2)         X(3)
  
   0  Start point      15.0000      -2.00000       0.00000    
  
 PITCON 7.0
   University of Pittsburgh continuation code
  
   Modified:        12 November 1999
   Linear algebra:  LAPACK
   Precision:       double
  
   1  Corrected        15.0000      -2.00000       0.00000    
 TRYSTP - DEBUG - about to return.
   2  Continuation     14.7407      -1.94890      0.578825E-01
 TRYSTP - DEBUG - about to return.
   3  Continuation     14.2963      -1.77723      0.227735    
 TRYSTP - DEBUG - about to return.
   4  Continuation     15.8338      -1.33743      0.504570    
 TRYSTP - DEBUG - about to return.
   5  Continuation     18.6039      -1.04741      0.578437    
 TRYSTP - DEBUG - about to return.
   6  Continuation     29.5771     -0.327668      0.476341    
   7  Limit point      20.4859     -0.896799      0.587587    
 TRYSTP - DEBUG - about to return.
   8  Continuation     40.0071      0.227762      0.211791    
 TRYSTP - DEBUG - about to return.
   9  Continuation     49.5624      0.752615     -0.101854    
 TRYSTP - DEBUG - about to return.
  10  Continuation     57.9488       1.34333     -0.437080    
 TRYSTP - DEBUG - about to return.
  11  Continuation     60.8565       2.25880     -0.686030    
  12  Limit point      61.0203       2.23014     -0.686353    
 TRYSTP - DEBUG - about to return.
  13  Continuation     55.6686       2.70414     -0.589661    
 TRYSTP - DEBUG - about to return.
  14  Continuation     40.7916       3.27084     -0.169093    
 TRYSTP - DEBUG - about to return.
  15  Continuation     29.3855       3.55244      0.189747    
 TRYSTP - DEBUG - about to return.
  16  Continuation     18.4706       3.77123      0.546943    
 TRYSTP - DEBUG - about to return.
  17  Continuation     7.78375       3.95562      0.905475    
 TRYSTP - DEBUG - about to return.
  18  Continuation    -2.76866       4.11736       1.26589    
  19  Target point     5.00000       4.00000       1.00000    
  
 We have reached the point we wanted.
 The code may stop now.
  
 All tests completed.