test_simeq_newton3.adb running test_case_1 A( 1, 1)= 1.000000E+00 A( 1, 2)= 1.000000E+00 A( 1, 3)= 1.000000E+00 A( 1, 4)= 0.000000E+00 A( 2, 1)= 1.000000E+00 A( 2, 2)= 1.000000E+00 A( 2, 3)= 0.000000E+00 A( 2, 4)= 1.000000E+00 X_soln( 1)= 1.000000E+00 X_soln( 2)= 2.000000E+00 Y( 1)= 4.000000E+00 Y( 2)= 7.000000E+00 X( 1)= 5.000000E-01 X( 2)= 5.000000E-01 simeq_newton3.adb running var1=( 1, 2, 1, 2) var2=( 0, 0, 1, 2) var3=( 0, 0, 0, 0) X_prev( 1)= 5.000000E-01 X_prev( 2)= 5.000000E-01 X_prev( 3)= 2.500000E-01 X_prev( 4)= 2.500000E-01 residual X_resid( 1)= -2.750000E+00 residual X_resid( 2)= -5.750000E+00 itr 1, prev= 1.00000000000000E+12, residual= 8.50000000000000E+00 Ja( 1, 1)= 2.000000E+00 Ja( 1, 2)= 1.000000E+00 Ja( 2, 1)= 1.000000E+00 Ja( 2, 2)= 2.000000E+00 JaInv( 1, 1)= 6.666667E-01 JaInv( 1, 2)= -3.333333E-01 JaInv( 2, 1)= -3.333333E-01 JaInv( 2, 2)= 6.666667E-01 change X_tmp2( 1)= 8.333333E-02 change X_tmp2( 2)= -2.916667E+00 X_next( 1)= 4.166667E-01 X_next( 2)= 3.416667E+00 X_next( 3)= 1.736111E-01 X_next( 4)= 1.167361E+01 residual X_resid( 1)= 6.944444E-03 residual X_resid( 2)= 8.506944E+00 itr 2, prev= 8.50000000000000E+00, residual= 8.51388888888889E+00 Ja( 1, 1)= 1.833333E+00 Ja( 1, 2)= 1.000000E+00 Ja( 2, 1)= 1.000000E+00 Ja( 2, 2)= 7.833333E+00 JaInv( 1, 1)= 5.862786E-01 JaInv( 1, 2)= -7.484407E-02 JaInv( 2, 1)= -7.484407E-02 JaInv( 2, 2)= 1.372141E-01 change X_tmp2( 1)= -6.326230E-01 change X_tmp2( 2)= 1.166753E+00 X_next( 1)= 1.049290E+00 X_next( 2)= 2.249913E+00 X_next( 3)= 1.101009E+00 X_next( 4)= 5.062110E+00 residual X_resid( 1)= 4.002119E-01 residual X_resid( 2)= 1.361313E+00 itr 3, prev= 8.51388888888889E+00, residual= 1.76152511359112E+00 Ja( 1, 1)= 3.098579E+00 Ja( 1, 2)= 1.000000E+00 Ja( 2, 1)= 1.000000E+00 Ja( 2, 2)= 5.499827E+00 JaInv( 1, 1)= 3.428467E-01 JaInv( 1, 2)= -6.233773E-02 JaInv( 2, 1)= -6.233773E-02 JaInv( 2, 2)= 1.931584E-01 change X_tmp2( 1)= 5.235015E-02 change X_tmp2( 2)= 2.380008E-01 X_next( 1)= 9.969395E-01 X_next( 2)= 2.011913E+00 X_next( 3)= 9.938884E-01 X_next( 4)= 4.047792E+00 residual X_resid( 1)= 2.740538E-03 residual X_resid( 2)= 5.664437E-02 itr 4, prev= 1.76152511359112E+00, residual= 5.93849126638721E-02 Ja( 1, 1)= 2.993879E+00 Ja( 1, 2)= 1.000000E+00 Ja( 2, 1)= 1.000000E+00 Ja( 2, 2)= 5.023825E+00 JaInv( 1, 1)= 3.578038E-01 JaInv( 1, 2)= -7.122139E-02 JaInv( 2, 1)= -7.122139E-02 JaInv( 2, 2)= 2.132282E-01 change X_tmp2( 1)= -3.053716E-03 change X_tmp2( 2)= 1.188300E-02 X_next( 1)= 9.999932E-01 X_next( 2)= 2.000030E+00 X_next( 3)= 9.999865E-01 X_next( 4)= 4.000118E+00 residual X_resid( 1)= 9.325184E-06 residual X_resid( 2)= 1.412056E-04 itr 5, prev= 5.93849126638721E-02, residual= 1.50530762461187E-04 Ja( 1, 1)= 2.999986E+00 Ja( 1, 2)= 1.000000E+00 Ja( 2, 1)= 1.000000E+00 Ja( 2, 2)= 5.000059E+00 JaInv( 1, 1)= 3.571443E-01 JaInv( 1, 2)= -7.142801E-02 JaInv( 2, 1)= -7.142801E-02 JaInv( 2, 2)= 2.142831E-01 change X_tmp2( 1)= -6.755597E-06 change X_tmp2( 2)= 2.959188E-05 X_next( 1)= 1.000000E+00 X_next( 2)= 2.000000E+00 X_next( 3)= 1.000000E+00 X_next( 4)= 4.000000E+00 residual X_resid( 1)= 4.563816E-11 residual X_resid( 2)= 8.756800E-10 itr 6, prev= 1.50530762461187E-04, residual= 9.21318132895976E-10 simeq_newton3 found solution test case1 solution X( 1)= 9.99999999953751E-01 X( 2)= 2.00000000018439E+00 sum of errors= 9.21318132895976E-10 test_case_2 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( 1, 5)= 1.223651E-01 A( 1, 6)= 5.898427E-01 A( 1, 7)= 4.864346E-01 A( 2, 1)= 5.060496E-01 A( 2, 2)= 1.750984E-01 A( 2, 3)= 8.782248E-01 A( 2, 4)= 3.246958E-01 A( 2, 5)= 1.620386E-01 A( 2, 6)= 3.833710E-01 A( 2, 7)= 3.158238E-01 A( 3, 1)= 5.133576E-02 A( 3, 2)= 8.000549E-01 A( 3, 3)= 5.223173E-01 A( 3, 4)= 5.863344E-01 A( 3, 5)= 5.226468E-01 A( 3, 6)= 1.255568E-01 A( 3, 7)= 2.337342E-01 A( 4, 1)= 3.708149E-01 A( 4, 2)= 2.857049E-01 A( 4, 3)= 8.418347E-01 A( 4, 4)= 7.152572E-01 A( 4, 5)= 3.275851E-01 A( 4, 6)= 7.227268E-01 A( 4, 7)= 8.689656E-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.814518E+00 Y( 2)= 4.412558E+00 Y( 3)= 4.318721E+00 Y( 4)= 7.591667E+00 X( 1)= 1.000000E+00 X( 2)= 1.000000E+00 X( 3)= 1.000000E+00 X( 4)= 1.000000E+00 simeq_newton3.adb running var1=( 1, 2, 3, 4, 2, 3, 4) var2=( 0, 0, 0, 0, 1, 1, 4) var3=( 0, 0, 0, 0, 0, 1, 4) X_prev( 1)= 1.000000E+00 X_prev( 2)= 1.000000E+00 X_prev( 3)= 1.000000E+00 X_prev( 4)= 1.000000E+00 X_prev( 5)= 1.000000E+00 X_prev( 6)= 1.000000E+00 X_prev( 7)= 1.000000E+00 residual X_resid( 1)= -2.556949E+00 residual X_resid( 2)= -1.667256E+00 residual X_resid( 3)= -1.476741E+00 residual X_resid( 4)= -3.458778E+00 itr 1, prev= 1.00000000000000E+12, residual= 9.15972339520032E+00 Ja( 1, 1)= 2.093128E+00 Ja( 1, 2)= 7.666601E-01 Ja( 1, 3)= 1.255694E+00 Ja( 1, 4)= 2.417007E+00 Ja( 2, 1)= 1.434830E+00 Ja( 2, 2)= 3.371370E-01 Ja( 2, 3)= 1.261596E+00 Ja( 2, 4)= 1.272167E+00 Ja( 3, 1)= 8.250963E-01 Ja( 3, 2)= 1.322702E+00 Ja( 3, 3)= 6.478741E-01 Ja( 3, 4)= 1.287537E+00 Ja( 4, 1)= 2.143854E+00 Ja( 4, 2)= 6.132900E-01 Ja( 4, 3)= 1.564561E+00 Ja( 4, 4)= 3.322154E+00 JaInv( 1, 1)= 2.288384E+00 JaInv( 1, 2)= -3.653749E-01 JaInv( 1, 3)= -6.414436E-01 JaInv( 1, 4)= -1.276383E+00 JaInv( 2, 1)= 9.386802E-02 JaInv( 2, 2)= -7.308212E-02 JaInv( 2, 3)= 9.008118E-01 JaInv( 2, 4)= -3.894267E-01 JaInv( 3, 1)= -2.135010E+00 JaInv( 3, 2)= 1.859356E+00 JaInv( 3, 3)= 4.553021E-01 JaInv( 3, 4)= 6.648410E-01 JaInv( 4, 1)= -4.885911E-01 JaInv( 4, 2)= -6.263844E-01 JaInv( 4, 3)= 3.321766E-02 JaInv( 4, 4)= 8.834704E-01 change X_tmp2( 1)= 1.198634E-01 change X_tmp2( 2)= -1.014942E-01 change X_tmp2( 3)= -6.128118E-01 change X_tmp2( 4)= -8.111365E-01 X_next( 1)= 8.801366E-01 X_next( 2)= 1.101494E+00 X_next( 3)= 1.612812E+00 X_next( 4)= 1.811137E+00 X_next( 5)= 9.694653E-01 X_next( 6)= 1.249349E+00 X_next( 7)= 5.940918E+00 residual X_resid( 1)= 1.145266E+00 residual X_resid( 2)= 7.425230E-01 residual X_resid( 3)= 5.641964E-01 residual X_resid( 4)= 2.085526E+00 itr 2, prev= 9.15972339520032E+00, residual= 4.53751158658076E+00 Ja( 1, 1)= 2.600419E+00 Ja( 1, 2)= 7.519930E-01 Ja( 1, 3)= 1.122767E+00 Ja( 1, 4)= 5.744534E+00 Ja( 2, 1)= 1.772920E+00 Ja( 2, 2)= 3.177145E-01 Ja( 2, 3)= 1.175199E+00 Ja( 2, 4)= 3.432607E+00 Ja( 3, 1)= 9.834827E-01 Ja( 3, 2)= 1.260055E+00 Ja( 3, 3)= 6.195787E-01 Ja( 3, 4)= 2.886430E+00 Ja( 4, 1)= 2.783461E+00 Ja( 4, 2)= 5.740245E-01 Ja( 4, 3)= 1.401688E+00 Ja( 4, 4)= 9.266441E+00 JaInv( 1, 1)= 1.617090E+00 JaInv( 1, 2)= -4.732991E-01 JaInv( 1, 3)= -5.464590E-01 JaInv( 1, 4)= -6.569366E-01 JaInv( 2, 1)= 1.390557E-01 JaInv( 2, 2)= -3.058482E-01 JaInv( 2, 3)= 9.036311E-01 JaInv( 2, 4)= -2.543827E-01 JaInv( 3, 1)= -1.851038E+00 JaInv( 3, 2)= 2.108714E+00 JaInv( 3, 3)= 4.732409E-01 JaInv( 3, 4)= 2.189613E-01 JaInv( 4, 1)= -2.143596E-01 JaInv( 4, 2)= -1.578583E-01 JaInv( 4, 3)= 3.658418E-02 JaInv( 4, 4)= 2.878844E-01 change X_tmp2( 1)= -1.778066E-01 change X_tmp2( 2)= -8.853991E-02 change X_tmp2( 3)= 1.694883E-01 change X_tmp2( 4)= 2.583191E-01 X_next( 1)= 1.057943E+00 X_next( 2)= 1.190034E+00 X_next( 3)= 1.443323E+00 X_next( 4)= 1.552817E+00 X_next( 5)= 1.258988E+00 X_next( 6)= 1.615431E+00 X_next( 7)= 3.744219E+00 residual X_resid( 1)= 1.655308E-01 residual X_resid( 2)= 1.087702E-01 residual X_resid( 3)= 8.801160E-02 residual X_resid( 4)= 2.998744E-01 itr 3, prev= 4.53751158658076E+00, residual= 6.62186950662226E-01 Ja( 1, 1)= 2.738022E+00 Ja( 1, 2)= 7.737503E-01 Ja( 1, 3)= 1.326029E+00 Ja( 1, 4)= 4.476438E+00 Ja( 2, 1)= 1.869661E+00 Ja( 2, 2)= 3.465260E-01 Ja( 2, 3)= 1.307310E+00 Ja( 2, 4)= 2.609279E+00 Ja( 3, 1)= 1.056742E+00 Ja( 3, 2)= 1.352986E+00 Ja( 3, 3)= 6.628460E-01 Ja( 3, 4)= 2.277104E+00 Ja( 4, 1)= 2.967794E+00 Ja( 4, 2)= 6.322713E-01 Ja( 4, 3)= 1.650742E+00 Ja( 4, 4)= 7.001117E+00 JaInv( 1, 1)= 1.799615E+00 JaInv( 1, 2)= -6.087784E-01 JaInv( 1, 3)= -5.207037E-01 JaInv( 1, 4)= -7.544074E-01 JaInv( 2, 1)= 7.581056E-02 JaInv( 2, 2)= -1.870929E-01 JaInv( 2, 3)= 8.652486E-01 JaInv( 2, 4)= -2.601649E-01 JaInv( 3, 1)= -1.997660E+00 JaInv( 3, 2)= 2.146544E+00 JaInv( 3, 3)= 4.358675E-01 JaInv( 3, 4)= 3.355114E-01 JaInv( 4, 1)= -2.986947E-01 JaInv( 4, 2)= -2.311585E-01 JaInv( 4, 3)= 3.981714E-02 JaInv( 4, 4)= 4.070176E-01 change X_tmp2( 1)= -4.038077E-02 change X_tmp2( 2)= -9.666014E-03 change X_tmp2( 3)= 4.177840E-02 change X_tmp2( 4)= 5.097221E-02 X_next( 1)= 1.098324E+00 X_next( 2)= 1.199700E+00 X_next( 3)= 1.401545E+00 X_next( 4)= 1.501845E+00 X_next( 5)= 1.317659E+00 X_next( 6)= 1.690706E+00 X_next( 7)= 3.387471E+00 residual X_resid( 1)= 5.113380E-03 residual X_resid( 2)= 3.351640E-03 residual X_resid( 3)= 2.840789E-03 residual X_resid( 4)= 9.602108E-03 itr 4, prev= 6.62186950662226E-01, residual= 2.09079174881861E-02 Ja( 1, 1)= 2.753829E+00 Ja( 1, 2)= 7.786915E-01 Ja( 1, 3)= 1.377387E+00 Ja( 1, 4)= 4.249220E+00 Ja( 2, 1)= 1.880732E+00 Ja( 2, 2)= 3.530693E-01 Ja( 2, 3)= 1.340691E+00 Ja( 2, 4)= 2.461755E+00 Ja( 3, 1)= 1.064907E+00 Ja( 3, 2)= 1.374090E+00 Ja( 3, 3)= 6.737784E-01 Ja( 3, 4)= 2.167924E+00 Ja( 4, 1)= 2.988878E+00 Ja( 4, 2)= 6.454994E-01 Ja( 4, 3)= 1.713671E+00 Ja( 4, 4)= 6.595215E+00 JaInv( 1, 1)= 1.850116E+00 JaInv( 1, 2)= -6.399326E-01 JaInv( 1, 3)= -5.159405E-01 JaInv( 1, 4)= -7.835491E-01 JaInv( 2, 1)= 6.166943E-02 JaInv( 2, 2)= -1.590529E-01 JaInv( 2, 3)= 8.567468E-01 JaInv( 2, 4)= -2.619869E-01 JaInv( 3, 1)= -2.029020E+00 JaInv( 3, 2)= 2.150294E+00 JaInv( 3, 3)= 4.260482E-01 JaInv( 3, 4)= 3.646002E-01 JaInv( 4, 1)= -3.172766E-01 JaInv( 4, 2)= -2.531452E-01 JaInv( 4, 3)= 3.926287E-02 JaInv( 4, 4)= 4.376264E-01 change X_tmp2( 1)= -1.673877E-03 change X_tmp2( 2)= -2.995384E-04 change X_tmp2( 3)= 1.543104E-03 change X_tmp2( 4)= 1.842867E-03 X_next( 1)= 1.099998E+00 X_next( 2)= 1.200000E+00 X_next( 3)= 1.400002E+00 X_next( 4)= 1.500002E+00 X_next( 5)= 1.319997E+00 X_next( 6)= 1.693996E+00 X_next( 7)= 3.375016E+00 residual X_resid( 1)= 6.468528E-06 residual X_resid( 2)= 4.240470E-06 residual X_resid( 3)= 3.617191E-06 residual X_resid( 4)= 1.218962E-05 itr 5, prev= 2.09079174881861E-02, residual= 2.65158073409211E-05 Ja( 1, 1)= 2.754630E+00 Ja( 1, 2)= 7.788963E-01 Ja( 1, 3)= 1.379558E+00 Ja( 1, 4)= 4.241147E+00 Ja( 2, 1)= 1.881278E+00 Ja( 2, 2)= 3.533405E-01 Ja( 2, 3)= 1.342102E+00 Ja( 2, 4)= 2.456513E+00 Ja( 3, 1)= 1.065227E+00 Ja( 3, 2)= 1.374965E+00 Ja( 3, 3)= 6.742404E-01 Ja( 3, 4)= 2.164045E+00 Ja( 4, 1)= 2.989914E+00 Ja( 4, 2)= 6.460478E-01 Ja( 4, 3)= 1.716331E+00 Ja( 4, 4)= 6.580794E+00 JaInv( 1, 1)= 1.852017E+00 JaInv( 1, 2)= -6.412016E-01 JaInv( 1, 3)= -5.156844E-01 JaInv( 1, 4)= -7.846460E-01 JaInv( 2, 1)= 6.115585E-02 JaInv( 2, 2)= -1.579419E-01 JaInv( 2, 3)= 8.563716E-01 JaInv( 2, 4)= -2.620674E-01 JaInv( 3, 1)= -2.030165E+00 JaInv( 3, 2)= 2.150431E+00 JaInv( 3, 3)= 4.256030E-01 JaInv( 3, 4)= 3.657068E-01 JaInv( 4, 1)= -3.179630E-01 JaInv( 4, 2)= -2.540232E-01 JaInv( 4, 3)= 3.922319E-02 JaInv( 4, 4)= 4.388010E-01 change X_tmp2( 1)= -2.169034E-06 change X_tmp2( 2)= -3.710010E-07 change X_tmp2( 3)= 1.983974E-06 change X_tmp2( 4)= 2.356763E-06 X_next( 1)= 1.100000E+00 X_next( 2)= 1.200000E+00 X_next( 3)= 1.400000E+00 X_next( 4)= 1.500000E+00 X_next( 5)= 1.320000E+00 X_next( 6)= 1.694000E+00 X_next( 7)= 3.375000E+00 residual X_resid( 1)= 1.055867E-11 residual X_resid( 2)= 6.918910E-12 residual X_resid( 3)= 5.901946E-12 residual X_resid( 4)= 1.990141E-11 itr 6, prev= 2.65158073409211E-05, residual= 4.32809343919871E-11 simeq_newton3 found solution test case2 solution X( 1)= 1.09999999999646E+00 X( 2)= 1.19999999999939E+00 X( 3)= 1.40000000000323E+00 X( 4)= 1.50000000000385E+00 sum of errors= 4.32809343919871E-11 test_case_3 A( 1, 1)= 7.055294E-01 A( 1, 2)= 8.333382E-01 A( 1, 3)= 9.146897E-01 A( 1, 4)= 1.905640E-01 A( 1, 5)= 8.099251E-01 A( 1, 6)= 4.112297E-01 A( 1, 7)= 5.375608E-01 A( 1, 8)= 7.839017E-01 A( 1, 9)= 3.662632E-02 A( 1,10)= 5.785426E-01 A( 2, 1)= 5.660214E-01 A( 2, 2)= 1.217688E-01 A( 2, 3)= 5.680125E-01 A( 2, 4)= 5.864237E-01 A( 2, 5)= 2.326706E-02 A( 2, 6)= 4.951762E-02 A( 2, 7)= 2.426066E-01 A( 2, 8)= 4.884598E-01 A( 2, 9)= 5.430390E-01 A( 2,10)= 8.572620E-01 A( 3, 1)= 2.925539E-03 A( 3, 2)= 1.695327E-01 A( 3, 3)= 3.362855E-01 A( 3, 4)= 9.501404E-01 A( 3, 5)= 1.036527E-02 A( 3, 6)= 2.091021E-01 A( 3, 7)= 3.791880E-01 A( 3, 8)= 1.246447E-02 A( 3, 9)= 4.903669E-01 A( 3,10)= 5.965996E-01 A( 4, 1)= 4.931436E-02 A( 4, 2)= 8.264277E-01 A( 4, 3)= 7.700864E-01 A( 4, 4)= 8.423222E-01 A( 4, 5)= 9.092780E-01 A( 4, 6)= 2.348641E-01 A( 4, 7)= 3.605372E-01 A( 4, 8)= 5.494863E-01 A( 4, 9)= 2.166863E-01 A( 4,10)= 8.464032E-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)= 1.077746E+01 Y( 2)= 8.038235E+00 Y( 3)= 5.969981E+00 Y( 4)= 1.038086E+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 simeq_newton3.adb running var1=( 1, 2, 3, 4, 1, 2, 3, 4, 3, 4) var2=( 0, 0, 0, 0, 1, 2, 3, 4, 2, 3) var3=( 0, 0, 0, 0, 1, 2, 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)= 5.120000E-01 X_prev( 8)= 6.141250E-01 X_prev( 9)= 4.200000E-01 X_prev(10)= 5.100000E-01 residual X_resid( 1)= -7.246481E+00 residual X_resid( 2)= -5.479482E+00 residual X_resid( 3)= -3.960346E+00 residual X_resid( 4)= -6.938789E+00 itr 1, prev= 1.00000000000000E+12, residual= 2.36250984366996E+01 Ja( 1, 1)= 1.918095E+00 Ja( 1, 2)= 1.941208E+00 Ja( 1, 3)= 2.334856E+00 Ja( 1, 4)= 2.236797E+00 Ja( 2, 1)= 9.260474E-01 Ja( 2, 2)= 1.092370E+00 Ja( 2, 3)= 1.865417E+00 Ja( 2, 4)= 2.159517E+00 Ja( 3, 1)= 3.123826E-01 Ja( 3, 2)= 1.202686E+00 Ja( 3, 3)= 1.702101E+00 Ja( 3, 4)= 1.335117E+00 Ja( 4, 1)= 1.515965E+00 Ja( 4, 2)= 1.919659E+00 Ja( 4, 3)= 2.115660E+00 Ja( 4, 4)= 2.541176E+00 JaInv( 1, 1)= 1.088675E+00 JaInv( 1, 2)= 1.865141E-01 JaInv( 1, 3)= -8.925732E-01 JaInv( 1, 4)= -6.478244E-01 JaInv( 2, 1)= -5.953872E-01 JaInv( 2, 2)= -1.884008E+00 JaInv( 2, 3)= 6.918468E-01 JaInv( 2, 4)= 1.761630E+00 JaInv( 3, 1)= 1.088131E+00 JaInv( 3, 2)= 7.721551E-01 JaInv( 3, 3)= 7.342605E-01 JaInv( 3, 4)= -1.999757E+00 JaInv( 4, 1)= -1.105617E+00 JaInv( 4, 2)= 6.690943E-01 JaInv( 4, 3)= -6.014721E-01 JaInv( 4, 4)= 1.114112E+00 change X_tmp2( 1)= -8.810486E-01 change X_tmp2( 2)= -3.256820E-01 change X_tmp2( 3)= -1.148166E+00 change X_tmp2( 4)= -1.003009E+00 X_next( 1)= 1.581049E+00 X_next( 2)= 1.075682E+00 X_next( 3)= 1.948166E+00 X_next( 4)= 1.853009E+00 X_next( 5)= 3.952171E+00 X_next( 6)= 1.244663E+00 X_next( 7)= 7.393971E+00 X_next( 8)= 6.362573E+00 X_next( 9)= 3.313256E+00 X_next(10)= 3.883179E+00 residual X_resid( 1)= 8.412586E+00 residual X_resid( 2)= 5.364292E+00 residual X_resid( 3)= 3.758419E+00 residual X_resid( 4)= 7.699737E+00 itr 2, prev= 2.36250984366996E+01, residual= 2.52350348986442E+01 Ja( 1, 1)= 6.856029E+00 Ja( 1, 2)= 4.462165E+00 Ja( 1, 3)= 8.250853E+00 Ja( 1, 4)= 9.477879E+00 Ja( 2, 1)= 1.878501E+00 Ja( 2, 2)= 5.060987E+00 Ja( 2, 3)= 5.962629E+00 Ja( 2, 4)= 7.414497E+00 Ja( 3, 1)= 1.108272E+00 Ja( 3, 2)= 4.559491E+00 Ja( 3, 3)= 6.676880E+00 Ja( 3, 4)= 2.328774E+00 Ja( 4, 1)= 7.322210E+00 Ja( 4, 2)= 5.364620E+00 Ja( 4, 3)= 6.930793E+00 Ja( 4, 4)= 8.276271E+00 JaInv( 1, 1)= -4.349274E-02 JaInv( 1, 2)= -1.978971E-01 JaInv( 1, 3)= -7.432618E-03 JaInv( 1, 4)= 2.291896E-01 JaInv( 2, 1)= -4.969901E-01 JaInv( 2, 2)= 1.540123E-01 JaInv( 2, 3)= 4.102593E-02 JaInv( 2, 4)= 4.196272E-01 JaInv( 3, 1)= 3.119326E-01 JaInv( 3, 2)= -1.392399E-01 JaInv( 3, 3)= 1.835944E-01 JaInv( 3, 4)= -2.841395E-01 JaInv( 4, 1)= 9.940294E-02 JaInv( 4, 2)= 1.918581E-01 JaInv( 4, 3)= -1.737643E-01 JaInv( 4, 4)= -1.159946E-01 change X_tmp2( 1)= 3.093007E-01 change X_tmp2( 2)= 3.040658E-02 change X_tmp2( 3)= 3.794615E-01 change X_tmp2( 4)= 3.192119E-01 X_next( 1)= 1.271748E+00 X_next( 2)= 1.045275E+00 X_next( 3)= 1.568704E+00 X_next( 4)= 1.533797E+00 X_next( 5)= 2.056852E+00 X_next( 6)= 1.142069E+00 X_next( 7)= 3.860320E+00 X_next( 8)= 3.608311E+00 X_next( 9)= 2.085321E+00 X_next(10)= 2.515011E+00 residual X_resid( 1)= 1.288710E+00 residual X_resid( 2)= 6.912794E-01 residual X_resid( 3)= 4.877215E-01 residual X_resid( 4)= 1.139251E+00 itr 3, prev= 2.52350348986442E+01, residual= 3.60696126115965E+00 Ja( 1, 1)= 4.695366E+00 Ja( 1, 2)= 3.646354E+00 Ja( 1, 3)= 5.859463E+00 Ja( 1, 4)= 6.671684E+00 Ja( 2, 1)= 1.569350E+00 Ja( 2, 2)= 3.430075E+00 Ja( 2, 3)= 4.455329E+00 Ja( 2, 4)= 5.439455E+00 Ja( 3, 1)= 8.572865E-01 Ja( 3, 2)= 3.268671E+00 Ja( 3, 3)= 4.743991E+00 Ja( 3, 4)= 2.016371E+00 Ja( 4, 1)= 4.816463E+00 Ja( 4, 2)= 4.065063E+00 Ja( 4, 3)= 5.076787E+00 Ja( 4, 4)= 6.108250E+00 JaInv( 1, 1)= 5.636691E-02 JaInv( 1, 2)= -3.756776E-01 JaInv( 1, 3)= -1.380819E-02 JaInv( 1, 4)= 2.775364E-01 JaInv( 2, 1)= -1.036854E+00 JaInv( 2, 2)= 1.858471E-01 JaInv( 2, 3)= 1.102186E-01 JaInv( 2, 4)= 9.306124E-01 JaInv( 3, 1)= 6.646026E-01 JaInv( 3, 2)= -2.064225E-01 JaInv( 3, 3)= 2.534189E-01 JaInv( 3, 4)= -6.257404E-01 JaInv( 4, 1)= 9.320858E-02 JaInv( 4, 2)= 3.441117E-01 JaInv( 4, 3)= -2.730885E-01 JaInv( 4, 4)= -1.543799E-01 change X_tmp2( 1)= 1.223914E-01 change X_tmp2( 2)= -9.377484E-02 change X_tmp2( 3)= 1.245069E-01 change X_tmp2( 4)= 4.892766E-02 X_next( 1)= 1.149357E+00 X_next( 2)= 1.139050E+00 X_next( 3)= 1.444197E+00 X_next( 4)= 1.484870E+00 X_next( 5)= 1.518323E+00 X_next( 6)= 1.477844E+00 X_next( 7)= 3.012171E+00 X_next( 8)= 3.273897E+00 X_next( 9)= 1.890707E+00 X_next(10)= 2.442631E+00 residual X_resid( 1)= 9.212726E-02 residual X_resid( 2)= 1.259478E-03 residual X_resid( 3)= 1.514260E-02 residual X_resid( 4)= 6.983057E-02 itr 4, prev= 3.60696126115965E+00, residual= 1.78359913564033E-01 Ja( 1, 1)= 3.975563E+00 Ja( 1, 2)= 3.735419E+00 Ja( 1, 3)= 5.304734E+00 Ja( 1, 4)= 6.327404E+00 Ja( 2, 1)= 1.551537E+00 Ja( 2, 2)= 3.054247E+00 Ja( 2, 3)= 4.246886E+00 Ja( 2, 4)= 5.227555E+00 Ja( 3, 1)= 8.506639E-01 Ja( 3, 2)= 3.076756E+00 Ja( 3, 3)= 4.359939E+00 Ja( 3, 4)= 2.014001E+00 Ja( 4, 1)= 4.009291E+00 Ja( 4, 2)= 3.915333E+00 Ja( 4, 3)= 4.741247E+00 Ja( 4, 4)= 5.869252E+00 JaInv( 1, 1)= 4.125068E-01 JaInv( 1, 2)= -5.597016E-01 JaInv( 1, 3)= -2.427209E-02 JaInv( 1, 4)= 6.213024E-02 JaInv( 2, 1)= -1.854177E+00 JaInv( 2, 2)= 2.848590E-01 JaInv( 2, 3)= 1.286987E-01 JaInv( 2, 4)= 1.701037E+00 JaInv( 3, 1)= 1.255148E+00 JaInv( 3, 2)= -2.881907E-01 JaInv( 3, 3)= 2.796152E-01 JaInv( 3, 4)= -1.192390E+00 JaInv( 4, 1)= -5.879860E-02 JaInv( 4, 2)= 4.251090E-01 JaInv( 4, 3)= -2.951499E-01 JaInv( 4, 4)= -4.358411E-02 change X_tmp2( 1)= 4.126924E-02 change X_tmp2( 2)= -4.972831E-02 change X_tmp2( 3)= 3.623916E-02 change X_tmp2( 4)= -1.239438E-02 X_next( 1)= 1.108087E+00 X_next( 2)= 1.188779E+00 X_next( 3)= 1.407958E+00 X_next( 4)= 1.497264E+00 X_next( 5)= 1.360573E+00 X_next( 6)= 1.679975E+00 X_next( 7)= 2.791061E+00 X_next( 8)= 3.356566E+00 X_next( 9)= 1.854662E+00 X_next(10)= 2.506047E+00 residual X_resid( 1)= 1.033465E-02 residual X_resid( 2)= -1.493832E-03 residual X_resid( 3)= 1.019692E-03 residual X_resid( 4)= 7.031560E-03 itr 5, prev= 1.78359913564033E-01, residual= 1.98797384507108E-02 Ja( 1, 1)= 3.750250E+00 Ja( 1, 2)= 3.853540E+00 Ja( 1, 3)= 5.189588E+00 Ja( 1, 4)= 6.430952E+00 Ja( 2, 1)= 1.560639E+00 Ja( 2, 2)= 2.986102E+00 Ja( 2, 3)= 4.251986E+00 Ja( 2, 4)= 5.306353E+00 Ja( 3, 1)= 8.618586E-01 Ja( 3, 2)= 3.078763E+00 Ja( 3, 3)= 4.299174E+00 Ja( 3, 4)= 2.032528E+00 Ja( 4, 1)= 3.761384E+00 Ja( 4, 2)= 3.944505E+00 Ja( 4, 3)= 4.706176E+00 Ja( 4, 4)= 5.954505E+00 JaInv( 1, 1)= 8.331515E-01 JaInv( 1, 2)= -6.991719E-01 JaInv( 1, 3)= -1.797918E-02 JaInv( 1, 4)= -2.706121E-01 JaInv( 2, 1)= -2.621533E+00 JaInv( 2, 2)= 4.344458E-01 JaInv( 2, 3)= 9.465276E-02 JaInv( 2, 4)= 2.411829E+00 JaInv( 3, 1)= 1.817104E+00 JaInv( 3, 2)= -3.890812E-01 JaInv( 3, 3)= 3.076558E-01 JaInv( 3, 4)= -1.720785E+00 JaInv( 4, 1)= -2.258397E-01 JaInv( 4, 2)= 4.613759E-01 JaInv( 4, 3)= -2.945021E-01 JaInv( 4, 4)= 1.012215E-01 change X_tmp2( 1)= 7.733620E-03 change X_tmp2( 2)= -1.068619E-02 change X_tmp2( 3)= 7.574271E-03 change X_tmp2( 4)= -2.611749E-03 X_next( 1)= 1.100354E+00 X_next( 2)= 1.199465E+00 X_next( 3)= 1.400384E+00 X_next( 4)= 1.499876E+00 X_next( 5)= 1.332284E+00 X_next( 6)= 1.725689E+00 X_next( 7)= 2.746258E+00 X_next( 8)= 3.374162E+00 X_next( 9)= 1.848276E+00 X_next(10)= 2.519358E+00 residual X_resid( 1)= 4.166147E-04 residual X_resid( 2)= -6.578165E-05 residual X_resid( 3)= 4.997389E-05 residual X_resid( 4)= 2.615332E-04 itr 6, prev= 1.98797384507108E-02, residual= 7.93903515246441E-04 Ja( 1, 1)= 3.708970E+00 Ja( 1, 2)= 3.879876E+00 Ja( 1, 3)= 5.166448E+00 Ja( 1, 4)= 6.452810E+00 Ja( 2, 1)= 1.562684E+00 Ja( 2, 2)= 2.972867E+00 Ja( 2, 3)= 4.254300E+00 Ja( 2, 4)= 5.322934E+00 Ja( 3, 1)= 8.642504E-01 Ja( 3, 2)= 3.080762E+00 Ja( 3, 3)= 4.287649E+00 Ja( 3, 4)= 2.036377E+00 Ja( 4, 1)= 3.716087E+00 Ja( 4, 2)= 3.951818E+00 Ja( 4, 3)= 4.699919E+00 Ja( 4, 4)= 5.972454E+00 JaInv( 1, 1)= 9.788219E-01 JaInv( 1, 2)= -7.396154E-01 JaInv( 1, 3)= -1.422457E-02 JaInv( 1, 4)= -3.935168E-01 JaInv( 2, 1)= -2.878605E+00 JaInv( 2, 2)= 4.844178E-01 JaInv( 2, 3)= 8.308070E-02 JaInv( 2, 4)= 2.650063E+00 JaInv( 3, 1)= 2.005050E+00 JaInv( 3, 2)= -4.236526E-01 JaInv( 3, 3)= 3.166529E-01 JaInv( 3, 4)= -1.896700E+00 JaInv( 4, 1)= -2.821679E-01 JaInv( 4, 2)= 4.730512E-01 JaInv( 4, 3)= -2.953063E-01 JaInv( 4, 4)= 1.513800E-01 change X_tmp2( 1)= 3.528162E-04 change X_tmp2( 2)= -5.339036E-04 change X_tmp2( 3)= 3.829762E-04 change X_tmp2( 4)= -1.238401E-04 X_next( 1)= 1.100001E+00 X_next( 2)= 1.199999E+00 X_next( 3)= 1.400001E+00 X_next( 4)= 1.500000E+00 X_next( 5)= 1.331003E+00 X_next( 6)= 1.727994E+00 X_next( 7)= 2.744006E+00 X_next( 8)= 3.374998E+00 X_next( 9)= 1.848001E+00 X_next(10)= 2.519998E+00 simeq_newton3 did not converge < 1.00000000000000E-06 test case3 solution X( 1)= 1.10000084984356E+00 X( 2)= 1.19999865021177E+00 X( 3)= 1.40000097469786E+00 X( 4)= 1.49999969397520E+00 sum of errors= 1.84320317497821E-06 end test_simeq_newton3