test_simeq_newton4.c running first A generated A[ 0][ 0]=0.000000 A[ 0][ 1]=0.000985 A[ 0][ 2]=0.041631 A[ 0][ 3]=0.176643 A[ 0][ 4]=0.364602 A[ 0][ 5]=0.091331 A[ 0][ 6]=0.092298 A[ 0][ 7]=0.487217 A[ 1][ 0]=0.526750 A[ 1][ 1]=0.454433 A[ 1][ 2]=0.233178 A[ 1][ 3]=0.831292 A[ 1][ 4]=0.931731 A[ 1][ 5]=0.568060 A[ 1][ 6]=0.556094 A[ 1][ 7]=0.050832 solve system of equations A X = Y for X = variables X0, X1, terms: X0 * X0 X1 * X1 X0 * X0 * X0 X1 * X1 * X1 X0 * X0 * X0 * X0 X1 * X1 * X1 * X1 X_soln[ 0]=1.100000 X_soln[ 1]=2.300000 X_soln[ 2]=1.210000 X_soln[ 3]=5.290000 X_soln[ 4]=1.331000 X_soln[ 5]=12.167000 X_soln[ 6]=1.464100 X_soln[ 7]=27.984100 Y[ 0]=16.353053 Y[ 1]=16.692681 simeq_newton4 running 2 equations, 2 variables, with 6 nonlinear terms X[ 0]=1.000000 X[ 1]=1.000000 X[ 2]=1.000000 X[ 3]=1.000000 X[ 4]=1.000000 X[ 5]=1.000000 X[ 6]=1.000000 X[ 7]=1.000000 residual X_resid[ 0]=-15.0983 residual X_resid[ 1]=-12.5403 simeq_newton4 itr=0, prev=1e+12, residual=27.6387 Ja computed Ja[ 0][ 0]=1.546259 Ja[ 0][ 1]=2.577131 Ja[ 1][ 0]=6.012679 Ja[ 1][ 1]=4.024523 Ja inverted invJa[ 0][ 0]=-0.434028 invJa[ 0][ 1]=0.277933 invJa[ 1][ 0]=0.648442 invJa[ 1][ 1]=-0.166757 change X_tmp2[ 0]=3.067738 change X_tmp2[ 1]=-7.699206 X_next[ 0]=-2.067738 X_next[ 1]=8.699206 X_next[ 2]=4.275542 X_next[ 3]=75.676184 X_next[ 4]=-8.840702 X_next[ 5]=658.322710 X_next[ 6]=18.280259 X_next[ 7]=5726.884836 residual X_resid[ 0]=2846.03 residual X_resid[ 1]=717.081 simeq_newton4 itr=1, prev=27.6387, residual=3563.11 Ja computed Ja[ 0][ 0]=1.240549 Ja[ 0][ 1]=1306.793597 Ja[ 1][ 0]=-8.151640 Ja[ 1][ 1]=277.738556 Ja inverted invJa[ 0][ 0]=0.025256 invJa[ 0][ 1]=-0.118831 invJa[ 1][ 0]=0.000741 invJa[ 1][ 1]=0.000113 change X_tmp2[ 0]=-13.333120 change X_tmp2[ 1]=2.190528 X_next[ 0]=11.265382 X_next[ 1]=6.508678 X_next[ 2]=126.908824 X_next[ 3]=42.362895 X_next[ 4]=1429.676332 X_next[ 5]=275.726456 X_next[ 6]=16105.849506 X_next[ 7]=1794.614831 residual X_resid[ 0]=2903.76 residual X_resid[ 1]=10593.3 simeq_newton4 itr=2, prev=3563.11, residual=13497.1 Ja computed Ja[ 0][ 0]=667.574755 Ja[ 0][ 1]=551.262208 Ja[ 1][ 0]=3540.654899 Ja[ 1][ 1]=139.532418 Ja inverted invJa[ 0][ 0]=-0.000075 invJa[ 0][ 1]=0.000297 invJa[ 1][ 0]=0.001905 invJa[ 1][ 1]=-0.000359 change X_tmp2[ 0]=2.923859 change X_tmp2[ 1]=1.726711 X_next[ 0]=8.341523 X_next[ 1]=4.781968 X_next[ 2]=69.581007 X_next[ 3]=22.867215 X_next[ 4]=580.411577 X_next[ 5]=109.350281 X_next[ 6]=4841.516560 X_next[ 7]=522.909506 residual X_resid[ 0]=913.825 residual X_resid[ 1]=3346.93 simeq_newton4 itr=3, prev=13497.1, residual=4260.76 Ja computed Ja[ 0][ 0]=291.085200 Ja[ 0][ 1]=221.065175 Ja[ 1][ 0]=1489.963675 Ja[ 1][ 1]=69.608614 Ja inverted invJa[ 0][ 0]=-0.000225 invJa[ 0][ 1]=0.000715 invJa[ 1][ 0]=0.004820 invJa[ 1][ 1]=-0.000942 change X_tmp2[ 0]=2.187781 change X_tmp2[ 1]=1.252999 X_next[ 0]=6.153742 X_next[ 1]=3.528968 X_next[ 2]=37.868539 X_next[ 3]=12.453618 X_next[ 4]=233.033215 X_next[ 5]=43.948427 X_next[ 6]=1434.026251 X_next[ 7]=155.092612 residual X_resid[ 0]=284.326 residual X_resid[ 1]=1054.76 simeq_newton4 itr=4, prev=4260.76, residual=1339.09 Ja computed Ja[ 0][ 0]=127.966907 Ja[ 0][ 1]=90.309630 Ja[ 1][ 0]=627.600320 Ja[ 1][ 1]=36.480762 Ja inverted invJa[ 0][ 0]=-0.000701 invJa[ 0][ 1]=0.001736 invJa[ 1][ 0]=0.012067 invJa[ 1][ 1]=-0.002460 change X_tmp2[ 0]=1.632046 change X_tmp2[ 1]=0.835771 X_next[ 0]=4.521695 X_next[ 1]=2.693198 X_next[ 2]=20.445729 X_next[ 3]=7.253314 X_next[ 4]=92.449361 X_next[ 5]=19.534609 X_next[ 6]=418.027850 X_next[ 7]=52.610565 residual X_resid[ 0]=85.4891 residual X_resid[ 1]=330.082 simeq_newton4 itr=5, prev=1339.09, residual=415.571 Ja computed Ja[ 0][ 0]=56.871596 Ja[ 0][ 1]=41.010194 Ja[ 1][ 0]=265.427525 Ja[ 1][ 1]=21.264970 Ja inverted invJa[ 0][ 0]=-0.002198 invJa[ 0][ 1]=0.004238 invJa[ 1][ 0]=0.027432 invJa[ 1][ 1]=-0.005878 change X_tmp2[ 0]=1.211139 change X_tmp2[ 1]=0.405015 X_next[ 0]=3.310557 X_next[ 1]=2.288183 X_next[ 2]=10.959787 X_next[ 3]=5.235779 X_next[ 4]=36.282997 X_next[ 5]=11.980419 X_next[ 6]=120.116926 X_next[ 7]=27.413384 residual X_resid[ 0]=23.7962 residual X_resid[ 1]=101.8 simeq_newton4 itr=6, prev=415.571, residual=125.597 Ja computed Ja[ 0][ 0]=25.658874 Ja[ 0][ 1]=25.592192 Ja[ 1][ 0]=113.412463 Ja[ 1][ 1]=15.617383 Ja inverted invJa[ 0][ 0]=-0.006243 invJa[ 0][ 1]=0.010230 invJa[ 1][ 0]=0.045333 invJa[ 1][ 1]=-0.010256 change X_tmp2[ 0]=0.892841 change X_tmp2[ 1]=0.034654 X_next[ 0]=2.417716 X_next[ 1]=2.253529 X_next[ 2]=5.845352 X_next[ 3]=5.078391 X_next[ 4]=14.132404 X_next[ 5]=11.444298 X_next[ 6]=34.168144 X_next[ 7]=25.790052 residual X_resid[ 0]=6.70649 residual X_resid[ 1]=31.1699 simeq_newton4 itr=7, prev=125.597, residual=37.8764 Ja computed Ja[ 0][ 0]=11.812540 Ja[ 0][ 1]=24.491998 Ja[ 1][ 0]=49.428964 Ja[ 1][ 1]=15.182541 Ja inverted invJa[ 0][ 0]=-0.014722 invJa[ 0][ 1]=0.023749 invJa[ 1][ 0]=0.047930 invJa[ 1][ 1]=-0.011454 change X_tmp2[ 0]=0.641531 change X_tmp2[ 1]=-0.035588 X_next[ 0]=1.776186 X_next[ 1]=2.289116 X_next[ 2]=3.154836 X_next[ 3]=5.240053 X_next[ 4]=5.603574 X_next[ 5]=11.995089 X_next[ 6]=9.952988 X_next[ 7]=27.458153 residual X_resid[ 0]=2.14148 residual X_resid[ 1]=9.34033 simeq_newton4 itr=8, prev=37.8764, residual=11.4818 Ja computed Ja[ 0][ 0]=5.667456 Ja[ 0][ 1]=25.622285 Ja[ 1][ 0]=22.637929 Ja[ 1][ 1]=15.629201 Ja inverted invJa[ 0][ 0]=-0.031802 invJa[ 0][ 1]=0.052135 invJa[ 1][ 0]=0.046063 invJa[ 1][ 1]=-0.011532 change X_tmp2[ 0]=0.418858 change X_tmp2[ 1]=-0.009070 X_next[ 0]=1.357327 X_next[ 1]=2.298186 X_next[ 2]=1.842338 X_next[ 3]=5.281658 X_next[ 4]=2.500656 X_next[ 5]=12.138232 X_next[ 6]=3.394209 X_next[ 7]=27.895912 residual X_resid[ 0]=0.583859 residual X_resid[ 1]=2.41754 simeq_newton4 itr=9, prev=11.4818, residual=3.0014 Ja computed Ja[ 0][ 0]=3.051394 Ja[ 0][ 1]=25.915854 Ja[ 1][ 0]=11.871843 Ja[ 1][ 1]=15.744288 Ja inverted invJa[ 0][ 0]=-0.060642 invJa[ 0][ 1]=0.099820 invJa[ 1][ 0]=0.045727 invJa[ 1][ 1]=-0.011753 change X_tmp2[ 0]=0.205911 change X_tmp2[ 1]=-0.001715 X_next[ 0]=1.151416 X_next[ 1]=2.299901 X_next[ 2]=1.325759 X_next[ 3]=5.289546 X_next[ 4]=1.526500 X_next[ 5]=12.165433 X_next[ 6]=1.757637 X_next[ 7]=27.979296 residual X_resid[ 0]=0.100628 residual X_resid[ 1]=0.397908 simeq_newton4 itr=10, prev=3.0014, residual=0.498535 Ja computed Ja[ 0][ 0]=2.109563 Ja[ 0][ 1]=25.971634 Ja[ 1][ 0]=8.164988 Ja[ 1][ 1]=15.766113 Ja inverted invJa[ 0][ 0]=-0.088178 invJa[ 0][ 1]=0.145256 invJa[ 1][ 0]=0.045666 invJa[ 1][ 1]=-0.011799 change X_tmp2[ 0]=0.048925 change X_tmp2[ 1]=-0.000099 X_next[ 0]=1.102491 X_next[ 1]=2.300001 X_next[ 2]=1.215486 X_next[ 3]=5.290003 X_next[ 4]=1.340061 X_next[ 5]=12.167012 X_next[ 6]=1.477405 X_next[ 7]=27.984137 residual X_resid[ 0]=0.00477998 residual X_resid[ 1]=0.0184448 simeq_newton4 itr=11, prev=0.498535, residual=0.0232248 Ja computed Ja[ 0][ 0]=1.916040 Ja[ 0][ 1]=25.974871 Ja[ 1][ 0]=7.419225 Ja[ 1][ 1]=15.767379 Ja inverted invJa[ 0][ 0]=-0.097029 invJa[ 0][ 1]=0.159843 invJa[ 1][ 0]=0.045656 invJa[ 1][ 1]=-0.011791 change X_tmp2[ 0]=0.002484 change X_tmp2[ 1]=0.000001 X_next[ 0]=1.100006 X_next[ 1]=2.300000 X_next[ 2]=1.210014 X_next[ 3]=5.290000 X_next[ 4]=1.331022 X_next[ 5]=12.167000 X_next[ 6]=1.464133 X_next[ 7]=27.984100 residual X_resid[ 0]=1.18437e-05 residual X_resid[ 1]=4.5443e-05 simeq_newton4 itr=12, prev=0.0232248, residual=5.72867e-05 Ja computed Ja[ 0][ 0]=1.906511 Ja[ 0][ 1]=25.974846 Ja[ 1][ 0]=7.382665 Ja[ 1][ 1]=15.767369 Ja inverted invJa[ 0][ 0]=-0.097508 invJa[ 0][ 1]=0.160633 invJa[ 1][ 0]=0.045656 invJa[ 1][ 1]=-0.011790 change X_tmp2[ 0]=0.000006 change X_tmp2[ 1]=0.000000 X_next[ 0]=1.100000 X_next[ 1]=2.300000 X_next[ 2]=1.210000 X_next[ 3]=5.290000 X_next[ 4]=1.331000 X_next[ 5]=12.167000 X_next[ 6]=1.464100 X_next[ 7]=27.984100 residual X_resid[ 0]=7.23013e-11 residual X_resid[ 1]=2.77339e-10 simeq_newton4 itr=13, prev=5.72867e-05, residual=3.4964e-10 simeq_newton4 finished Returned solution X[ 0]=1.100000, err=3.75e-11 X[ 1]=2.300000, err=3.10862e-14 first test finished