Метод простых итераций (для трехмерной матрицы) - C++ - Ответ 2056294

#include "stdafx.h"
#include "iostream"
#include "conio.h"
using namespace std;
 
 
int _tmain(int argc, _TCHAR* argv[])
{
    int i, j, n = 3;
    float a[3][3], al[3][3], b[3], bt[3], x0[3], x1[3];
    float e = 0.001, k = 0;
 
    cout << "Massiv a: " << endl;
 
    for (i = 0; i < n; i++)
        for (j = 0; j < n; j++)
        {
            cout << "a[" << i << "]" << "[" << j << "] = ";
            cin >> a[i][j];
        }
        cout << endl;
    
    cout << "Massiv b: " << endl;
 
    for (i = 0; i < n; i++)
    {
        cout << "b[" << i <<"] = ";
        cin >> b[i];
    }
    cout << endl;
 
    for (i = 0; i < n; i++)
        for (j = 0; j < n; j++)
        {
            if (i = j)
                al[i][j] = 0;
            else
                al[i][j] = (-1) * (a[i][j] / a[i][i]);
        }
 
        for (i = 0; i < n; i++)
        {
            bt[i] = b[i] / a[i][i];
            x0[i] = bt[i];
        }
 
        for (i = 0; i < n; i++)
            for (j = 0; j < n; j++)
            {
                x1[i] = bt[i] + (al[i][j] * x0[i]);
                x0[i] = x1[i];
            }
            
        for (i = 0; i < n; i++)
            k = k + sqrt(x1[i] - x0[i]);
 
    for (i = 0; i < n; i++)
    {
        for (j = 0; j < n; j++)
        {
            do
            {
                x1[i] = bt[i] + (al[i][j] * x0[i]);
                x0[i] = x1[i];
            }
            while (k > e);
        }
    
        cout << "x1["<< i << "] = " << x1[i] << endl;
    }
    
    
    getch();
    return 0;
}