pcprb2.out.txt

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

  
 PCPRB2:
   PITCON sample program.
   The aircraft stability problem
  
   Number of equations is  7
   Number of variables is  8
  
 using option  1
 fix variable  6  at  -5.00000000000000028E-2
 fix variable  8  at  0.E+0
 rheinboldt version of barray
  
 step  type of point     
  
   0  start point     0.000  0.000  0.000  0.000  0.000 -0.050  0.000  0.000
  
 PITCON 7.0
   University of Pittsburgh continuation code
  
   Modified:        12 November 1999
   Linear algebra:  LAPACK
   Precision:       double
  
   1  corrected       0.001  0.051  0.000  0.060  0.000 -0.050  0.000  0.000
 TRYSTP - DEBUG - about to return.
   2  continuation    0.100  0.051  0.006  0.060  0.001 -0.050 -0.010  0.000
 TRYSTP - DEBUG - about to return.
   3  continuation    0.241  0.052  0.014  0.060  0.002 -0.050 -0.025  0.000
 TRYSTP - DEBUG - about to return.
   4  continuation    0.466  0.054  0.027  0.060  0.005 -0.050 -0.049  0.000
 TRYSTP - DEBUG - about to return.
   5  continuation    0.855  0.061  0.051  0.061  0.009 -0.050 -0.091  0.000
 TRYSTP - DEBUG - about to return.
   6  continuation    1.237  0.074  0.075  0.063  0.016 -0.050 -0.135  0.000
 TRYSTP - DEBUG - about to return.
   7  continuation    1.609  0.096  0.100  0.065  0.025 -0.050 -0.182  0.000
 TRYSTP - DEBUG - about to return.
   8  continuation    1.966  0.136  0.125  0.067  0.039 -0.050 -0.237  0.000
 TRYSTP - DEBUG - about to return.
   9  continuation    2.303  0.207  0.146  0.068  0.064 -0.050 -0.303  0.000
 TRYSTP - DEBUG - about to return.
  10  continuation    2.572  0.315  0.155  0.067  0.100 -0.050 -0.374  0.000
 TRYSTP - DEBUG - about to return.
  11  continuation    2.801  0.499  0.140  0.060  0.160 -0.050 -0.454  0.000
 TRYSTP - DEBUG - about to return.
  12  continuation    2.954  0.786  0.083  0.044  0.254 -0.050 -0.504  0.000
 TRYSTP - DEBUG - about to return.
  13  continuation    3.008  1.141 -0.010  0.020  0.376 -0.050 -0.462  0.000
  14  limit point     2.965  0.826  0.074  0.041  0.267 -0.050 -0.505  0.000
 TRYSTP - DEBUG - about to return.
  15  continuation    2.997  1.485 -0.108 -0.004  0.500 -0.050 -0.327  0.000
 TRYSTP - DEBUG - about to return.
  16  continuation    2.965  1.764 -0.187 -0.023  0.606 -0.050 -0.160  0.000
 TRYSTP - DEBUG - about to return.
  17  continuation    2.922  2.036 -0.262 -0.042  0.714 -0.050  0.050  0.000
 TRYSTP - DEBUG - about to return.
  18  continuation    2.874  2.297 -0.332 -0.059  0.823 -0.050  0.290  0.000
 TRYSTP - DEBUG - about to return.
  19  continuation    2.825  2.546 -0.395 -0.074  0.931 -0.050  0.553  0.000
 TRYSTP - DEBUG - about to return.
  20  continuation    2.775  2.782 -0.454 -0.088  1.037 -0.050  0.834  0.000
 TRYSTP - DEBUG - about to return.
  21  continuation    2.727  3.007 -0.509 -0.101  1.143 -0.050  1.128  0.000
 TRYSTP - DEBUG - about to return.
  22  continuation    2.681  3.222 -0.559 -0.114  1.247 -0.050  1.434  0.000
 TRYSTP - DEBUG - about to return.
  23  continuation    2.636  3.427 -0.606 -0.125  1.351 -0.050  1.748  0.000
 TRYSTP - DEBUG - about to return.
  24  continuation    2.593  3.624 -0.651 -0.135  1.453 -0.050  2.070  0.000
 TRYSTP - DEBUG - about to return.
  25  continuation    2.551  3.812 -0.692 -0.145  1.555 -0.050  2.399  0.000
 TRYSTP - DEBUG - about to return.
  26  continuation    2.511  3.993 -0.732 -0.155  1.655 -0.050  2.733  0.000
 TRYSTP - DEBUG - about to return.
  27  continuation    2.472  4.166 -0.770 -0.163  1.755 -0.050  3.072  0.000
 TRYSTP - DEBUG - about to return.
  28  continuation    2.435  4.334 -0.806 -0.172  1.854 -0.050  3.415  0.000
 TRYSTP - DEBUG - about to return.
  29  continuation    2.399  4.495 -0.840 -0.180  1.952 -0.050  3.762  0.000
 TRYSTP - DEBUG - about to return.
  30  continuation    2.364  4.650 -0.873 -0.187  2.050 -0.050  4.112  0.000
 TRYSTP - DEBUG - about to return.
  31  continuation    2.330  4.801 -0.905 -0.194  2.147 -0.050  4.465  0.000
 TRYSTP - DEBUG - about to return.
  32  continuation    2.298  4.946 -0.936 -0.201  2.243 -0.050  4.821  0.000
 TRYSTP - DEBUG - about to return.
  33  continuation    2.266  5.086 -0.965 -0.207  2.340 -0.050  5.178  0.000
 TRYSTP - DEBUG - about to return.
  34  continuation    2.236  5.222 -0.994 -0.214  2.435 -0.050  5.539  0.000
 TRYSTP - DEBUG - about to return.
  35  continuation    2.206  5.354 -1.022 -0.220  2.531 -0.050  5.901  0.000
 TRYSTP - DEBUG - about to return.
  36  continuation    2.177  5.482 -1.050 -0.225  2.625 -0.050  6.264  0.000
 TRYSTP - DEBUG - about to return.
  37  continuation    2.149  5.606 -1.076 -0.231  2.720 -0.050  6.630  0.000
 TRYSTP - DEBUG - about to return.
  38  continuation    2.121  5.726 -1.102 -0.236  2.814 -0.050  6.997  0.000
 TRYSTP - DEBUG - about to return.
  39  continuation    2.095  5.842 -1.128 -0.241  2.908 -0.050  7.365  0.000
 TRYSTP - DEBUG - about to return.
  40  continuation    2.069  5.956 -1.152 -0.246  3.002 -0.050  7.735  0.000
 TRYSTP - DEBUG - about to return.
  41  continuation    2.043  6.066 -1.177 -0.251  3.095 -0.050  8.105  0.000
 TRYSTP - DEBUG - about to return.
  42  continuation    2.019  6.173 -1.201 -0.255  3.189 -0.050  8.477  0.000
 TRYSTP - DEBUG - about to return.
  43  continuation    1.994  6.277 -1.224 -0.260  3.282 -0.050  8.850  0.000
 TRYSTP - DEBUG - about to return.
  44  continuation    1.971  6.378 -1.247 -0.264  3.375 -0.050  9.223  0.000
 TRYSTP - DEBUG - about to return.
  45  continuation    1.948  6.477 -1.270 -0.268  3.467 -0.050  9.598  0.000
 TRYSTP - DEBUG - about to return.
  46  continuation    1.925  6.573 -1.292 -0.272  3.560 -0.050  9.973  0.000
 TRYSTP - DEBUG - about to return.
  47  continuation    1.903  6.666 -1.314 -0.276  3.652 -0.050 10.349  0.000
 TRYSTP - DEBUG - about to return.
  48  continuation    1.882  6.757 -1.336 -0.280  3.744 -0.050 10.726  0.000
 TRYSTP - DEBUG - about to return.
  49  continuation    1.860  6.846 -1.358 -0.283  3.836 -0.050 11.104  0.000
 TRYSTP - DEBUG - about to return.
  50  continuation    1.840  6.932 -1.379 -0.287  3.928 -0.050 11.482  0.000
  
 jacobians:      153
 factorizations: 153
 solves:         153
 functions:      154