Помогите исправить программу решения системы уравнений

# include <stdio.h>
#include <iostream>
# include <math.h>
#include <conio.h>
#include <locale.h>
 
using namespace std;
 
void Iteration( double **A, double *b, double *x1, int n);
int F(int n, double *x1, double *x2);
void Full_1(double *A, int n);
void Full_2(double **A,int n);
 
int main(){
 
    // Значения для проверки:
 
    //const int n = 3;
    //double A[n][n] = { { 6.0, -1.0, -1.0 }, { -1.0, 6.0, -1.0 }, { -1.0, -1.0, 6.0 } };
    //double b[n] = { 11.0, 32.0, 42.0 };
    //double x1[n] = { 4.0, 4.0, 4.0 };
 
    int n;
    cout << "Enter n:" << endl;
    cin >> n;
 
    double **A = new double*[n];
    for (int i=0 ; i < n ; i++) 
        A[i] = new double[n];
 
    double *b = new double[n];
    double *x1 = new double[n];
 
    cout << "Enter A:" << endl;
    Full_2(A, n);
    cout << "Enter b:" << endl;
    Full_1(b, n);
    cout << "Enter x0:" << endl;
    Full_1(x1, n);
    
 
    Iteration( A, b, x1, n);
 
    _getch();
    return 0;
}
 
void Iteration( double **A, double *b, double *x1, int n){
 
    double *x2 = new double[n];
    for (int i=0 ; i < n ; i++) 
        x2[i] = 0;
 
    for (int i = 0; i < n; i++){
        double s = 0;
        for (int j = 0; j < n; j++)
            s = s + A[i][j] * x1[j];
        s = s - A[i][i] * x1[i];
        x2[i] = (1.0 / A[i][i]) * (b[i] - s);
    }
 
    int f = F(n, x1, x2);
 
    while ( f == 1){
        for (int i = 0; i < n; i++){
            double t = x2[i];
            double s = 0;
            for (int j = 0; j < n; j++)
                s = s + A[i][j] * x1[j];
            s = s - A[i][i] * x1[i];
            x2[i] = (1.0 / A[i][i]) * (b[i] - s);
            x1[i] = t;
        }
        f = F(n, x1, x2);
    }
 
    cout << "х:" << endl;
    for (int i = 0; i < n; i++)
        cout << x2[i] << endl;
 
}
 
// Условие выполнения итерации
 
int F(int n, double *x1, double *x2){
    double eps = 0.001;
    double max = fabs(x2[0] - x1[0]);
    for (int i = 1; i < n; i++){
        if (fabs(x2[i] - x1[i]) > max)
            max = fabs(x2[i] - x1[i]);
    }
    if (max > eps)
        return 1;
    else
        return 0;
 
}
 
//Инициализация массива:
 
void Full_1(double *A, int n){
    for (int i = 0; i < n; i++)
        scanf_s("%f", &A[i]);
}
 
void Full_2(double **A,int n){
    for (int i = 0 ; i < n ; i++)
        for (int j = 0 ; j < n ; j++)
            scanf_s("%f", &A[i][j]);
}

while( F(n, x1, x2) ){
...

 for (int i = 0; i < n; i++){
            double t = x2[i];
            double s = 0;
            for (int j = 0; j < n; j++)
                s = s + A[i][j] * x1[j];//используешь значения x1[i], но при вычислении x1[i+1] будет использоваться 
 
//значения x не от предыдущей итерации, а наполовину от предыдущей - наполовину от текущей!
            s = s - A[i][i] * x1[i];
            x2[i] = (1.0 / A[i][i]) * (b[i] - s);
            x1[i] = t;//Вот! тут же их меняешь, не закончив толком текущую итерацию
 
        }

    do{
        for (int i = 0; i < n; i++){
            double s = 0;
            for (int j = 0; j < n; j++)
                s = s + A[i][j] * x1[j];
            s = s - A[i][i] * x1[i];
            x2[i] = (1.0 / A[i][i]) * (b[i] - s);
        } 
        f=F(n, x1, x2);
        for (int i = 0; i < n; i++)
            x1[i]=x2[i];
    }while (f);
  cout << "х:" << endl;
    for (int i = 0; i < n; i++)
        cout << x1[i] << endl;

# include <stdio.h>
#include <iostream>
# include <math.h>
#include <conio.h>
#include <locale.h>
 
const int n = 3;
 
int F(int n, double x1[], double x2[]);
 
int main(){
    double A[n][n] = { { 6.0, -1.0, -1.0 }, { -1.0, 6.0, -1.0 }, { -1.0, -1.0, 6.0 } };
    double b[n] = { 11.0, 32.0, 42.0 };
    double x1[n] = { 4.0, 4.0, 4.0 };
    double x2[n] = { 0 };
 
    for (int i = 0; i < n; i++){
        double s = 0;
        for (int j = 0; j < n; j++)
            s = s + A[i][j] * x1[j];
        s = s - A[i][i] * x1[i];
        x2[i] = (1.0 / A[i][i]) * (b[i] - s);
    }
 
    int f = F(n, x1, x2);
 
 
    while ( f == 1){
        for (int i = 0; i < n; i++){
            double t = x2[i];
            double s = 0;
            for (int j = 0; j < n; j++)
                s = s + A[i][j] * x1[j];
            s = s - A[i][i] * x1[i];
            x2[i] = (1.0 / A[i][i]) * (b[i] - s);
            x1[i] = t;
        }
        f = F(n, x1, x2);
    }
 
    for (int i = 0; i < n; i++)
        std::cout << x2[i] << std::endl;
 
    _getch();
    return 0;
}
 
int F(int n, double x1[], double x2[]){
    double eps = 0.001;
    double max = fabs(x2[0] - x1[0]);
    for (int i = 1; i < n; i++){
        if (fabs(x2[i] - x1[i]) > max)
            max = fabs(x2[i] - x1[i]);
    }
    if (max > eps)
        return 1;
    else
        return 0;
 
}