Приветствую, подскажите в чем ошибка.
| C++ |
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
| int gauss(int n, double **A, double **&E)
{
//создание единичной матрицы
for (int i = 0; i < n; i++)
for (int j = 0; j < n; j++)
{
if (i == j)
E[i][j] = 1;
else
E[i][j] = 0;
}
//приведение к верхнетреугольному виду и нахождение определителя
for (int step = 0; step<n - 1; step++)
for (int row = step + 1; row < n; row++)
{
double coeff = -A[row][step] / A[step][step];
for (int col = step; col < n; col++)
{
A[row][col] += A[step][col] * coeff;
E[row][col] += E[step][col] * coeff;
}
}
double e;
for (int i = 0; i < n; i++)
det *= A[i][i];
cout << "Определитель А = " << det << endl << endl;
//получаем единицы в главной диагонали
if (det != 0)
{
for (int k = 0; k < n; k++)
{
double t = A[k][k];
for (int j = 0; j < n; j++)
{
A[k][j] /= t;
E[k][j] /= t;
}
}
for (int step = n-1; step > 0; step--)
for (int row = step - 1; row >=0; row--)
{
double coeff = -A[row][step];
for (int col = n - 1; col >= 0; col--)
{
A[row][col] += A[step][col] * coeff;
E[row][col] += E[step][col] * coeff;
}
}
cout << "Обратная матрица А" << endl;
print_matrix(n, n, E);
}
return 1;
} |
|
Приведение к верхнетреугольному виду и нахождение определителя работает нормально, проблемы начинаются при получении нулей ниже глав. диагонали. Никак не могу понять, в чем причина некорректных значений.
Добавлено через 1 час 0 минут
Пробовал по разному делать - та же проблема - не все значения корректны.
(
