Q = -0.4888 0.6382 -0.1493 0.5743 -0.0396 -0.5844 -0.7385 -0.1753 0.2815 0.0553 -0.3614 0.1276 0.4569 -0.2788 0.7527 -0.3123 0.1696 -0.6880 -0.6327 0.0045 -0.4374 0.0469 0.5148 -0.3358 -0.6548 D = 0.0574 0 0 0 0 0 0.2910 0 0 0 0 0 0.6761 0 0 0 0 0 0.9522 0 0 0 0 0 0.0403 A = 0.4614 0.0508 -0.1659 -0.2364 -0.2136 0.0508 0.2747 -0.1425 -0.1140 -0.1479 -0.1659 -0.1425 0.2502 -0.0316 0.2391 -0.2364 -0.1140 -0.0316 0.7152 -0.0270 -0.2136 -0.1479 0.2391 -0.0270 0.3154 z = 0.8216 0.2236 0.9997 0.2605 0.2303 b = 0.1139 -0.1031 0.1287 -0.0712 0.0961 x0 = 0.8217 0.2237 0.9998 0.2606 0.2303 Running the method of Jacobi... norm(dx) = 4.3533e-05 norm(dx) = 5.2145e-05 norm(dx) = 8.4970e-05 norm(dx) = 1.3410e-04 norm(dx) = 2.2022e-04 norm(dx) = 3.5636e-04 norm(dx) = 5.8190e-04 norm(dx) = 9.4612e-04 norm(dx) = 1.5418e-03 norm(dx) = 2.5097e-03 norm(dx) = 4.0877e-03 norm(dx) = 6.6558e-03 norm(dx) = 1.0839e-02 norm(dx) = 1.7650e-02 norm(dx) = 2.8742e-02 norm(dx) = 4.6803e-02 norm(dx) = 7.6215e-02 norm(dx) = 1.2411e-01 norm(dx) = 2.0210e-01 norm(dx) = 3.2911e-01 norm(dx) = 5.3592e-01 norm(dx) = 8.7270e-01 norm(dx) = 1.4211e+00 norm(dx) = 2.3142e+00 norm(dx) = 3.7684e+00 norm(dx) = 6.1366e+00 norm(dx) = 9.9929e+00 norm(dx) = 1.6273e+01 norm(dx) = 2.6498e+01 norm(dx) = 4.3150e+01 norm(dx) = 7.0267e+01 norm(dx) = 1.1442e+02 norm(dx) = 1.8633e+02 norm(dx) = 3.0342e+02 norm(dx) = 4.9409e+02 norm(dx) = 8.0459e+02 norm(dx) = 1.3102e+03 norm(dx) = 2.1336e+03 norm(dx) = 3.4743e+03 norm(dx) = 5.6576e+03 norm(dx) = 9.2129e+03 norm(dx) = 1.5002e+04 norm(dx) = 2.4430e+04 norm(dx) = 3.9783e+04 norm(dx) = 6.4782e+04 norm(dx) = 1.0549e+05 norm(dx) = 1.7179e+05 norm(dx) = 2.7974e+05 norm(dx) = 4.5553e+05 norm(dx) = 7.4179e+05 Failed to reached accuracy requirement in 50 steps. x = 1.0e+05 * -1.5999 -2.0154 2.8145 0.3867 2.5353 jdx = 1.0e+05 * -2.5824 -3.2530 4.5429 0.6242 4.0923 Running the method of Gauss-Seidel... norm(dx) = 3.6645e-05 norm(dx) = 2.7888e-05 norm(dx) = 2.0578e-05 norm(dx) = 1.5186e-05 norm(dx) = 1.1254e-05 norm(dx) = 8.3932e-06 norm(dx) = 6.3109e-06 norm(dx) = 4.7886e-06 norm(dx) = 3.6669e-06 norm(dx) = 2.8313e-06 norm(dx) = 2.2009e-06 norm(dx) = 1.7190e-06 norm(dx) = 1.3463e-06 norm(dx) = 1.0551e-06 norm(dx) = 8.2618e-07 norm(dx) = 6.4549e-07 norm(dx) = 5.0269e-07 norm(dx) = 3.8992e-07 norm(dx) = 3.0109e-07 norm(dx) = 2.3136e-07 norm(dx) = 1.7687e-07 norm(dx) = 1.3451e-07 norm(dx) = 1.0176e-07 norm(dx) = 7.6578e-08 norm(dx) = 5.7340e-08 norm(dx) = 4.2731e-08 norm(dx) = 3.1707e-08 norm(dx) = 2.3439e-08 norm(dx) = 1.7277e-08 norm(dx) = 1.2713e-08 norm(dx) = 9.3514e-09 Succeeded in 31 steps x = 0.8216 0.2236 0.9997 0.2605 0.2303 gdx = 1.0e-08 * 0.4003 0.5136 0.6236 0.2407 0.0597 Running the method of Jacobi... norm(dx) = 4.3533e-05 norm(dx) = 5.2145e-05 norm(dx) = 8.4970e-05 norm(dx) = 1.3410e-04 norm(dx) = 2.2022e-04 norm(dx) = 3.5636e-04 norm(dx) = 5.8190e-04 norm(dx) = 9.4612e-04 norm(dx) = 1.5418e-03 norm(dx) = 2.5097e-03 norm(dx) = 4.0877e-03 norm(dx) = 6.6558e-03 norm(dx) = 1.0839e-02 norm(dx) = 1.7650e-02 norm(dx) = 2.8742e-02 norm(dx) = 4.6803e-02 norm(dx) = 7.6215e-02 norm(dx) = 1.2411e-01 norm(dx) = 2.0210e-01 norm(dx) = 3.2911e-01 norm(dx) = 5.3592e-01 norm(dx) = 8.7270e-01 norm(dx) = 1.4211e+00 norm(dx) = 2.3142e+00 norm(dx) = 3.7684e+00 norm(dx) = 6.1366e+00 norm(dx) = 9.9929e+00 norm(dx) = 1.6273e+01 norm(dx) = 2.6498e+01 norm(dx) = 4.3150e+01 norm(dx) = 7.0267e+01 norm(dx) = 1.1442e+02 norm(dx) = 1.8633e+02 norm(dx) = 3.0342e+02 norm(dx) = 4.9409e+02 norm(dx) = 8.0459e+02 norm(dx) = 1.3102e+03 norm(dx) = 2.1336e+03 norm(dx) = 3.4743e+03 norm(dx) = 5.6576e+03 norm(dx) = 9.2129e+03 norm(dx) = 1.5002e+04 norm(dx) = 2.4430e+04 norm(dx) = 3.9783e+04 norm(dx) = 6.4782e+04 norm(dx) = 1.0549e+05 norm(dx) = 1.7179e+05 norm(dx) = 2.7974e+05 norm(dx) = 4.5553e+05 norm(dx) = 7.4179e+05 norm(dx) = 1.2079e+06 norm(dx) = 1.9670e+06 norm(dx) = 3.2031e+06 norm(dx) = 5.2160e+06 norm(dx) = 8.4939e+06 norm(dx) = 1.3832e+07 norm(dx) = 2.2524e+07 norm(dx) = 3.6678e+07 norm(dx) = 5.9726e+07 norm(dx) = 9.7259e+07 norm(dx) = 1.5838e+08 norm(dx) = 2.5791e+08 norm(dx) = 4.1998e+08 norm(dx) = 6.8390e+08 norm(dx) = 1.1137e+09 norm(dx) = 1.8135e+09 norm(dx) = 2.9531e+09 norm(dx) = 4.8089e+09 norm(dx) = 7.8309e+09 norm(dx) = 1.2752e+10 norm(dx) = 2.0766e+10 norm(dx) = 3.3815e+10 norm(dx) = 5.5065e+10 norm(dx) = 8.9668e+10 norm(dx) = 1.4602e+11 norm(dx) = 2.3778e+11 norm(dx) = 3.8720e+11 norm(dx) = 6.3052e+11 norm(dx) = 1.0267e+12 norm(dx) = 1.6720e+12 norm(dx) = 2.7227e+12 norm(dx) = 4.4336e+12 norm(dx) = 7.2197e+12 norm(dx) = 1.1757e+13 norm(dx) = 1.9145e+13 norm(dx) = 3.1176e+13 norm(dx) = 5.0767e+13 norm(dx) = 8.2670e+13 norm(dx) = 1.3462e+14 norm(dx) = 2.1922e+14 norm(dx) = 3.5698e+14 norm(dx) = 5.8131e+14 norm(dx) = 9.4661e+14 norm(dx) = 1.5415e+15 norm(dx) = 2.5102e+15 norm(dx) = 4.0876e+15 norm(dx) = 6.6563e+15 norm(dx) = 1.0839e+16 norm(dx) = 1.7651e+16 norm(dx) = 2.8743e+16 norm(dx) = 4.6805e+16 norm(dx) = 7.6217e+16 norm(dx) = 1.2411e+17 norm(dx) = 2.0211e+17 norm(dx) = 3.2912e+17 norm(dx) = 5.3594e+17 norm(dx) = 8.7273e+17 norm(dx) = 1.4212e+18 norm(dx) = 2.3142e+18 norm(dx) = 3.7685e+18 norm(dx) = 6.1367e+18 norm(dx) = 9.9932e+18 norm(dx) = 1.6273e+19 norm(dx) = 2.6499e+19 norm(dx) = 4.3152e+19 norm(dx) = 7.0269e+19 norm(dx) = 1.1443e+20 norm(dx) = 1.8633e+20 norm(dx) = 3.0343e+20 norm(dx) = 4.9411e+20 norm(dx) = 8.0461e+20 norm(dx) = 1.3102e+21 norm(dx) = 2.1336e+21 norm(dx) = 3.4744e+21 norm(dx) = 5.6578e+21 norm(dx) = 9.2132e+21 norm(dx) = 1.5003e+22 norm(dx) = 2.4431e+22 norm(dx) = 3.9784e+22 norm(dx) = 6.4784e+22 norm(dx) = 1.0550e+23 norm(dx) = 1.7179e+23 norm(dx) = 2.7975e+23 norm(dx) = 4.5554e+23 norm(dx) = 7.4181e+23 norm(dx) = 1.2080e+24 norm(dx) = 1.9671e+24 norm(dx) = 3.2032e+24 norm(dx) = 5.2162e+24 norm(dx) = 8.4941e+24 norm(dx) = 1.3832e+25 norm(dx) = 2.2524e+25 norm(dx) = 3.6679e+25 norm(dx) = 5.9728e+25 norm(dx) = 9.7262e+25 norm(dx) = 1.5838e+26 norm(dx) = 2.5791e+26 norm(dx) = 4.1999e+26 norm(dx) = 6.8392e+26 norm(dx) = 1.1137e+27 norm(dx) = 1.8136e+27 norm(dx) = 2.9532e+27 norm(dx) = 4.8091e+27 norm(dx) = 7.8312e+27 norm(dx) = 1.2752e+28 norm(dx) = 2.0766e+28 norm(dx) = 3.3816e+28 norm(dx) = 5.5066e+28 norm(dx) = 8.9671e+28 norm(dx) = 1.4602e+29 norm(dx) = 2.3778e+29 norm(dx) = 3.8721e+29 norm(dx) = 6.3054e+29 norm(dx) = 1.0268e+30 norm(dx) = 1.6720e+30 norm(dx) = 2.7227e+30 norm(dx) = 4.4337e+30 norm(dx) = 7.2200e+30 norm(dx) = 1.1757e+31 norm(dx) = 1.9145e+31 norm(dx) = 3.1177e+31 norm(dx) = 5.0768e+31 norm(dx) = 8.2672e+31 norm(dx) = 1.3462e+32 norm(dx) = 2.1922e+32 norm(dx) = 3.5699e+32 norm(dx) = 5.8132e+32 norm(dx) = 9.4664e+32 norm(dx) = 1.5415e+33 norm(dx) = 2.5102e+33 norm(dx) = 4.0877e+33 norm(dx) = 6.6564e+33 norm(dx) = 1.0839e+34 norm(dx) = 1.7651e+34 norm(dx) = 2.8743e+34 norm(dx) = 4.6806e+34 norm(dx) = 7.6220e+34 norm(dx) = 1.2412e+35 norm(dx) = 2.0211e+35 norm(dx) = 3.2913e+35 norm(dx) = 5.3595e+35 norm(dx) = 8.7275e+35 norm(dx) = 1.4212e+36 norm(dx) = 2.3143e+36 norm(dx) = 3.7686e+36 norm(dx) = 6.1369e+36 norm(dx) = 9.9935e+36 norm(dx) = 1.6273e+37 norm(dx) = 2.6500e+37 norm(dx) = 4.3153e+37 norm(dx) = 7.0271e+37 norm(dx) = 1.1443e+38 norm(dx) = 1.8634e+38 norm(dx) = 3.0344e+38 norm(dx) = 4.9412e+38 norm(dx) = 8.0464e+38 norm(dx) = 1.3103e+39 norm(dx) = 2.1337e+39 norm(dx) = 3.4745e+39 norm(dx) = 5.6579e+39 norm(dx) = 9.2135e+39 norm(dx) = 1.5003e+40 norm(dx) = 2.4432e+40 norm(dx) = 3.9785e+40 norm(dx) = 6.4786e+40 norm(dx) = 1.0550e+41 norm(dx) = 1.7180e+41 norm(dx) = 2.7975e+41 norm(dx) = 4.5556e+41 norm(dx) = 7.4183e+41 norm(dx) = 1.2080e+42 norm(dx) = 1.9671e+42 norm(dx) = 3.2033e+42 norm(dx) = 5.2163e+42 norm(dx) = 8.4944e+42 norm(dx) = 1.3832e+43 norm(dx) = 2.2525e+43 norm(dx) = 3.6680e+43 norm(dx) = 5.9730e+43 norm(dx) = 9.7265e+43 norm(dx) = 1.5839e+44 norm(dx) = 2.5792e+44 norm(dx) = 4.2000e+44 norm(dx) = 6.8394e+44 norm(dx) = 1.1137e+45 norm(dx) = 1.8136e+45 norm(dx) = 2.9533e+45 norm(dx) = 4.8092e+45 norm(dx) = 7.8314e+45 norm(dx) = 1.2753e+46 norm(dx) = 2.0767e+46 norm(dx) = 3.3817e+46 norm(dx) = 5.5068e+46 norm(dx) = 8.9673e+46 norm(dx) = 1.4603e+47 norm(dx) = 2.3779e+47 norm(dx) = 3.8722e+47 norm(dx) = 6.3056e+47 norm(dx) = 1.0268e+48 norm(dx) = 1.6721e+48 norm(dx) = 2.7228e+48 norm(dx) = 4.4339e+48 norm(dx) = 7.2202e+48 norm(dx) = 1.1757e+49 norm(dx) = 1.9146e+49 norm(dx) = 3.1178e+49 norm(dx) = 5.0770e+49 norm(dx) = 8.2675e+49 norm(dx) = 1.3463e+50 norm(dx) = 2.1923e+50 norm(dx) = 3.5700e+50 norm(dx) = 5.8134e+50 norm(dx) = 9.4666e+50 norm(dx) = 1.5416e+51 norm(dx) = 2.5103e+51 norm(dx) = 4.0878e+51 norm(dx) = 6.6566e+51 norm(dx) = 1.0840e+52 norm(dx) = 1.7652e+52 norm(dx) = 2.8744e+52 norm(dx) = 4.6807e+52 norm(dx) = 7.6222e+52 norm(dx) = 1.2412e+53 norm(dx) = 2.0212e+53 norm(dx) = 3.2914e+53 norm(dx) = 5.3597e+53 norm(dx) = 8.7278e+53 norm(dx) = 1.4212e+54 norm(dx) = 2.3144e+54 norm(dx) = 3.7688e+54 norm(dx) = 6.1371e+54 norm(dx) = 9.9937e+54 norm(dx) = 1.6274e+55 norm(dx) = 2.6501e+55 norm(dx) = 4.3154e+55 norm(dx) = 7.0273e+55 norm(dx) = 1.1443e+56 norm(dx) = 1.8634e+56 norm(dx) = 3.0345e+56 norm(dx) = 4.9414e+56 norm(dx) = 8.0466e+56 norm(dx) = 1.3103e+57 norm(dx) = 2.1337e+57 norm(dx) = 3.4746e+57 norm(dx) = 5.6581e+57 norm(dx) = 9.2137e+57 norm(dx) = 1.5004e+58 norm(dx) = 2.4432e+58 norm(dx) = 3.9786e+58 norm(dx) = 6.4788e+58 norm(dx) = 1.0550e+59 norm(dx) = 1.7180e+59 norm(dx) = 2.7976e+59 norm(dx) = 4.5557e+59 norm(dx) = 7.4186e+59 norm(dx) = 1.2080e+60 norm(dx) = 1.9672e+60 norm(dx) = 3.2034e+60 norm(dx) = 5.2165e+60 norm(dx) = 8.4946e+60 norm(dx) = 1.3833e+61 norm(dx) = 2.2525e+61 norm(dx) = 3.6681e+61 norm(dx) = 5.9731e+61 norm(dx) = 9.7268e+61 norm(dx) = 1.5839e+62 norm(dx) = 2.5793e+62 norm(dx) = 4.2001e+62 norm(dx) = 6.8396e+62 norm(dx) = 1.1138e+63 norm(dx) = 1.8137e+63 norm(dx) = 2.9534e+63 norm(dx) = 4.8094e+63 norm(dx) = 7.8316e+63 norm(dx) = 1.2753e+64 norm(dx) = 2.0767e+64 norm(dx) = 3.3818e+64 norm(dx) = 5.5070e+64 norm(dx) = 8.9676e+64 norm(dx) = 1.4603e+65 norm(dx) = 2.3780e+65 norm(dx) = 3.8723e+65 norm(dx) = 6.3057e+65 norm(dx) = 1.0268e+66 norm(dx) = 1.6721e+66 norm(dx) = 2.7229e+66 norm(dx) = 4.4340e+66 norm(dx) = 7.2204e+66 norm(dx) = 1.1758e+67 norm(dx) = 1.9147e+67 norm(dx) = 3.1178e+67 norm(dx) = 5.0771e+67 norm(dx) = 8.2677e+67 norm(dx) = 1.3463e+68 norm(dx) = 2.1924e+68 norm(dx) = 3.5701e+68 norm(dx) = 5.8136e+68 norm(dx) = 9.4669e+68 norm(dx) = 1.5416e+69 norm(dx) = 2.5104e+69 norm(dx) = 4.0879e+69 norm(dx) = 6.6568e+69 norm(dx) = 1.0840e+70 norm(dx) = 1.7652e+70 norm(dx) = 2.8745e+70 norm(dx) = 4.6809e+70 norm(dx) = 7.6224e+70 norm(dx) = 1.2412e+71 norm(dx) = 2.0213e+71 norm(dx) = 3.2914e+71 norm(dx) = 5.3598e+71 norm(dx) = 8.7280e+71 norm(dx) = 1.4213e+72 norm(dx) = 2.3144e+72 norm(dx) = 3.7689e+72 norm(dx) = 6.1373e+72 norm(dx) = 9.9940e+72 norm(dx) = 1.6274e+73 norm(dx) = 2.6501e+73 norm(dx) = 4.3155e+73 norm(dx) = 7.0275e+73 norm(dx) = 1.1444e+74 norm(dx) = 1.8635e+74 norm(dx) = 3.0346e+74 norm(dx) = 4.9415e+74 norm(dx) = 8.0468e+74 norm(dx) = 1.3104e+75 norm(dx) = 2.1338e+75 norm(dx) = 3.4747e+75 norm(dx) = 5.6583e+75 norm(dx) = 9.2140e+75 norm(dx) = 1.5004e+76 norm(dx) = 2.4433e+76 norm(dx) = 3.9787e+76 norm(dx) = 6.4790e+76 norm(dx) = 1.0550e+77 norm(dx) = 1.7181e+77 norm(dx) = 2.7977e+77 norm(dx) = 4.5558e+77 norm(dx) = 7.4188e+77 norm(dx) = 1.2081e+78 norm(dx) = 1.9673e+78 norm(dx) = 3.2035e+78 norm(dx) = 5.2167e+78 norm(dx) = 8.4949e+78 norm(dx) = 1.3833e+79 norm(dx) = 2.2526e+79 norm(dx) = 3.6682e+79 norm(dx) = 5.9733e+79 norm(dx) = 9.7270e+79 norm(dx) = 1.5840e+80 norm(dx) = 2.5794e+80 norm(dx) = 4.2003e+80 norm(dx) = 6.8398e+80 norm(dx) = 1.1138e+81 norm(dx) = 1.8137e+81 norm(dx) = 2.9535e+81 norm(dx) = 4.8095e+81 norm(dx) = 7.8319e+81 norm(dx) = 1.2754e+82 norm(dx) = 2.0768e+82 norm(dx) = 3.3819e+82 norm(dx) = 5.5071e+82 norm(dx) = 8.9679e+82 norm(dx) = 1.4603e+83 norm(dx) = 2.3780e+83 norm(dx) = 3.8724e+83 norm(dx) = 6.3059e+83 norm(dx) = 1.0269e+84 norm(dx) = 1.6722e+84 norm(dx) = 2.7230e+84 norm(dx) = 4.4341e+84 norm(dx) = 7.2206e+84 norm(dx) = 1.1758e+85 norm(dx) = 1.9147e+85 norm(dx) = 3.1179e+85 norm(dx) = 5.0773e+85 norm(dx) = 8.2679e+85 norm(dx) = 1.3464e+86 norm(dx) = 2.1924e+86 norm(dx) = 3.5702e+86 norm(dx) = 5.8137e+86 norm(dx) = 9.4672e+86 norm(dx) = 1.5417e+87 norm(dx) = 2.5104e+87 norm(dx) = 4.0880e+87 norm(dx) = 6.6570e+87 norm(dx) = 1.0840e+88 norm(dx) = 1.7653e+88 norm(dx) = 2.8746e+88 norm(dx) = 4.6810e+88 norm(dx) = 7.6226e+88 norm(dx) = 1.2413e+89 norm(dx) = 2.0213e+89 norm(dx) = 3.2915e+89 norm(dx) = 5.3600e+89 norm(dx) = 8.7283e+89 norm(dx) = 1.4213e+90 norm(dx) = 2.3145e+90 norm(dx) = 3.7690e+90 norm(dx) = 6.1375e+90 norm(dx) = 9.9943e+90 norm(dx) = 1.6275e+91 norm(dx) = 2.6502e+91 norm(dx) = 4.3157e+91 norm(dx) = 7.0277e+91 norm(dx) = 1.1444e+92 norm(dx) = 1.8636e+92 norm(dx) = 3.0346e+92 norm(dx) = 4.9417e+92 norm(dx) = 8.0471e+92 norm(dx) = 1.3104e+93 norm(dx) = 2.1339e+93 norm(dx) = 3.4748e+93 norm(dx) = 5.6584e+93 norm(dx) = 9.2143e+93 norm(dx) = 1.5005e+94 norm(dx) = 2.4434e+94 norm(dx) = 3.9788e+94 norm(dx) = 6.4792e+94 norm(dx) = 1.0551e+95 norm(dx) = 1.7181e+95 norm(dx) = 2.7978e+95 norm(dx) = 4.5560e+95 norm(dx) = 7.4190e+95 norm(dx) = 1.2081e+96 norm(dx) = 1.9673e+96 norm(dx) = 3.2036e+96 norm(dx) = 5.2168e+96 norm(dx) = 8.4951e+96 norm(dx) = 1.3834e+97 norm(dx) = 2.2527e+97 norm(dx) = 3.6683e+97 norm(dx) = 5.9735e+97 norm(dx) = 9.7273e+97 norm(dx) = 1.5840e+98 norm(dx) = 2.5794e+98 norm(dx) = 4.2004e+98 norm(dx) = 6.8399e+98 norm(dx) = 1.1138e+99 norm(dx) = 1.8138e+99 norm(dx) = 2.9536e+99 norm(dx) = 4.8096e+99 norm(dx) = 7.8321e+99 norm(dx) = 1.2754e+100 norm(dx) = 2.0769e+100 norm(dx) = 3.3820e+100 norm(dx) = 5.5073e+100 norm(dx) = 8.9681e+100 norm(dx) = 1.4604e+101 Failed to reached accuracy requirement in 500 steps. x = 1.0e+100 * -3.1498 -3.9677 5.5410 0.7614 4.9914 jdx = 1.0e+100 * -5.0841 -6.4042 8.9438 1.2290 8.0566 found case diverging for Jacobi, Gauss-Seidel converges