-
Notifications
You must be signed in to change notification settings - Fork 1
/
terminalka_3.c
78 lines (76 loc) · 2.14 KB
/
terminalka_3.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
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
int main(){
int n;
int k;
printf("Введите количество координат у вектора:");
scanf("%d",&n);
printf("Введите количество векторов:");
scanf("%d",&k);
float** A = (float**)malloc(k*sizeof(float*));
for(int i = 0; i<k; i++){
A[i] = (float*)malloc(n*sizeof(float));
}
printf("Введите координаты векторов:\n");
for(int i = 0; i<k; i++){
for(int j = 0; j<n; j++){
scanf("%f",&A[i][j]);
}
}
int i = 0;
int j = 0;
while(i<n){
while(j<k){
//max в столбце
float max = 0.0;
int r = 0;
for(int q = j; q<k; q++){
if(fabsf(A[q][i]) > max){
max = A[q][i];
r = q;
}
}
//перестановка
if(r){
for(int q = 0; q<n; q++){
float p = A[j][q];
A[j][q] = A[r][q];
A[r][q] = p;
}
}
//
if(A[j][i] == 0){
break;
}
//зануление
for(int q = j+1; q<k; q++){
float koef=-(A[q][i]/(float)A[j][i]);
for(int l = 0; l<n; l++){
A[q][l] = A[q][l] + koef*A[j][l];
}
}
j += 1;
i += 1;
}
i += 1;
}
int egnar = 0;
for(int i = 0; i<k; i++){
int count = 0;
for(int j = 0; j<n; j++){
if(fabsf(A[i][j]) == 0){
count += 1;
}
}
if(count == n){
egnar += 1;
}
}
printf("\nРанг системы векторов равен %d",(k-egnar));
if(egnar == 0){
printf(" и система линейно независима\n");
}else{
printf(" и система линейно зависима\n");
}
}