Метод пополнения. Mathematica for student

[Onni]

Если вы разбираетесь в Mathematica for student,тогда ваш совет просто необходим...
имеется следущее в Mathematica for student:

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)

Кто знает где ошибка?

Оформляй код! // Модератор

[voproshaika2]

Ошибка в мозгу

[Z1000000]

Приложи промежуточные результаты для Математики, также как и для C++

[Onni]

Сообщение от Z1000000

Приложи промежуточные результаты для Математики, также как и для C++

ээээ.... не поняла немного....

Модератор, спасибо)))