数値解析 ヤコビ法・ガウス-ザイデル法実行例
連立一次方程式 A x = b の反復解法 Jacobi 法と Gauss-Seidel 法の比較 A : 3.0000000 -2.0000000 1.0000000 3.0000000 b : 1.0000000 4.0000000 ----- Jacobi 法 ----- 0-th step : 0.0000000 0.0000000 1-th step : 0.3333333 1.3333333 2-th step : 1.2222222 1.2222222 3-th step : 1.1481481 0.9259259 4-th step : 0.9506173 0.9506173 5-th step : 0.9670782 1.0164609 6-th step : 1.0109739 1.0109739 7-th step : 1.0073160 0.9963420 8-th step : 0.9975613 0.9975613 9-th step : 0.9983742 1.0008129 10-th step : 1.0005419 1.0005419 11-th step : 1.0003613 0.9998194 12-th step : 0.9998796 0.9998796 13-th step : 0.9999197 1.0000401 14-th step : 1.0000268 1.0000268 15-th step : 1.0000178 0.9999911 16-th step : 0.9999941 0.9999941 17-th step : 0.9999960 1.0000020 18-th step : 1.0000013 1.0000013 19-th step : 1.0000009 0.9999996 20-th step : 0.9999997 0.9999997 21-th step : 0.9999998 1.0000001 22-th step : 1.0000001 1.0000001 23-th step : 1.0000000 1.0000000 24-th step : 1.0000000 1.0000000 24-th step : 1.0000000 1.0000000 検算 Ax-b : -0.0000000 -0.0000001 ----- Gauss-Seidel 法 ----- 0-th step : 0.0000000 0.0000000 1-th step : 0.3333333 1.2222222 2-th step : 1.1481481 0.9506173 3-th step : 0.9670782 1.0109739 4-th step : 1.0073160 0.9975613 5-th step : 0.9983742 1.0005419 6-th step : 1.0003613 0.9998796 7-th step : 0.9999197 1.0000268 8-th step : 1.0000178 0.9999941 9-th step : 0.9999960 1.0000013 10-th step : 1.0000009 0.9999997 11-th step : 0.9999998 1.0000001 12-th step : 1.0000000 1.0000000 13-th step : 1.0000000 1.0000000 13-th step : 1.0000000 1.0000000 検算 Ax-b : -0.0000000 0.0000000