помогите пожалуйста решить СЛАУ LU - разложением.
| 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
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
| #include <iostream>
using namespace std;
int main ()
{setlocale(LC_ALL,"rus");
int n;
double sum = 0;
cout << "Введите порядок матрицы\n n = ";
cin >> n;
double A [ n ][ n ];
double L [ n ][ n ];
double U [ n ][ n ];
//задаем матрицу A[][] ...
for (int i = 0; i < n; i++)
{
for (int j = 0; j < n; j++)
{
cout << "\na[" << i << "][" << j << "] = ";
cin >> A [i][j];
L [i][j] = 0;
U [i][j] = 0;
if (i == j)
U [i][j] = 1;
}
}
//==============================================
//находим первый столбец L[][] и первую строку U[][]
for (int i = 0; i < n; i++)
{
L [i][0] = A [i][0];
U [0][i] = A [0][i] / L [0][0];
}
//дальше вычисляем L[][], U[][] по формуле
for (int i = 1; i < n; i++)
{
for (int j = 1; j < n; j++)
{
if (i >= j) //нижний треугольник
{
sum = 0;
for (int k = 0; k < j; k++)
sum += L [i][k] * U [k][j];
L [i][j] = A [i][j] - sum;
}
else // верхний
{
sum = 0;
for (int k = 0; k < i; k++)
sum += L [i][k] * U [k][j];
U [i][j] = (A [i][j] - sum) / L [i][i];
}
}
}
//====================================================
cout << "\n";
cout << "Полученные матрицы: " << endl;
cout << "\n";
cout << "Матрица L : " << endl << endl;
for (int i = 0; i < n; i++)
{
for (int j = 0; j < n; j++)
cout << " " << L [i][j] << " ";
cout << "\n\n";
}
cout << "\n\n";
cout << "Матрица U : " << endl << endl;
for (int i = 0; i < n; i++)
{
for (int j = 0; j < n; j++)
cout << " " << U [i][j] << " ";
cout << "\n\n";
}
return 0;
} |
|
Заранее спасибо.
(
