test_simeq_newton4.adb running test_case_0 A( 1, 1)= 4.163100E-02 A( 1, 2)= 1.766430E-01 A( 1, 3)= 3.646020E-01 A( 1, 4)= 9.133100E-02 A( 2, 1)= 9.229800E-02 A( 2, 2)= 4.872170E-01 A( 2, 3)= 5.267500E-01 A( 2, 4)= 4.544330E-01 X_soln( 1)= 1.100000E+00 X_soln( 2)= 2.300000E+00 Y( 1)= 3.541703E+00 Y( 2)= 1.471024E+01 X( 1)= 1.000000E+00 X( 2)= 1.000000E+00 terms: x1 x2 power terms: x1 * x1 * x1 * x1 x2 * x2 * x2 * x2 simeq_newton4.adb running var1=( 1, 2, 1, 2) var2=( 0, 0, 1, 2) var3=( 0, 0, 1, 2) var4=( 0, 0, 1, 2) residual X_resid( 1)= -2.867496E+00 residual X_resid( 2)= -1.314954E+01 itr 1, prev= 1.00000000000000E+12, residual= 1.60170377156000E+01 X_next( 1)= 8.287889E-01 X_next( 2)= 5.442051E+00 X_next( 3)= 4.718194E-01 X_next( 4)= 8.771026E+02 residual X_resid( 1)= 7.773278E+01 residual X_resid( 2)= 3.868506E+02 itr 2, prev= 1.60170377156000E+01, residual= 4.64583384350258E+02 X_next( 1)= 9.019763E-01 X_next( 2)= 4.947378E+00 X_next( 3)= 6.618818E-01 X_next( 4)= 5.991014E+02 residual X_resid( 1)= 5.232762E+01 residual X_resid( 2)= 2.603835E+02 itr 3, prev= 4.64583384350258E+02, residual= 3.12711170375205E+02 X_next( 1)= 9.869449E-01 X_next( 2)= 4.135280E+00 X_next( 3)= 9.487934E-01 X_next( 4)= 2.924284E+02 residual X_resid( 1)= 2.428356E+01 residual X_resid( 2)= 1.207845E+02 itr 4, prev= 3.12711170375205E+02, residual= 1.45068092483651E+02 X_next( 1)= 1.041013E+00 X_next( 2)= 3.344555E+00 X_next( 3)= 1.174421E+00 X_next( 4)= 1.251277E+02 residual X_resid( 1)= 8.948661E+00 residual X_resid( 2)= 4.449615E+01 itr 5, prev= 1.45068092483651E+02, residual= 5.34448067509087E+01 X_next( 1)= 1.076111E+00 X_next( 2)= 2.744397E+00 X_next( 3)= 1.340998E+00 X_next( 4)= 5.672672E+01 residual X_resid( 1)= 2.657714E+00 residual X_resid( 2)= 1.321106E+01 itr 6, prev= 5.34448067509087E+01, residual= 1.58687762412134E+01 X_next( 1)= 1.094301E+00 X_next( 2)= 2.409543E+00 X_next( 3)= 1.433994E+00 X_next( 4)= 3.370846E+01 residual X_resid( 1)= 5.309478E-01 residual X_resid( 2)= 2.638326E+00 itr 7, prev= 1.58687762412134E+01, residual= 3.16927428309605E+00 X_next( 1)= 1.099531E+00 X_next( 2)= 2.309155E+00 X_next( 3)= 1.461607E+00 X_next( 4)= 2.843231E+01 residual X_resid( 1)= 4.162394E-02 residual X_resid( 2)= 2.067849E-01 itr 8, prev= 3.16927428309605E+00, residual= 2.48408838689476E-01 X_next( 1)= 1.099993E+00 X_next( 2)= 2.300142E+00 X_next( 3)= 1.464061E+00 X_next( 4)= 2.799103E+01 residual X_resid( 1)= 6.437583E-04 residual X_resid( 2)= 3.197996E-03 itr 9, prev= 2.48408838689476E-01, residual= 3.84175404726506E-03 X_next( 1)= 1.100000E+00 X_next( 2)= 2.300001E+00 X_next( 3)= 1.464100E+00 X_next( 4)= 2.798413E+01 residual X_resid( 1)= 3.201694E-06 residual X_resid( 2)= 1.590502E-05 itr 10, prev= 3.84175404726506E-03, residual= 1.91067102415232E-05 X_next( 1)= 1.100000E+00 X_next( 2)= 2.300000E+00 X_next( 3)= 1.464100E+00 X_next( 4)= 2.798410E+01 residual X_resid( 1)= 7.818085E-09 residual X_resid( 2)= 3.883781E-08 itr 11, prev= 1.91067102415232E-05, residual= 4.66558964795638E-08 simeq_newton4 found solution test case0 solution X( 1)= 1.09999999991174E+00 X( 2)= 2.30000000172953E+00 sum of residuals= 4.66558964795638E-08 end test case0 test_case_1 A( 1, 1)= 7.910778E-01 A( 1, 2)= 6.442950E-01 A( 1, 3)= 6.658510E-01 A( 1, 4)= 9.577035E-01 A( 2, 1)= 1.223651E-01 A( 2, 2)= 5.898427E-01 A( 2, 3)= 4.864346E-01 A( 2, 4)= 5.060496E-01 X_soln( 1)= 1.100000E+00 X_soln( 2)= 2.300000E+00 Y( 1)= 3.003878E+01 Y( 2)= 1.630003E+01 X( 1)= 1.000000E+00 X( 2)= 1.000000E+00 terms: x1 x2 power terms: x1 * x1 * x1 x2 * x2 * x2 * x2 simeq_newton4.adb running var1=( 1, 2, 1, 2) var2=( 0, 0, 1, 2) var3=( 0, 0, 1, 2) var4=( 0, 0, 0, 2) residual X_resid( 1)= -2.697985E+01 residual X_resid( 2)= -1.459533E+01 itr 1, prev= 1.00000000000000E+12, residual= 4.15751884197385E+01 X_next( 1)= 1.948173E+01 X_next( 2)= -5.995099E+00 X_next( 3)= 7.394052E+03 X_next( 4)= 1.291770E+03 residual X_resid( 1)= 6.141980E+03 residual X_resid( 2)= 4.232970E+03 itr 2, prev= 4.15751884197385E+01, residual= 1.03749505251697E+04 X_next( 1)= 1.705972E+01 X_next( 2)= -5.431208E+00 X_next( 3)= 4.964960E+03 X_next( 4)= 8.701333E+02 residual X_resid( 1)= 4.119211E+03 residual X_resid( 2)= 2.838043E+03 itr 3, prev= 1.03749505251697E+04, residual= 6.95725366166110E+03 X_next( 1)= 1.317839E+01 X_next( 2)= -4.497495E+00 X_next( 3)= 2.288689E+03 X_next( 4)= 4.091503E+02 residual X_resid( 1)= 1.893259E+03 residual X_resid( 2)= 1.303007E+03 itr 4, prev= 6.95725366166110E+03, residual= 3.19626635326801E+03 X_next( 1)= 9.518731E+00 X_next( 2)= -3.563335E+00 X_next( 3)= 8.624563E+02 X_next( 4)= 1.612228E+02 residual X_resid( 1)= 7.038664E+02 residual X_resid( 2)= 4.838782E+02 itr 5, prev= 3.19626635326801E+03, residual= 1.18774467398540E+03 X_next( 1)= 6.669585E+00 X_next( 2)= -2.804975E+00 X_next( 3)= 2.966856E+02 X_next( 4)= 6.190361E+01 residual X_resid( 1)= 2.302638E+02 residual X_resid( 2)= 1.585060E+02 itr 6, prev= 1.18774467398540E+03, residual= 3.88769838149359E+02 X_next( 1)= 4.659137E+00 X_next( 2)= -2.315879E+00 X_next( 3)= 1.011385E+02 X_next( 4)= 2.876493E+01 residual X_resid( 1)= 6.704627E+01 residual X_resid( 2)= 4.665781E+01 itr 7, prev= 3.88769838149359E+02, residual= 1.13704079214585E+02 X_next( 1)= 3.353991E+00 X_next( 2)= -2.143082E+00 X_next( 3)= 3.772991E+01 X_next( 4)= 2.109382E+01 residual X_resid( 1)= 1.655784E+01 residual X_resid( 2)= 1.187396E+01 itr 8, prev= 1.13704079214585E+02, residual= 2.84317944617749E+01 X_next( 1)= 2.626217E+00 X_next( 2)= -2.157474E+00 X_next( 3)= 1.811306E+01 X_next( 4)= 2.166619E+01 residual X_resid( 1)= 3.459094E+00 residual X_resid( 2)= 2.523745E+00 itr 9, prev= 2.84317944617749E+01, residual= 5.98283948125383E+00 X_next( 1)= 2.349565E+00 X_next( 2)= -2.173022E+00 X_next( 3)= 1.297068E+01 X_next( 4)= 2.229753E+01 residual X_resid( 1)= 4.107978E-01 residual X_resid( 2)= 2.987772E-01 itr 10, prev= 5.98283948125383E+00, residual= 7.09575096300053E-01 X_next( 1)= 2.307691E+00 X_next( 2)= -2.175223E+00 X_next( 3)= 1.228947E+01 X_next( 4)= 2.238800E+01 residual X_resid( 1)= 9.315017E-03 residual X_resid( 2)= 6.775216E-03 itr 11, prev= 7.09575096300053E-01, residual= 1.60902327313046E-02 X_next( 1)= 2.306701E+00 X_next( 2)= -2.175275E+00 X_next( 3)= 1.227365E+01 X_next( 4)= 2.239014E+01 residual X_resid( 1)= 1.597041E-05 residual X_resid( 2)= 1.161577E-05 itr 12, prev= 1.60902327313046E-02, residual= 2.75861786604992E-05 X_next( 1)= 2.306699E+00 X_next( 2)= -2.175275E+00 X_next( 3)= 1.227362E+01 X_next( 4)= 2.239014E+01 simeq_newton4 did not converge < 1.00000000000000E-06 test case1 solution X( 1)= 2.30669888005092E+00 X( 2)= -2.17527528294027E+00 sum of residuals= 1.68606888450995E-08 end test case1 test_case_2 A( 1, 1)= 1.750984E-01 A( 1, 2)= 8.782248E-01 A( 1, 3)= 3.246958E-01 A( 1, 4)= 1.620386E-01 A( 1, 5)= 3.833710E-01 A( 1, 6)= 3.158238E-01 A( 1, 7)= 5.133576E-02 A( 2, 1)= 8.000549E-01 A( 2, 2)= 5.223173E-01 A( 2, 3)= 5.863344E-01 A( 2, 4)= 5.226468E-01 A( 2, 5)= 1.255568E-01 A( 2, 6)= 2.337342E-01 A( 2, 7)= 3.708149E-01 A( 3, 1)= 2.857049E-01 A( 3, 2)= 8.418347E-01 A( 3, 3)= 7.152572E-01 A( 3, 4)= 3.275851E-01 A( 3, 5)= 7.227268E-01 A( 3, 6)= 8.689656E-01 A( 3, 7)= 7.055294E-01 A( 4, 1)= 8.333382E-01 A( 4, 2)= 9.146897E-01 A( 4, 3)= 1.905640E-01 A( 4, 4)= 8.099251E-01 A( 4, 5)= 4.112297E-01 A( 4, 6)= 5.375608E-01 A( 4, 7)= 7.839017E-01 X_soln( 1)= 1.100000E+00 X_soln( 2)= 1.200000E+00 X_soln( 3)= 1.400000E+00 X_soln( 4)= 1.700000E+00 Y( 1)= 3.660337E+00 Y( 2)= 7.033351E+00 Y( 3)= 1.179022E+01 Y( 4)= 1.202289E+01 X( 1)= 1.000000E+00 X( 2)= 1.000000E+00 X( 3)= 1.000000E+00 X( 4)= 1.000000E+00 terms: x1 x2 x3 x4 power terms: x1 * x2 x1 * x1 * x3 * x3 x4 * x4 * x4 * x4 simeq_newton4.adb running var1=( 1, 2, 3, 4, 2, 3, 4) var2=( 0, 0, 0, 0, 1, 3, 4) var3=( 0, 0, 0, 0, 0, 1, 4) var4=( 0, 0, 0, 0, 0, 1, 4) residual X_resid( 1)= -1.369748E+00 residual X_resid( 2)= -3.871892E+00 residual X_resid( 3)= -7.322619E+00 residual X_resid( 4)= -7.541681E+00 itr 1, prev= 1.00000000000000E+12, residual= 2.01059396398578E+01 X_next( 1)= -4.164781E+00 X_next( 2)= 3.700338E+00 X_next( 3)= 3.931958E+00 X_next( 4)= 3.622100E+00 X_next( 5)= -1.541110E+01 X_next( 6)= 2.681651E+02 X_next( 7)= 1.721242E+02 residual X_resid( 1)= 8.834465E+01 residual X_resid( 2)= 1.203365E+02 residual X_resid( 3)= 3.374608E+02 residual X_resid( 4)= 2.643200E+02 itr 2, prev= 2.01059396398578E+01, residual= 8.10461949683850E+02 X_next( 1)= -3.849124E+00 X_next( 2)= 4.602595E+00 X_next( 3)= 2.933285E+00 X_next( 4)= 3.295893E+00 X_next( 5)= -1.771596E+01 X_next( 6)= 1.274771E+02 X_next( 7)= 1.180028E+02 residual X_resid( 1)= 4.072058E+01 residual X_resid( 2)= 6.706222E+01 residual X_resid( 3)= 1.753863E+02 residual X_resid( 4)= 1.459518E+02 itr 3, prev= 8.10461949683850E+02, residual= 4.29120899469680E+02 X_next( 1)= -3.597772E+00 X_next( 2)= 5.835093E+00 X_next( 3)= 1.829891E+00 X_next( 4)= 2.763649E+00 X_next( 5)= -2.099333E+01 X_next( 6)= 4.334289E+01 X_next( 7)= 5.833533E+01 residual X_resid( 1)= 1.051137E+01 residual X_resid( 2)= 2.477981E+01 residual X_resid( 3)= 5.795656E+01 residual X_resid( 4)= 5.329883E+01 itr 4, prev= 4.29120899469680E+02, residual= 1.46546577040454E+02 X_next( 1)= -3.559562E+00 X_next( 2)= 6.601499E+00 X_next( 3)= 1.381202E+00 X_next( 4)= 2.254325E+00 X_next( 5)= -2.349844E+01 X_next( 6)= 2.417172E+01 X_next( 7)= 2.582651E+01 residual X_resid( 1)= 2.278954E+00 residual X_resid( 2)= 7.831164E+00 residual X_resid( 3)= 1.671936E+01 residual X_resid( 4)= 1.671411E+01 itr 5, prev= 1.46546577040454E+02, residual= 4.35435926119187E+01 X_next( 1)= -3.539008E+00 X_next( 2)= 6.932928E+00 X_next( 3)= 1.267215E+00 X_next( 4)= 1.886273E+00 X_next( 5)= -2.453569E+01 X_next( 6)= 2.011239E+01 X_next( 7)= 1.265954E+01 residual X_resid( 1)= 1.213532E-01 residual X_resid( 2)= 1.799977E+00 residual X_resid( 3)= 3.235394E+00 residual X_resid( 4)= 3.784283E+00 itr 6, prev= 4.35435926119187E+01, residual= 8.94100768031901E+00 X_next( 1)= -3.513456E+00 X_next( 2)= 7.029175E+00 X_next( 3)= 1.255266E+00 X_next( 4)= 1.707432E+00 X_next( 5)= -2.469670E+01 X_next( 6)= 1.945094E+01 X_next( 7)= 8.499122E+00 residual X_resid( 1)= -3.067073E-01 residual X_resid( 2)= 5.265363E-02 residual X_resid( 3)= -3.698446E-01 residual X_resid( 4)= 6.335386E-02 itr 7, prev= 8.94100768031901E+00, residual= 7.92559371553069E-01 X_next( 1)= -3.496166E+00 X_next( 2)= 7.020577E+00 X_next( 3)= 1.267413E+00 X_next( 4)= 1.669440E+00 X_next( 5)= -2.454511E+01 X_next( 6)= 1.963452E+01 X_next( 7)= 7.767526E+00 residual X_resid( 1)= -2.349054E-01 residual X_resid( 2)= -1.600836E-01 residual X_resid( 3)= -6.229794E-01 residual X_resid( 4)= -3.710339E-01 itr 8, prev= 7.92559371553069E-01, residual= 1.38900235912510E+00 X_next( 1)= -3.491528E+00 X_next( 2)= 7.010049E+00 X_next( 3)= 1.272552E+00 X_next( 4)= 1.669135E+00 X_next( 5)= -2.447579E+01 X_next( 6)= 1.974159E+01 X_next( 7)= 7.761864E+00 residual X_resid( 1)= -1.816217E-01 residual X_resid( 2)= -1.273892E-01 residual X_resid( 3)= -4.878009E-01 residual X_resid( 4)= -2.944444E-01 itr 9, prev= 1.38900235912510E+00, residual= 1.09125625705411E+00 X_next( 1)= -3.486131E+00 X_next( 2)= 6.997474E+00 X_next( 3)= 1.278591E+00 X_next( 4)= 1.669193E+00 X_next( 5)= -2.439411E+01 X_next( 6)= 1.986784E+01 X_next( 7)= 7.762939E+00 residual X_resid( 1)= -1.185110E-01 residual X_resid( 2)= -8.590606E-02 residual X_resid( 3)= -3.230140E-01 residual X_resid( 4)= -1.979552E-01 itr 10, prev= 1.09125625705411E+00, residual= 7.25386265733362E-01 X_next( 1)= -3.482054E+00 X_next( 2)= 6.987682E+00 X_next( 3)= 1.283204E+00 X_next( 4)= 1.669612E+00 X_next( 5)= -2.433149E+01 X_next( 6)= 1.996469E+01 X_next( 7)= 7.770728E+00 residual X_resid( 1)= -6.983223E-02 residual X_resid( 2)= -5.144443E-02 residual X_resid( 3)= -1.917332E-01 residual X_resid( 4)= -1.183698E-01 itr 11, prev= 7.25386265733362E-01, residual= 4.31379643061316E-01 X_next( 1)= -3.479491E+00 X_next( 2)= 6.981435E+00 X_next( 3)= 1.286120E+00 X_next( 4)= 1.669994E+00 X_next( 5)= -2.429184E+01 X_next( 6)= 2.002600E+01 X_next( 7)= 7.777861E+00 residual X_resid( 1)= -3.893407E-02 residual X_resid( 2)= -2.879580E-02 residual X_resid( 3)= -1.070917E-01 residual X_resid( 4)= -6.623245E-02 itr 12, prev= 4.31379643061316E-01, residual= 2.41053997324212E-01 X_next( 1)= -3.478015E+00 X_next( 2)= 6.977826E+00 X_next( 3)= 1.287799E+00 X_next( 4)= 1.670232E+00 X_next( 5)= -2.426899E+01 X_next( 6)= 2.006132E+01 X_next( 7)= 7.782281E+00 simeq_newton4 did not converge < 1.00000000000000E-06 test case2 solution X( 1)= -3.47801515417990E+00 X( 2)= 6.97782615076833E+00 X( 3)= 1.28779948292732E+00 X( 4)= 1.67023169481777E+00 sum of residuals= 1.30794222345151E-01 end test case2 test_case_3 A( 1, 1)= 3.662632E-02 A( 1, 2)= 5.785426E-01 A( 1, 3)= 5.660214E-01 A( 1, 4)= 1.217688E-01 A( 1, 5)= 5.680125E-01 A( 1, 6)= 5.864237E-01 A( 1, 7)= 2.326706E-02 A( 1, 8)= 4.951762E-02 A( 1, 9)= 2.426066E-01 A( 1,10)= 4.884598E-01 A( 2, 1)= 5.430390E-01 A( 2, 2)= 8.572620E-01 A( 2, 3)= 2.925539E-03 A( 2, 4)= 1.695327E-01 A( 2, 5)= 3.362855E-01 A( 2, 6)= 9.501404E-01 A( 2, 7)= 1.036527E-02 A( 2, 8)= 2.091021E-01 A( 2, 9)= 3.791880E-01 A( 2,10)= 1.246447E-02 A( 3, 1)= 4.903669E-01 A( 3, 2)= 5.965996E-01 A( 3, 3)= 4.931436E-02 A( 3, 4)= 8.264277E-01 A( 3, 5)= 7.700864E-01 A( 3, 6)= 8.423222E-01 A( 3, 7)= 9.092780E-01 A( 3, 8)= 2.348641E-01 A( 3, 9)= 3.605372E-01 A( 3,10)= 5.494863E-01 A( 4, 1)= 2.166863E-01 A( 4, 2)= 8.464032E-01 A( 4, 3)= 4.987955E-01 A( 4, 4)= 2.563720E-01 A( 4, 5)= 8.437223E-01 A( 4, 6)= 4.406502E-01 A( 4, 7)= 8.642795E-03 A( 4, 8)= 2.594640E-01 A( 4, 9)= 8.108377E-01 A( 4,10)= 7.485342E-01 X_soln( 1)= 1.100000E+00 X_soln( 2)= 1.200000E+00 X_soln( 3)= 1.400000E+00 X_soln( 4)= 1.500000E+00 Y( 1)= 6.293103E+00 Y( 2)= 6.100440E+00 Y( 3)= 1.273644E+01 Y( 4)= 1.049535E+01 X_init( 1)= 7.000000E-01 X_init( 2)= 7.500000E-01 X_init( 3)= 8.000000E-01 X_init( 4)= 8.500000E-01 terms: x1 x2 x3 x4 power terms: x1 * x1 * x1 x2 * x2 * x2 x3 * x3 * x3 * x3 x4 * x4 * x4 * x4 x1 * x2 * x3 * x3 x2 * x3 * x4 * x4 simeq_newton4.adb running var1=( 1, 2, 3, 4, 1, 2, 3, 4, 3, 4) var2=( 0, 0, 0, 0, 1, 2, 3, 4, 3, 4) var3=( 0, 0, 0, 0, 1, 2, 3, 4, 2, 3) var4=( 0, 0, 0, 0, 0, 0, 3, 4, 1, 2) X_prev( 1)= 7.000000E-01 X_prev( 2)= 7.500000E-01 X_prev( 3)= 8.000000E-01 X_prev( 4)= 8.500000E-01 X_prev( 5)= 3.430000E-01 X_prev( 6)= 4.218750E-01 X_prev( 7)= 4.096000E-01 X_prev( 8)= 5.220062E-01 X_prev( 9)= 3.360000E-01 X_prev(10)= 4.335000E-01 residual X_resid( 1)= -4.506370E+00 residual X_resid( 2)= -4.168528E+00 residual X_resid( 3)= -9.729945E+00 residual X_resid( 4)= -7.880701E+00 itr 1, prev= 1.00000000000000E+12, residual= 2.62855439532732E+01 Ja( 1, 1)= 9.880559E-01 Ja( 1, 2)= 1.959150E+00 Ja( 1, 3)= 1.082146E+00 Ja( 1, 4)= 7.416378E-01 Ja( 2, 1)= 1.219389E+00 Ja( 2, 2)= 2.637705E+00 Ja( 2, 3)= 3.494257E-01 Ja( 2, 4)= 6.959058E-01 Ja( 3, 1)= 1.795452E+00 Ja( 3, 2)= 2.497142E+00 Ja( 3, 3)= 2.512120E+00 Ja( 3, 4)= 1.963847E+00 Ja( 4, 1)= 1.846160E+00 Ja( 4, 2)= 2.385909E+00 Ja( 4, 3)= 1.603212E+00 Ja( 4, 4)= 1.657250E+00 X_next( 1)= 5.381387E-01 X_next( 2)= 1.323471E+00 X_next( 3)= 1.484811E+00 X_next( 4)= 3.108682E+00 X_next( 5)= 1.558413E-01 X_next( 6)= 2.318160E+00 X_next( 7)= 4.860540E+00 X_next( 8)= 9.339107E+01 X_next( 9)= 1.570185E+00 X_next(10)= 1.899058E+01 residual X_resid( 1)= 1.155388E+01 residual X_resid( 2)= 1.852346E+01 residual X_resid( 3)= 3.038696E+01 residual X_resid( 4)= 3.319392E+01 itr 2, prev= 2.62855439532732E+01, residual= 9.36582169192537E+01 Ja( 1, 1)= 1.237983E+00 Ja( 1, 2)= 1.095682E+01 Ja( 1, 3)= 7.631144E+00 Ja( 1, 4)= 1.204009E+01 Ja( 2, 1)= 1.941595E+00 Ja( 2, 2)= 6.478718E+00 Ja( 2, 3)= 1.100050E+00 Ja( 2, 4)= 2.544921E+01 Ja( 3, 1)= 2.211381E+00 Ja( 3, 2)= 1.333514E+01 Ja( 3, 3)= 1.974584E+01 Ja( 3, 4)= 3.576308E+01 Ja( 4, 1)= 3.315564E+00 Ja( 4, 2)= 1.486466E+01 Ja( 4, 3)= 1.190056E+01 Ja( 4, 4)= 4.058107E+01 X_next( 1)= -4.378562E-01 X_next( 2)= 1.339628E+00 X_next( 3)= 1.363563E+00 X_next( 4)= 2.911324E+00 X_next( 5)= -8.394495E-02 X_next( 6)= 2.404102E+00 X_next( 7)= 3.457014E+00 X_next( 8)= 7.183933E+01 X_next( 9)= -1.090602E+00 X_next(10)= 1.548249E+01 residual X_resid( 1)= 7.890077E+00 residual X_resid( 2)= 1.240078E+01 residual X_resid( 3)= 2.041177E+01 residual X_resid( 4)= 2.233317E+01 itr 3, prev= 9.36582169192537E+01, residual= 6.30357982031984E+01 Ja( 1, 1)= 9.675999E-01 Ja( 1, 2)= 9.183505E+00 Ja( 1, 3)= 5.960076E+00 Ja( 1, 4)= 1.020460E+01 Ja( 2, 1)= 1.680928E+00 Ja( 2, 2)= 5.807994E+00 Ja( 2, 3)= -3.569947E-01 Ja( 2, 4)= 2.094117E+01 Ja( 3, 1)= 1.831303E+00 Ja( 3, 2)= 1.118857E+01 Ja( 3, 3)= 1.493279E+01 Ja( 3, 4)= 2.985265E+01 Ja( 4, 1)= 2.721572E+00 Ja( 4, 2)= 1.120971E+01 Ja( 4, 3)= 7.788579E+00 Ja( 4, 4)= 3.382777E+01 X_next( 1)= -2.268368E+00 X_next( 2)= 1.370272E+00 X_next( 3)= 1.169103E+00 X_next( 4)= 2.639325E+00 X_next( 5)= -1.167187E+01 X_next( 6)= 2.572883E+00 X_next( 7)= 1.868147E+00 X_next( 8)= 4.852567E+01 X_next( 9)= -4.248403E+00 X_next(10)= 1.115951E+01 residual X_resid( 1)= -2.854645E+00 residual X_resid( 2)= 1.507160E+00 residual X_resid( 3)= 8.233344E-02 residual X_resid( 4)= 2.339709E-01 itr 4, prev= 6.30357982031984E+01, residual= 4.67810890611528E+00 Ja( 1, 1)= 9.259115E+00 Ja( 1, 2)= 7.107672E+00 Ja( 1, 3)= 3.614048E+00 Ja( 1, 4)= 7.894001E+00 Ja( 2, 1)= 6.444280E+00 Ja( 2, 2)= 5.135212E+00 Ja( 2, 3)= -2.567707E+00 Ja( 2, 4)= 1.565284E+01 Ja( 3, 1)= 1.305304E+01 Ja( 3, 2)= 8.698556E+00 Ja( 3, 3)= 8.485906E+00 Ja( 3, 4)= 2.274557E+01 Ja( 4, 1)= 1.475940E+01 Ja( 4, 2)= 6.910702E+00 Ja( 4, 3)= 1.806064E+00 Ja( 4, 4)= 2.566787E+01 X_next( 1)= -1.797301E+00 X_next( 2)= 1.420711E+00 X_next( 3)= 1.167339E+00 X_next( 4)= 2.347308E+00 X_next( 5)= -5.805804E+00 X_next( 6)= 2.867588E+00 X_next( 7)= 1.856896E+00 X_next( 8)= 3.035850E+01 X_next( 9)= -3.479527E+00 X_next(10)= 9.137821E+00 residual X_resid( 1)= -1.040783E+00 residual X_resid( 2)= 4.768175E-01 residual X_resid( 3)= -2.431123E-01 residual X_resid( 4)= -2.215119E-01 itr 5, prev= 4.67810890611528E+00, residual= 1.98222497761508E+00 Ja( 1, 1)= 6.010842E+00 Ja( 1, 2)= 6.677017E+00 Ja( 1, 3)= 3.091393E+00 Ja( 1, 4)= 6.486523E+00 Ja( 2, 1)= 4.536036E+00 Ja( 2, 2)= 5.762088E+00 Ja( 2, 3)= -2.094069E+00 Ja( 2, 4)= 1.108412E+01 Ja( 3, 1)= 8.651166E+00 Ja( 3, 2)= 8.348290E+00 Ja( 3, 3)= 7.986898E+00 Ja( 3, 4)= 1.725491E+01 Ja( 4, 1)= 9.962851E+00 Ja( 4, 2)= 6.343266E+00 Ja( 4, 3)= 1.579462E+00 Ja( 4, 4)= 1.950723E+01 X_next( 1)= -1.440216E+00 X_next( 2)= 1.388949E+00 X_next( 3)= 1.195876E+00 X_next( 4)= 2.183421E+00 X_next( 5)= -2.987327E+00 X_next( 6)= 2.679529E+00 X_next( 7)= 2.045245E+00 X_next( 8)= 2.272741E+01 X_next( 9)= -2.860793E+00 X_next(10)= 7.918583E+00 residual X_resid( 1)= -3.781662E-01 residual X_resid( 2)= 1.062579E-02 residual X_resid( 3)= -2.768120E-01 residual X_resid( 4)= -2.929696E-01 itr 6, prev= 1.98222497761508E+00, residual= 9.58573534460529E-01 Ja( 1, 1)= 4.053083E+00 Ja( 1, 2)= 6.257573E+00 Ja( 1, 3)= 2.798829E+00 Ja( 1, 4)= 5.726481E+00 Ja( 2, 1)= 3.388836E+00 Ja( 2, 2)= 5.646287E+00 Ja( 2, 3)= -1.657829E+00 Ja( 2, 4)= 8.966188E+00 Ja( 3, 1)= 5.998515E+00 Ja( 3, 2)= 7.861672E+00 Ja( 3, 3)= 8.183175E+00 Ja( 3, 4)= 1.459094E+01 Ja( 4, 1)= 7.077506E+00 Ja( 4, 2)= 5.994108E+00 Ja( 4, 3)= 1.635000E+00 Ja( 4, 4)= 1.648890E+01 X_next( 1)= -1.254446E+00 X_next( 2)= 1.373161E+00 X_next( 3)= 1.211666E+00 X_next( 4)= 2.124931E+00 X_next( 5)= -1.974039E+00 X_next( 6)= 2.589192E+00 X_next( 7)= 2.155415E+00 X_next( 8)= 2.038821E+01 X_next( 9)= -2.528942E+00 X_next(10)= 7.512659E+00 residual X_resid( 1)= -8.713475E-02 residual X_resid( 2)= -2.419461E-02 residual X_resid( 3)= -9.109493E-02 residual X_resid( 4)= -9.883182E-02 itr 7, prev= 9.58573534460529E-01, residual= 3.01256116895347E-01 Ja( 1, 1)= 3.207249E+00 Ja( 1, 2)= 6.121363E+00 Ja( 1, 3)= 2.747445E+00 Ja( 1, 4)= 5.476091E+00 Ja( 2, 1)= 2.895047E+00 Ja( 2, 2)= 5.601779E+00 Ja( 2, 3)= -1.428890E+00 Ja( 2, 4)= 8.282809E+00 Ja( 3, 1)= 4.852708E+00 Ja( 3, 2)= 7.703652E+00 Ja( 3, 3)= 8.421291E+00 Ja( 3, 4)= 1.372569E+01 Ja( 4, 1)= 5.834453E+00 Ja( 4, 2)= 5.941003E+00 Ja( 4, 3)= 1.816713E+00 Ja( 4, 4)= 1.550721E+01 X_next( 1)= -1.205637E+00 X_next( 2)= 1.369929E+00 X_next( 3)= 1.215577E+00 X_next( 4)= 2.113595E+00 X_next( 5)= -1.752464E+00 X_next( 6)= 2.570951E+00 X_next( 7)= 2.183385E+00 X_next( 8)= 1.995664E+01 X_next( 9)= -2.440504E+00 X_next(10)= 7.439165E+00 residual X_resid( 1)= -6.384790E-03 residual X_resid( 2)= -2.522283E-03 residual X_resid( 3)= -7.422441E-03 residual X_resid( 4)= -8.074334E-03 itr 8, prev= 3.01256116895347E-01, residual= 2.44038488562168E-02 Ja( 1, 1)= 3.004641E+00 Ja( 1, 2)= 6.100473E+00 Ja( 1, 3)= 2.748336E+00 Ja( 1, 4)= 5.430394E+00 Ja( 2, 1)= 2.777041E+00 Ja( 2, 2)= 5.598830E+00 Ja( 2, 3)= -1.368907E+00 Ja( 2, 4)= 8.154671E+00 Ja( 3, 1)= 4.578282E+00 Ja( 3, 2)= 7.680571E+00 Ja( 3, 3)= 8.497276E+00 Ja( 3, 4)= 1.356483E+01 Ja( 4, 1)= 5.537222E+00 Ja( 4, 2)= 5.947609E+00 Ja( 4, 3)= 1.885994E+00 Ja( 4, 4)= 1.532503E+01 X_next( 1)= -1.202076E+00 X_next( 2)= 1.369727E+00 X_next( 3)= 1.215857E+00 X_next( 4)= 2.112874E+00 X_next( 5)= -1.736984E+00 X_next( 6)= 2.569817E+00 X_next( 7)= 2.185395E+00 X_next( 8)= 1.992942E+01 X_next( 9)= -2.434058E+00 X_next(10)= 7.434708E+00 residual X_resid( 1)= -8.691596E-05 residual X_resid( 2)= -3.608578E-05 residual X_resid( 3)= -1.027846E-04 residual X_resid( 4)= -1.118020E-04 itr 9, prev= 2.44038488562168E-02, residual= 3.37588375579401E-04 Ja( 1, 1)= 2.990187E+00 Ja( 1, 2)= 6.099381E+00 Ja( 1, 3)= 2.748769E+00 Ja( 1, 4)= 5.427593E+00 Ja( 2, 1)= 2.768633E+00 Ja( 2, 2)= 5.598912E+00 Ja( 2, 3)= -1.364549E+00 Ja( 2, 4)= 8.146569E+00 Ja( 3, 1)= 4.558705E+00 Ja( 3, 2)= 7.679429E+00 Ja( 3, 3)= 8.503153E+00 Ja( 3, 4)= 1.355476E+01 Ja( 4, 1)= 5.516037E+00 Ja( 4, 2)= 5.948645E+00 Ja( 4, 3)= 1.891585E+00 Ja( 4, 4)= 1.531365E+01 X_next( 1)= -1.202028E+00 X_next( 2)= 1.369725E+00 X_next( 3)= 1.215861E+00 X_next( 4)= 2.112865E+00 X_next( 5)= -1.736776E+00 X_next( 6)= 2.569803E+00 X_next( 7)= 2.185422E+00 X_next( 8)= 1.992906E+01 X_next( 9)= -2.433972E+00 X_next(10)= 7.434650E+00 residual X_resid( 1)= -4.288657E-07 residual X_resid( 2)= -1.783061E-07 residual X_resid( 3)= -5.074145E-07 residual X_resid( 4)= -5.519276E-07 itr 10, prev= 3.37588375579401E-04, residual= 1.66651397570661E-06 simeq_newton4 found solution test case3 solution X( 1)= -1.20202804077262E+00 X( 2)= 1.36972460795337E+00 X( 3)= 1.21586069402545E+00 X( 4)= 2.11286485539314E+00 sum of residuals= 1.66651397837114E-06 end test case3 test_case_big terms: x1 x2 x3 x4 power terms: x1 * x1 x1 * x2 x1 * x3 x1 * x4 x2 * x2 x2 * x3 x2 * x4 x3 * x3 x3 * x4 x4 * x4 x1 * x1 * x1 x1 * x1 * x2 x1 * x1 * x3 x1 * x1 * x4 x1 * x2 * x2 x1 * x2 * x3 x1 * x2 * x4 x1 * x3 * x3 x1 * x3 * x4 x1 * x4 * x4 x2 * x2 * x2 x2 * x2 * x3 x2 * x2 * x4 x2 * x3 * x3 x2 * x3 * x4 x2 * x4 * x4 x3 * x3 * x3 x3 * x3 * x4 x3 * x4 * x4 x4 * x4 * x4 x1 * x1 * x1 * x1 x1 * x1 * x1 * x2 x1 * x1 * x1 * x3 x1 * x1 * x1 * x4 x1 * x1 * x2 * x2 x1 * x1 * x2 * x3 x1 * x1 * x2 * x4 x1 * x1 * x3 * x3 x1 * x1 * x3 * x4 x1 * x1 * x4 * x4 x1 * x2 * x2 * x2 x1 * x2 * x2 * x3 x1 * x2 * x2 * x4 x1 * x2 * x3 * x3 x1 * x2 * x3 * x4 x1 * x2 * x4 * x4 x1 * x3 * x3 * x3 x1 * x3 * x3 * x4 x1 * x3 * x4 * x4 x1 * x4 * x4 * x4 x2 * x2 * x2 * x2 x2 * x2 * x2 * x3 x2 * x2 * x2 * x4 x2 * x2 * x3 * x3 x2 * x2 * x3 * x4 x2 * x2 * x4 * x4 x2 * x3 * x3 * x3 x2 * x3 * x3 * x4 x2 * x3 * x4 * x4 x2 * x4 * x4 * x4 x3 * x3 * x3 * x3 x3 * x3 * x3 * x4 x3 * x3 * x4 * x4 x3 * x4 * x4 * x4 x4 * x4 * x4 * x4 simeq_newton4.adb running itr 1, prev= 1.00000000000000E+12, residual= 1.42108547152020E-14 simeq_newton4 found solution test case_big solution init = solution X( 1)= 1.10000000000000E+00, err= 0.00000000000000E+00 X( 2)= 1.20000000000000E+00, err= 0.00000000000000E+00 X( 3)= 1.40000000000000E+00, err= 0.00000000000000E+00 X( 4)= 1.50000000000000E+00, err= 0.00000000000000E+00 sum of residuals= 3.16777670766878E+02 simeq_newton4.adb running itr 1, prev= 1.00000000000000E+12, residual= 1.95296594562390E+02 itr 2, prev= 1.95296594562390E+02, residual= 4.97926054345928E+01 itr 3, prev= 4.97926054345928E+01, residual= 1.73632740039519E+01 itr 4, prev= 1.73632740039519E+01, residual= 4.97593672407086E+00 itr 5, prev= 4.97593672407086E+00, residual= 1.09383918418345E-01 itr 6, prev= 1.09383918418345E-01, residual= 1.70645636382858E-02 itr 7, prev= 1.70645636382858E-02, residual= 5.36766188858451E-04 itr 8, prev= 5.36766188858451E-04, residual= 5.78791274392643E-05 itr 9, prev= 5.78791274392643E-05, residual= 4.64364819663388E-06 itr 10, prev= 4.64364819663388E-06, residual= 4.08742749868907E-07 simeq_newton4 found solution test case_big solution, init 1.0 X( 1)= 1.09999997678052E+00, err=-2.32194758975623E-08 X( 2)= 1.20000005056627E+00, err= 5.05662725025502E-08 X( 3)= 1.39999999611008E+00, err=-3.88992105193608E-09 X( 4)= 1.49999998414348E+00, err=-1.58565247598830E-08 sum of residuals= 4.08742749868907E-07 end test case_big end test_simeq_newton4.adb