Ошибка в решении СЛАУ методом Гаусса

#include<math.h>
#include<conio.h>
#include<stdio.h>
#include<iostream.h>
 
main()
{
int n,v,p,c,s,d,h;
 double **a,
        **a1,
        *x,
        *b,
         l,
         sum;
n=4; c=n-1;h=0;
    b = new double[n];
    x = new double[n];
    a = new double*[n];
    a1 = new double*[n];
    for(int i=0; i < n; i++){
        a[i] = new double[n];
        a1[i]= new double[n];
    }
cout<<"Koeff."<<endl;
    for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
cin>>*(*(a+i)+j);
cout<<endl<<"Svobod._chleni"<<endl;
for(int i=0;i<n;i++)
cin>>*(b+i);
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
a1[i][j]=a[i][j];
v=0;p=0;s=0;d=0;                    //прямой ход, начало
 
 for(int j=0;j<n-1;j++)
           { l=a[s][d];
 
 
        for (int i=0; i<n;i++)
            {a[h][i]=a[h][i]/l;}b[h]=b[h]/l;
   for(int g=0;g<c;g++)
   { for(int i=0;i<n;i++)
    {a[v+1][i]=a[v+1][i]-(a[v+1][i]*a[p][i]);}
    b[v+1]=b[v+1]-(b[v+1]*b[p]); v++;}p++;c--;s++;d++;h++;v=s;}
 
 
     for(int i=0;i<n;i++){cout<<endl;
for(int j=0;j<n;j++)
cout<<*(*(a+i)+j)<<" ";}
cout<<endl<<endl;
for(int i=0;i<n;i++)
cout<<" "<<b[i]<<" ";
cout<<endl<<endl;                               //прямой ход, конец
 
for(int m=n-1; m >= 0; m--){                //обратный ход
        double sum=0;
        for(int i=n-1; i > m; i--){
            sum += x[i] * a[m][i];
        }
        x[m] = (b[m] - sum)/a[m][m];
    }                                                 //обратный ход
 
    cout << "Reshenie" << endl;
    for(int i=0; i < n; i++){
        cout << "x[" << i+1 << "]=" << x[i] << endl;
    }
   
cout<<endl<<"Proverka"<<endl;
    for(int i=0;i<n;i++)                  //
    {sum=0;                               //не сходятся ответы
    for(int j=0;j<n;j++)                //
    {sum+=a1[i][j]*x[j];}             //
    cout<<i+1<<"="<<sum<<endl;}// 
getch();
}

#include <iostream>
#include <cstdlib>
#include <ctime>
 
// printing a matrix;
void print_extended_matrix(double ** matrix, double * b, unsigned int size)
{
    for (unsigned int i = 0u; i < size; ++i)
    {
    for (unsigned int j = 0u; j < size; ++j)
        std::cout << matrix[i][j] << " ";
 
    std::cout << " = " << b[i] << std::endl;
    }
}
 
int main()
{
    // Entering the data;
    unsigned int size;
 
    std::cout << "Enter the size of the matrix" << std::endl;
    std::cin >> size;
 
    double ** A = new double*[size];
    double * b = new double[size];
 
    std::cout << "Enter the matrix A" << std::endl;
    for (unsigned int i = 0u; i < size; ++i)
    {
    A[i] = new double[size];
    for (unsigned int j = 0u; j < size; ++j)
        std::cin >> A[i][j];
    }
 
    std::cout << "Enter the rigth vector B" << std::endl;
    for (unsigned int i = 0u; i < size; ++i)
    std::cin >> b[i];
 
    print_extended_matrix(A, b, size);
 
    // Gausse method - calculating the vector X, where A * X = B;
    for (unsigned int i = 0u; i < size; ++i)
    {
    // norming the current row, which will be subtracted from others;
    double value = A[i][i];
    for (unsigned int j = 0u; j < size; ++j)
        A[i][j] /= value;
    b[i] /= value;
 
    // the Gausse iteration;
    for (unsigned int j = 0u; j < i; ++j)
    {
        value = A[j][i];
        for (unsigned int k = i; k < size; ++k)
        A[j][k] -= A[i][k] * value;
 
        b[j] -= b[i] * value;
    }
 
    for (unsigned int j = i + 1; j < size; ++j)
    {
        value = A[j][i];
        for (unsigned int k = i; k < size; ++k)
        A[j][k] -= A[i][k] * value;
 
        b[j] -= b[i] * value;
    }
    }
 
    std::cout << "Vector X: " << std::endl;
    for (unsigned int i = 0u; i < size; ++i)
    {
    std::cout << b[i] << " ";
    delete[]A[i];
    }
 
    delete[]A;
    delete[]b;
 
    return EXIT_SUCCESS;
}