Метод Гаусса и Гаусса-Зейделя

#include "stdafx.h"
#include "iostream"
#include "math.h"
#include "conio.h"
using namespace std;
double E = 0.001;
int main()
{
    setlocale(LC_ALL, "RUSSIAN");
 
    int i, j, k;
    double a[5][5] = { 3, 1.2, -1, 2.2, -1,
        1.4, 5, 3, 2.6, 4,
        3, 3.4, 1, 1.6, 1,
        3.5, 2, 4, 3.8, 1.5,
        1, 1, 2.2, 1, 8 };
 
    double b[5] = { 5, 11, 2.5, 7, -8 };
 
    double e = 0.0001, m, x[5], s[5] = { 0 };
 
 
    cout << "Матрица А имеет вид: " << endl;
    for (i = 0; i < 5; i++)
    {
        for (j = 0; j < 5; j++)
        {
            printf("%5.1f", a[i][j]);
        }
        cout << endl;
    }
 
    cout << endl;
    cout << "Матрица В имеет вид: " << endl;
    for (j = 0; j < 5; j++)
    {
        printf("%5.1f", b[j]);
    }
 
    for (k = 0; k < 4; k++)
    {
        for (i = k + 1; i < 5; i++)
        {
            m = -a[i][k] / a[k][k];
            b[i] += m * b[k];
            for (j = k; j < 5; j++)
                a[i][j] = a[i][j] + m * a[k][j];
        }
    }
    x[4] = b[4] / a[4][4];
    cout << endl;
    cout << endl;
    cout << "Верхняя треугольная матрица: " << endl;
    for (i = 0; i < 5; i++)
    {
        for (j = 0; j < 5; j++)
        {
            printf("%5.1f", a[i][j]);
        }
        cout << endl;
    }
    cout << endl;
 
 
 
    for (k = 3; k >= 0; k--)
    {
        for (j = k + 1; j < 5; j++)
        {
            s[k] += a[k][j] * x[j];
        }
        x[k] = (b[k] - s[k]) / a[k][k];
    }
    cout << endl;
    cout << "Корни уравнения методом Гаусса: " << endl;
    for (k = 0; k < 5; k++)
    {
        cout << "x[" << k << "] = " << x[k] << ";" << endl;
    }
 
 
 
    double z[5] = { 0 }, d[5] = { 0 }; //вспомагательные
    double max = 0;
    int  l = 0;
    
    do
    {
        d[0] = z[0];
 
        z[0] = (b[0] - a[0][1] * z[1] - a[0][2] * z[2] - z[3] * a[0][3] - a[0][4] * z[4]) / a[0][0];
 
 
        d[1] = z[1];
 
        z[1] = (b[1] - a[1][0] * z[0] - a[1][2] * z[2] - z[3] * a[1][3] - a[1][4] * z[4]) / a[1][1];
 
 
        d[2] = z[2];
 
        z[2] = (b[2] - a[2][0] * z[0] - a[2][1] * z[1] - z[3] * a[2][3] - a[2][4] * z[4]) / a[2][2];
 
 
        d[3] = z[3];
 
        z[3] = (b[3] - a[3][0] * z[0] - a[3][1] * z[1] - z[2] * a[3][2] - a[3][4] * z[4]) / a[3][3];
 
 
        d[4] = z[4];
 
        z[4] = (b[4] - a[4][0] * z[0] - a[4][1] * z[1] - z[2] * a[4][2] - a[4][3] * z[3]) / a[4][4];
 
 
        for (i = 1; i < 5; i++)
        {
            if (z[i] > max)
            {
                max = z[i];
                l = i;
            }
        }
 
        j++;
    } while ((max - d[l]) >= E);
 
    cout << endl;
    cout << "\tРезультат метода Гаусса - Зейделя: " << endl;
    cout << endl;
    for (int i = 0; i < 5; i++)
    {
        cout << "   x[ " << i + 1 << "] = " << z[i] << endl;
    }
 
    _getche();
    return 0;
}