pcprb1.out.txt

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

  
 PCPRB1:
   PITCON test problem
   Freudenstein-Roth function
  
   Number of equations is  2
   Number of variables is  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
  
 CHECKW - Machine epsilon =  2.22044604925031308E-16
 CHECKW - Jacobian finite difference increment =  1.220703125E-4
 START  - Checking the initial point.
   Fixing variable number 2
  
 Step      FX             X             DX
  
    0   0.00000       15.1327    
              1             1
 COQUAL - Corrector convergence quality factor =  8.
  1  Corrected        15.0000      -2.00000       0.00000    
 TANPAR - Continuation index: First choice =  1
                              Second choice =  3
 TANPAR - Tangent vector has limit component =  -0.96514602486622125
 TRYSTP - Predictor using stepsize  0.30000001192092896
 TRYSTP - Corrector is fixing index  1
  
 Step      FX             X             DX
  
    0  0.313698E-01   14.8391    
              1             1
                                   0.673116E-02
                                          2
    1  0.537876E-03   14.8383    
              1             1
                                   0.118851E-03
                                          2
    2  0.167436E-06   14.8382    
              1             1
 TRYSTP - DEBUG - about to return.
 COQUAL - Corrector convergence quality factor =  8.
  2  Continuation     14.7105      -1.94205      0.653814E-01
 TANPAR - A limit point may be coming
 in the preferable index  1
 We'll try second-best index.
 TANPAR - Continuation index: First choice =  3
                              Second choice =  1
 TANPAR - Tangent vector has limit component =  -0.94636179587137015
 SETSTP - Next stepsize HTAN =  0.85405264522127078
 TRYSTP - Predictor using stepsize  0.85405264522127078
 TRYSTP - Corrector is fixing index  3
  
 Step      FX             X             DX
  
    0  0.403853       14.0152    
              1             1
                                   0.375401    
                                          1
    1  0.742820E-02   14.3845    
              1             1
                                   0.700336E-02
                                          1
    2  0.265002E-05   14.3914    
              1             1
                                   0.249236E-05
                                          1
    3  0.341413E-12   14.3914    
              1             1
 TRYSTP - DEBUG - about to return.
 COQUAL - Corrector convergence quality factor =  8.
  3  Continuation     14.2846      -1.72915      0.268744    
 TANPAR - Continuation index: First choice =  2
                              Second choice =  3
 TANPAR - Tangent vector has limit component =  0.18928168682171298
 LIMIT  - Attempt correction of approximate limit point.
 LIMIT - For S =  0.E+0  Tan(X(S))(Lim)= -0.94636179587137015
 LIMIT - For S =  1.  Tan(X(S))(Lim)= 0.18928168682171298
  
 Step      FX             X             DX
  
    0  0.713931E-01   14.4656    
              1             1
                                   0.669412E-01
                                          1
    1  0.716072E-14   14.3992    
              2             1
                                   0.732747E-14
                                          1
    2  0.874754E-14   14.3992    
              2             1
 LIMIT - For S =  0.83332648872087267  Tan(X(S))(Lim)= -0.34063910414628051
  
 Step      FX             X             DX
  
    0  0.320282E-02   14.3942    
              1             1
                                   0.296802E-02
                                          1
    1   0.00000       14.3912    
              1             1
 LIMIT - For S =  0.94046611512301181  Tan(X(S))(Lim)= -6.98064563124521816E-3
  
 Step      FX             X             DX
  
    0  0.608814E-04   14.3912    
              1             1
                                   0.563383E-04
                                          1
    1  0.397205E-14   14.3912    
              1             1
 LIMIT - For S =  0.94258361239405442  Tan(X(S))(Lim)= 9.59027887242148232E-5
  
 Step      FX             X             DX
  
    0  0.138854E-04   14.3912    
              1             1
                                   0.130114E-04
                                          1
    1  0.502430E-14   14.3912    
              1             1
 LIMIT - For S =  0.94255491565076133  Tan(X(S))(Lim)= -1.61440507147034498E-8
  4  Limit point      14.2831      -1.74138      0.258578    
 SETSTP - Next stepsize HTAN =  0.68579738650057587
 TRYSTP - Predictor using stepsize  0.68579738650057587
 TRYSTP - Corrector is fixing index  2
  
 Step      FX             X             DX
  
    0   2.79306       14.4819    
              1             1
                                    2.49163    
                                          1
    1  0.888178E-15   16.9581    
              2             1
 TRYSTP - DEBUG - about to return.
 COQUAL - Corrector convergence quality factor =  8.
  5  Continuation     16.9061      -1.20942      0.546844    
 TANPAR - Continuation index: First choice =  1
                              Second choice =  2
 TANPAR - Tangent vector has limit component =  0.99391377975399309
 SETSTP - Next stepsize HTAN =  8.0606778005602493
 TRYSTP - Predictor using stepsize  8.0606778005602493
 TRYSTP - Corrector is fixing index  1
  
 Step      FX             X             DX
  
    0   5.86072       24.9323    
              1             1
                                   0.236730    
                                          2
    1  0.353014       24.9311    
              1             1
                                   0.131102E-01
                                          3
    2  0.742281E-03   24.9311    
              1             1
                                   0.287064E-04
                                          3
    3  0.473759E-08   24.9311    
              1             1
 TRYSTP - DEBUG - about to return.
 COQUAL - Corrector convergence quality factor =  8.
  6  Continuation     24.9177     -0.599075      0.555139    
 TANPAR - Continuation index: First choice =  1
                              Second choice =  2
 TANPAR - Tangent vector has limit component =  0.99802289059366012
 SETSTP - Next stepsize HTAN =  10.
 TRYSTP - Predictor using stepsize  10.
 TRYSTP - Corrector is fixing index  1
  
 Step      FX             X             DX
  
    0   2.33947       34.9005    
              1             1
                                   0.718304E-01
                                          3
    1  0.162190E-01   34.8997    
              1             1
                                   0.476056E-03
                                          3
    2  0.680609E-07   34.8997    
              1             1
                                   0.201656E-08
                                          3
    3  0.502430E-14   34.8997    
              1             1
 TRYSTP - DEBUG - about to return.
 COQUAL - Corrector convergence quality factor =  8.
  7  Continuation     34.8979     -0.405252E-01  0.353317    
 TANPAR - Continuation index: First choice =  1
                              Second choice =  2
 TANPAR - Tangent vector has limit component =  0.99828072813109192
 SETSTP - Next stepsize HTAN =  10.
 TRYSTP - Predictor using stepsize  10.
 TRYSTP - Corrector is fixing index  1
  
 Step      FX             X             DX
  
    0   1.33519       44.8835    
              1             1
                                   0.376262E-01
                                          3
    1  0.584960E-05   44.8834    
              1             1
                                   0.150724E-06
                                          3
    2  0.966443E-13   44.8834    
              1             1
 TRYSTP - DEBUG - about to return.
 COQUAL - Corrector convergence quality factor =  8.
  8  Continuation     44.8807      0.487801      0.594475E-01
 TANPAR - Continuation index: First choice =  1
                              Second choice =  2
 TANPAR - Tangent vector has limit component =  0.99797927558197019
 SETSTP - Next stepsize HTAN =  10.
 TRYSTP - Predictor using stepsize  10.
 TRYSTP - Corrector is fixing index  1
  
 Step      FX             X             DX
  
    0   1.25604       54.8709    
              2             1
                                   0.614508E-01
                                          2
    1  0.171716E-01   54.8722    
              2             1
                                   0.139305E-02
                                          2
    2  0.895293E-05   54.8723    
              2             1
                                   0.748551E-06
                                          2
    3  0.258611E-11   54.8723    
              2             1
 TRYSTP - DEBUG - about to return.
 COQUAL - Corrector convergence quality factor =  8.
  9  Continuation     54.8605       1.09475     -0.304452    
 TANPAR - Continuation index: First choice =  1
                              Second choice =  2
 TANPAR - Tangent vector has limit component =  0.99658696009041658
 SETSTP - Next stepsize HTAN =  10.
 TRYSTP - Predictor using stepsize  10.
 TRYSTP - Corrector is fixing index  1
  
 Step      FX             X             DX
  
    0   2.63438       64.8556    
              2             1
                                    1.03403    
                                          2
    1   8.14547       64.8954    
              2             1
 CORRECTOR - Warning!
   The residual FX is not decreasing.
 TRYSTP - Predictor using stepsize  10.
 TRYSTP - Corrector is fixing index  2
  
 Step      FX             X             DX
  
    0   2.63438       64.8556    
              2             1
                                    3.44549    
                                          1
    1  0.355271E-14   61.4108    
              2             1
                                   0.503301E-14
                                          1
    2  0.100486E-13   61.4108    
              1             1
 TRYSTP - DEBUG - about to return.
 COQUAL - Corrector convergence quality factor =  8.
 10  Continuation     61.3809       1.81337     -0.624493    
 TANPAR - A limit point may be coming
 in the preferable index  1
 We'll try second-best index.
 TANPAR - Continuation index: First choice =  2
                              Second choice =  1
 TANPAR - Tangent vector has limit component =  0.95465732033638262
 SETSTP - Next stepsize HTAN =  10.
 TRYSTP - Predictor using stepsize  10.
 TRYSTP - Corrector is fixing index  2
  
 Step      FX             X             DX
  
    0   81.0557       71.0960    
              2             1
                                    119.510    
                                          1
    1  0.142109E-13   48.8921    
              2             1
                                   0.201320E-13
                                          1
    2  0.710543E-14   48.8921    
              1             1
 TRYSTP - DEBUG - about to return.
 COQUAL - Corrector convergence quality factor =  8.
 11  Continuation    -48.5825       4.67841       2.87942    
  
 TARGET - A target point has been detected.
  
 Step      FX             X             DX
  
    0   23.6765       10.9095    
              2             1
                                    9.18494    
                                          1
    1   26.2709       4.81676    
              2             2
                                    2.97641    
                                          1
    2   3.51787       5.94022    
              2             1
                                   0.784750    
                                          1
    3  0.112048       6.46378    
              2             1
                                   0.243219E-01
                                          1
    4  0.125657E-03   6.48072    
              2             1
                                   0.268572E-04
                                          1
    5  0.158203E-09   6.48074    
              2             1
 12  Target point     5.00000       4.00000       1.00000    
  
 PITCON reached the target point.
  
 The computation succeeded.