Метод Гаусса

[Maray]

Добрый день!

Подскажите, пожалуйста, написала программу, которая реализует метод Гаусса. такой вот вопрос, я использую цикл не до конца, то есть, цикл у меня проходит один круг и принудительно завершается. Программа работает, но вот можно ли так реализовывать?

#include<iostream>
#include<math.h>
#include <locale>
#define n 3
 
using namespace std;
 
int Vivod(double A[n][n], double B[n]){
    int i, j;
    for (i = 0; i < n; i++){
 
        for (j = 0; j < n; j++){
            cout << A[i][j] << "X" << i + 1 << "+" << A[i][j + 1] << "X" << i + 2 << "+" << A[i][j + 2] << "X" << i + 3 << "=" << B[i] << endl;
            cout << A[i + 1][j] << "X" << i + 1 << "+" << A[i + 1][j + 1] << "X" << i + 2 << "+" << A[i + 1][j + 2] << "X" << i + 3 << "=" << B[i + 1] << endl;
            cout << A[i + 2][j] << "X" << i + 1 << "+" << A[i + 2][j + 1] << "X" << i + 2 << "+" << A[i + 2][j + 2] << "X" << i + 3 << "=" << B[i + 2] << endl;
            cout << endl;
            return 0;
        }
 
    }
 
}
 
int main(){
    setlocale(LC_ALL, "Russian");
    double A[n][n], B[n], X[n], k = 0;
    int i, j;
    bool flag1, flag2, flag3;
 
    flag1 = false;
    flag2 = false;
    flag3 = false;
    cout << "Введите матрицу А" << endl;
    for (i = 0; i < n; i++){
 
        for (j = 0; j < n; j++){
            cout << "A[" << i + 1 << "][" << j + 1 << "]= ";
            cin >> A[i][j];
        }
 
    }
    cout << endl;
    cout << "Введите матрицу B" << endl;
    for (i = 0; i < n; i++){
        cout << "B[" << i + 1 << "]= ";
        cin >> B[i];
    }
    cout << endl;
    Vivod(A, B);
 
    if (A[1][0] == 0)
    {
    }
    else
    {
        k = A[1][0] / A[0][0];
        A[1][0] = A[1][0] - k*A[0][0];
        A[1][1] = A[1][1] - k*A[0][1];
        A[1][2] = A[1][2] - k*A[0][2];
        B[1] = B[1] - k*B[0];
        k = 0;
        Vivod(A, B);
    }
    if (A[2][0] == 0)
    {
    }
    else
    {
        k = A[2][0] / A[0][0];
        A[2][0] = A[2][0] - k*A[0][0];
        A[2][1] = A[2][1] - k*A[0][1];
        A[2][2] = A[2][2] - k*A[0][2];
        B[2] = B[2] - k*B[0];
        k = 0;
        Vivod(A, B);
    }
    if (A[2][1] == 0)
    {
    }
    else
    {
        if (A[1][1] / A[2][1] == A[1][2] / A[2][2])
        {
            cout << "Матрица вырожденна!" << endl;
            system("pause");
            return 0;
        }
        k = A[2][1] / A[1][1];
        A[2][0] = A[2][0] - k*A[1][0];
        A[2][1] = A[2][1] - k*A[1][1];
        A[2][2] = A[2][2] - k*A[1][2];
        B[2] = B[2] - k*B[1];
        Vivod(A, B);
    }
 
    X[2] = B[2] / A[2][2];
    X[1] = (B[1] - A[1][2] * X[2]) / A[1][1];
    X[0] = (B[0] - A[0][2] * X[2] - A[0][1] * X[1]) / A[0][0];
 
    cout << endl;
    cout << "Элементы массива X:" << endl;
    for (i = 0; i < n; i++){
        cout << "X" << i + 1 << "=" << X[i] << endl;
 
    }
 
    system("pause");
}