Решения матрицы методом Гаусса

#include "stdafx.h"
#include "iostream"
#include "math.h"
#include "stdio.h"
 
using namespace std;
 
int main(void)
 
{
    int n = 4;
    double A[4][4] = { 
    { 8.3, 2.62 + 0.02, 4.1, 1.9 }, 
    { 3.92, 8.45, 7.78 - 0.02, 2.45 }, 
    { 3.77, 7.21 + 0.02, 8.04, 2.28 },
    { 2.21, 3.65 - 0.02, 1.69, 6.69 } 
    };
    double B[4] = { -10.65 + 0.02, 12.21, 15.45 - 0.02, -8.35 };
    double V[4][4], C[4][4], P[4], X[4], Y[4];
    for (int i = 0; i < n; i++)
    for (int j = 0; j < n; j++)
    {
        V[i][j] = A[i][j]; //это для того чтоб потом вывести самую первую матрицу//
        P[i] = B[i];
    }
 
 
    //прямой ход//
    for (int k = 0; k < n; k++)
    {
        Y[k] = P[k] / V[k][k];
        for (int i = k + 1; i < n; i++)
        {
            P[i] += -V[i][k] * Y[k];
            for (int j = k + 1; j < n; j++)
            {
                C[k][j] = V[k][j] / V[k][k];
                V[i][j] += -V[i][k] * C[k][j];
            }
        }
        }
        
    //обратный ход//
for (int i = 0; i < n; i++)
    X[i] = 0;
   X[n - 1] = Y[n - 1];
   for (int i = n - 2; i >= 0; i--)
   {
       for (int j = i + 1; j < n; j++)
           X[i] += C[i][j] * X[j];
       X[i] = Y[i] - X[i];
   }
 
   //вывод//
   for (int i = 0; i < n; i++)
   {
       for (int j = 0; j < n; j++)
           cout << V[i][j] << " ";
       cout << endl;
   }
    system("pause");
}