Метод пополнения. Метод Ершова - C++ - Ответ 5448900

Clear[a, A, dop];
A = a = ( {
     {1, 0.42, 0.54, 0.66},
     {0.42, 1, 0.32, 0.44},
     {0.54, 0.32, 1, 0.22},
     {0.66, 0.44, 0.22, 1}
    } );
dop = Table[0, {4}, {4}];
For[k = 1, k <= 4,
  For[i = 1, i <= Length[a],
   For[j = 1, j <= Length[a],
    If[i == k && j == k, dop[[i, j]] = 1/A[[k, k]]];
    If[i != k && j == k, dop[[i, j]] = A[[i, k]]/A[[k, k]]];
    If[i == k && j != k, dop[[i, j]] = -( A[[k, j]]/A[[k, k]])];
    If[i != k && j != k, 
     dop[[i, j]] = A[[i, j]] - A[[i, k]]/A[[k, k]]*A[[k, j]]];
    Print["k=", k, " ;i=", i, " ; j=", j, " ; a=", MatrixForm[dop]]
     j++];
   i++];
  A = dop;
  k++];
Print["Inverse Matrix=", MatrixForm[A]];

#include "stdafx.h"
#include <iostream>
using namespace std;
int main()
{
    const int n = 4;
//Clear[a, A, dop];
 
    double a[n][n] = {1.0, 0.42, 0.54, 0.66,
                      0.42, 1.0, 0.32, 0.44, 
                      0.54, 0.32, 1.0, 0.22,
                      0.66, 0.44, 0.22, 1.0}, 
    dop[n][n];
    for(int k = 0; k < n;k++){
        for(int i = 0; i <n;i++){
            for(int j = 0; j <n;j++){
    if((i == k) && (j == k))
        dop[i][j] = 1/(a[k][k]);
    if((i != k) && (j == k))
        dop[i][j] = (a[i][k])/(a[k][k]);
    if((i == k) && (j != k))
        dop[i][j] = -((a[k][j])/(a[k][k]));
    if((i != k) && (j != k)) 
     dop[i][j] = a[i][j] - (((a[i][k])/(a[k][k]))*(a[k][j]));
    cout<<"k= "<< k<< " i="<< i<< "  j= "<< j<< "  a= "<< dop[i][j]<<endl;
        a[i][j] = dop[i][j];    }}
        //a[i][j] = dop[i][j];
    }
    cout<<"Inverse Matrix = "<<endl;
    for(int i = 0; i < n; i++)
    {
    for(int j = 0; j < n; j++)
        cout<<a[i][j]<<"  ";
    }
    cin.get();
    cin.get();
 
    return 0;
}

My Inverse Matrix=(2.50759    -0.123039 -1.01149    -1.37834
                         -0.123039   1.33221    -0.261427   -0.447454
                          -1.01149   -0.261427  1.53183 0.445609
                          -1.37834  -0.447454 0.445609   2.00855)