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
89
90
91
92
93
| #include <stdio.h>
#include <math.h>
#include <Windows.h>
void Showmatrix(int n,double a[n][n],double *b)
{
int i,j;
for ( i = 1; i <= n; i++ )
{
for ( j = 1; j <= n ; j++ )
printf( "%.2f ",a[i][j]);
printf(" ");
printf( "%.2f",b[i]);
printf("\n");
}
}
void PryamoiHod(int n, double a[n][n], double *b)
{
int k,j,i;
double d;
for (k = 1; k <= n; k++) // прямой ход
{
if ( fabs( a[k] [k] ) == 0 )
{
printf( "Система не имеет единственного решения \n" );
}
for (j = k + 1; j <= n; j++)
{
d = a[j][k] / a[k][k]; // формула (1)
for (i = 1; i <= n; i++)
{
a[j][i] = a[j][i] - d * a[k][i]; // формула (2)
}
b[j] = b[j] - d * b[k]; // формула (3)
}
}
}
void ObratniHod(int n, double a[n][n], double *b,double *x)
{
int k,j;
double d,s;
for (k = n; k >= 1; k--) // обратный ход
{
d = 0;
for (j = k + 1; j <= n; j++)
{
s = a[k][j] * x[j]; // формула (4)
d = d + s; // формула (4)
}
x[k] = (b[k] - d) / a[k][k]; // формула (4)
}
}
int main()
{
SetConsoleCP(1251);
SetConsoleOutputCP(1251);
int n, i, j;
printf( "Введите число уравнений системы : \n" );
scanf( "%d", & n );
printf("\n");
double a[n][n];
double b[n];
double x[n];
printf("Введите матрицу коэфициентов (a[]) и свободных членов (b[])\n");
for ( i = 1; i <= n; i++ )
{
printf("\n");
for ( j = 1; j <= n ; j++ )
{
printf( "a[%d][%d] = ",i,j);
a[i][j]=-10+rand()%10;
//scanf( "%lf", & a[i][j] );
}
printf( "b[%d] = ",i);
scanf("%lf",&b[i]);
}
printf("\n");
printf("Получившаяся матрица : \n\n");
Showmatrix(n,a,b);
printf("\n");
PryamoiHod(n,a,b);
printf("Полученная треугольная матрица: \n\n");
Showmatrix(n,a,b);
ObratniHod(n,a,b,x);
printf("Корни системы: \n");
for( i = 1; i <= n; i++)
printf("x[%d] = %.2f\n",i,x[i]);
return 0;
} |
(
