teratail
回答率
85.30%
teratail header banner
teratail header banner
質問するログイン新規登録
トップに関する質問C言語 ガウスの消去法

編集履歴

質問編集履歴

7

題名変更

2020/09/02 03:08

投稿

pocpocpoc
pocpocpoc

スコア26

title CHANGED
@@ -1,1 +1,1 @@
1
- C言語 部分ピボット選択付ガウスの消去法
1
+ C言語 ガウスの消去法
body CHANGED
File without changes

6

間違っている回答やコードの部分を削除しました。

2020/09/02 03:08

投稿

pocpocpoc
pocpocpoc

スコア26

title CHANGED
File without changes
body CHANGED
@@ -9,7 +9,7 @@
9
9
  #include <stdlib.h>
10
10
  #include <math.h>
11
11
 
12
- void PrintAb(double *A, double *b, int n){/*行列A・ベクトルbを表示する関数*/
12
+ void PrintAb(double *A, double *b, int n){
13
13
  int i,j;
14
14
  double *ai;
15
15
  ai= A;
@@ -27,27 +27,24 @@
27
27
  int main(void) {
28
28
  int i, j, k, m, N;
29
29
  double r;
30
- double *ai, *ak; /* 各行の先頭へのポインタを表すポインタを用意する*/
30
+ double *ai, *ak;
31
- double *A, *b; /*計算結果を代入するための行列を用意 */
31
+ double *A, *b;
32
32
  double mmax, temp;
33
33
  printf("N= ");
34
34
  scanf("%d",&N);
35
35
 
36
- /*行列Aのメモリ確保*/
37
36
  A = (double *)malloc(N*N*sizeof(double));
38
37
  if(A==NULL){
39
38
  printf("Can’t allocate memory. Y=n");
40
39
  exit(1);
41
40
  }
42
41
 
43
- /*ベクトルbのメモリ確保*/
44
- b = (double *)malloc(N*sizeof(double));
42
+ b = (double *)malloc(N*sizeof(double));
45
43
  if(b==NULL){
46
44
  printf("Can’t allocate memory. Y=n");
47
45
  exit(1);
48
46
  }
49
47
 
50
- /*行列・ベクトルの入力*/
51
48
  ai = A;
52
49
  ak = A;
53
50
  for(i = 0 ; i < N ; i++ ){
@@ -66,31 +63,28 @@
66
63
 
67
64
  PrintAb(A, b, N);
68
65
 
69
- /* 前進消去（ピボット選択）*/
70
- for(k = 0; k < N-1; k++){ /* 第ｋステップ */
66
+ for(k = 0; k < N-1; k++){
71
67
  m = k;
72
68
  mmax = fabs(*(ak+k));
73
- for(i = k+1; i < N; i++){ /* ピボット選択 */
69
+ for(i = k+1; i < N; i++){
74
70
  if(fabs(*(ai+i*N+k)) > mmax){
75
71
  m = i;
76
72
  mmax = fabs(*(ai+i*N+k));
77
73
  }
78
74
  }
79
75
  printf("column %d: row %d is max\n", k, m);
80
- if(m != k){ /* 第ｋ行と第m行を交換　*/
76
+ if(m != k){
81
- for(i = k; i < N; i++){/* 係数行列　*/
77
+ for(i = k; i < N; i++){
82
78
  temp = *(ak+i);
83
79
  *(ak+i) = *(ak+m*N+i);
84
80
  *(ak+m*N+i) = temp;
85
81
  }
86
- /* 既知ベクトル */
87
82
  temp = b[k];
88
83
  b[k] = b[m];
89
84
  b[m] = temp;
90
85
  }
91
86
  ai = A;
92
87
  ak = A;
93
- /* 第k列消去 */
94
88
  for (i = k + 1; i < N; i++){
95
89
  r = *(ai+i*N+k) / *(ak+k);
96
90
  for (j = k; j < N; j++){
@@ -102,7 +96,6 @@
102
96
  ak+=N;
103
97
  }
104
98
 
105
- /*後退代入*/
106
99
  ai = ak;
107
100
  for (i = N - 1; i >= 0; i--) {
108
101
  for (j = i + 1; j < N; j++) {
@@ -164,47 +157,4 @@
164
157
  A, b =
165
158
  3.0 2.0 -1.0 | 1.0
166
159
  0.0 -2.3 1.7 | -1.0
167
- 0.0 -0.0 2.6 | 2.0
160
+ 0.0 -0.0 2.6 | 2.0
168
-
169
-
170
- =私の誤った出力=
171
- N= 3
172
- A[0][0] = 1
173
- A[0][1] = 1
174
- A[0][2] = 2
175
- A[1][0] = 3
176
- A[1][1] = 2
177
- A[1][2] = -1
178
- A[2][0] = 2
179
- A[2][1] = -1
180
- A[2][2] = 1
181
- b[0] = 4
182
- b[1] = -1
183
- b[2] = 5
184
-
185
- A, b =
186
- 1.0 1.0 2.0 | 4.0
187
- 3.0 2.0 -1.0 | -1.0
188
- 2.0 -1.0 1.0 | 5.0
189
- column 0: row 1 is max
190
- A, b =
191
- 3.0 2.0 -1.0 | -1.0
192
- 0.0 0.3 2.3 | 4.3
193
- 0.0 -2.3 1.7 | 5.7
194
- column 1: row 2 is max ←ここまでは正しいです。
195
- A, b =
196
- 3.0 2.0 -1.0 | 0.3
197
- 0.0 -1.0 4.3 | 5.7
198
- 0.0 0.0 0.5 | 8.9
199
- A, b =
200
- 3.0 2.0 -1.0 | 0.3
201
- 0.0 -1.0 4.3 | 5.7
202
- 0.0 0.0 0.5 | 2.1
203
- A, b =
204
- 3.0 2.0 -1.0 | 0.3
205
- 0.0 -1.0 4.3 | 3.9
206
- 0.0 0.0 0.5 | 2.1
207
- A, b =
208
- 3.0 2.0 -1.0 | 0.0
209
- 0.0 -1.0 4.3 | 3.9
210
- 0.0 0.0 0.5 | 2.1

5

ソースコードの更新

2020/09/02 03:07

投稿

pocpocpoc
pocpocpoc

スコア26

title CHANGED
File without changes
body CHANGED
@@ -16,9 +16,9 @@
16
16
  printf("A, b = \n");
17
17
  for( i = 0 ; i < n ; i++){
18
18
  for( j = 0 ; j < n ; j++){
19
- printf("%6.2f\t", *(ai+j));
19
+ printf("%6.1f\t", *(ai+j));
20
20
  }
21
- printf("|%6.2f\n",*(b+i));
21
+ printf("|%6.1f\n",*(b+i));
22
22
  ai+=n;
23
23
  }
24
24
  }

4

答えに一番近づいたソースコード

2020/09/01 15:22

投稿

pocpocpoc
pocpocpoc

スコア26

title CHANGED
File without changes
body CHANGED
@@ -2,6 +2,7 @@
2
2
 
3
3
  キーボードから次元 N を入力した後，N × N 行列 A と N 次元ベクトル b のデータを要素 ごとにキーボードから入力し，その行列・ベクトルを使って，連立 1 次方程式 Ax = b を 部分ピボット選択付きガウスの消去法を用いて解け.
4
4
  以上のような問いに対し、以下のソースコードを書きました。
5
+ 質問をしてから、何度もやり直しや修正を図り少しづつ答えには近づいてきたものの、まだ正確にはうまくいきません。よろしくお願いいたします。
5
6
 
6
7
  ```
7
8
  #include <stdio.h>
@@ -15,9 +16,9 @@
15
16
  printf("A, b = \n");
16
17
  for( i = 0 ; i < n ; i++){
17
18
  for( j = 0 ; j < n ; j++){
18
- printf("%6.1f\t", *(ai+j));
19
+ printf("%6.2f\t", *(ai+j));
19
20
  }
20
- printf("|%6.1f\n",*(b+i));
21
+ printf("|%6.2f\n",*(b+i));
21
22
  ai+=n;
22
23
  }
23
24
  }
@@ -68,9 +69,9 @@
68
69
  /* 前進消去（ピボット選択）*/
69
70
  for(k = 0; k < N-1; k++){ /* 第ｋステップ */
70
71
  m = k;
71
- mmax = fabs(*(ai+k*N+k));
72
+ mmax = fabs(*(ak+k));
72
73
  for(i = k+1; i < N; i++){ /* ピボット選択 */
73
- if(fabs(*A) > mmax){
74
+ if(fabs(*(ai+i*N+k)) > mmax){
74
75
  m = i;
75
76
  mmax = fabs(*(ai+i*N+k));
76
77
  }
@@ -78,7 +79,7 @@
78
79
  printf("column %d: row %d is max\n", k, m);
79
80
  if(m != k){ /* 第ｋ行と第m行を交換　*/
80
81
  for(i = k; i < N; i++){/* 係数行列　*/
81
- temp= *(ak+i);
82
+ temp = *(ak+i);
82
83
  *(ak+i) = *(ak+m*N+i);
83
84
  *(ak+m*N+i) = temp;
84
85
  }
@@ -90,17 +91,15 @@
90
91
  ai = A;
91
92
  ak = A;
92
93
  /* 第k列消去 */
93
- for (k = 0; k < N - 1; k++){
94
- for (i = k + 1; i < N; i++){
94
+ for (i = k + 1; i < N; i++){
95
- r = *(ai+i*N+k) / *(ak+k);
95
+ r = *(ai+i*N+k) / *(ak+k);
96
- for (j = k; j < N; j++){
96
+ for (j = k; j < N; j++){
97
- *(ai+i*N+j) -= *(ak+j) * r;
97
+ *(ai+i*N+j) -= *(ak+j) * r;
98
- }
99
- b[i] -= b[k] * r;
100
98
  }
101
- ak+=N;
99
+ b[i] -= b[k] * r;
102
100
  }
103
101
  PrintAb(A, b, N);
102
+ ak+=N;
104
103
  }
105
104
 
106
105
  /*後退代入*/
@@ -122,7 +121,8 @@
122
121
  }
123
122
  ```
124
123
 
125
- しかし、私自身が書いたソースコードを実行してもうまくいきません(エラーなどは出ていません)。
124
+ しかし、私自身が書いたソースコードを実行しても正しくいきません(エラーなどは出ていません)。
125
+ (第二回目の消去までは正しいのですが、それ以降が正しくありません。)
126
126
  どなたかわかる方がいらっしゃれば教えていただけると幸いです。
127
127
  =正しい出力=
128
128
  N= 3
@@ -140,31 +140,31 @@
140
140
  b[2] = 5
141
141
 
142
142
  A, b =
143
- 1.00 1.00 2.00 | 4.00
143
+ 1.0 1.0 2.0 | 4.0
144
- 3.00 2.00 -1.00 | -1.00
144
+ 3.0 2.0 -1.0 | -1.0
145
- 2.00 -1.00 1.00 | 5.00
145
+ 2.0 -1.0 1.0 | 5.0
146
146
  column 0: row 1 is max
147
147
  A, b =
148
- 3.00 2.00 -1.00 | -1.00
148
+ 3.0 2.0 -1.0 | -1.0
149
- 0.00 0.30 2.30 | 4.30
149
+ 0.0 0.3 2.3 | 4.3
150
- 0.00 -2.30 1.70 | 5.70
150
+ 0.0 -2.3 1.7 | 5.7
151
151
  column 1: row 2 is max
152
152
  A, b =
153
- 3.00 2.00 -1.00 | -1.00
153
+ 3.0 2.0 -1.0 | -1.0
154
- 0.00 -2.30 1.70 | 5.70
154
+ 0.0 -2.3 1.7 | 5.7
155
- 0.00 -0.00 2.60 | 5.10
155
+ 0.0 -0.0 2.6 | 5.1
156
156
  A, b =
157
- 3.00 2.00 -1.00 | -1.00
157
+ 3.0 2.0 -1.0 | -1.0
158
- 0.00 -2.30 1.70 | 5.70
158
+ 0.0 -2.3 1.7 | 5.7
159
- 0.00 -0.00 2.60 | 2.00
159
+ 0.0 -0.0 2.6 | 2.0
160
160
  A, b =
161
- 3.00 2.00 -1.00 | -1.00
161
+ 3.0 2.0 -1.0 | -1.0
162
- 0.00 -2.30 1.70 | -1.00
162
+ 0.0 -2.3 1.7 | -1.0
163
- 0.00 -0.00 2.60 | 2.00
163
+ 0.0 -0.0 2.6 | 2.0
164
164
  A, b =
165
- 3.00 2.00 -1.00 | 1.00
165
+ 3.0 2.0 -1.0 | 1.0
166
- 0.00 -2.30 1.70 | -1.00
166
+ 0.0 -2.3 1.7 | -1.0
167
- 0.00 -0.00 2.60 | 2.00
167
+ 0.0 -0.0 2.6 | 2.0
168
168
 
169
169
 
170
170
  =私の誤った出力=
@@ -183,23 +183,28 @@
183
183
  b[2] = 5
184
184
 
185
185
  A, b =
186
- 1.00 1.00 2.00 | 4.00
186
+ 1.0 1.0 2.0 | 4.0
187
- 3.00 2.00 -1.00 | -1.00
187
+ 3.0 2.0 -1.0 | -1.0
188
- 2.00 -1.00 1.00 | 5.00
188
+ 2.0 -1.0 1.0 | 5.0
189
189
  column 0: row 1 is max
190
190
  A, b =
191
- 3.00 2.00 -1.00 | -1.00
191
+ 3.0 2.0 -1.0 | -1.0
192
- 0.00 0.33 2.33 | 4.33
192
+ 0.0 0.3 2.3 | 4.3
193
- 0.00 0.00 18.00 | 36.00
193
+ 0.0 -2.3 1.7 | 5.7
194
+ column 1: row 2 is max ←ここまでは正しいです。
194
195
  A, b =
195
- 3.00 2.00 -1.00 | -1.00
196
+ 3.0 2.0 -1.0 | 0.3
196
- 0.00 0.33 2.33 | 4.33
197
+ 0.0 -1.0 4.3 | 5.7
197
- 0.00 0.00 18.00 | 2.00
198
+ 0.0 0.0 0.5 | 8.9
198
199
  A, b =
199
- 3.00 2.00 -1.00 | -1.00
200
+ 3.0 2.0 -1.0 | 0.3
200
- 0.00 0.33 2.33 | -1.00
201
+ 0.0 -1.0 4.3 | 5.7
201
- 0.00 0.00 18.00 | 2.00
202
+ 0.0 0.0 0.5 | 2.1
202
203
  A, b =
204
+ 3.0 2.0 -1.0 | 0.3
205
+ 0.0 -1.0 4.3 | 3.9
206
+ 0.0 0.0 0.5 | 2.1
207
+ A, b =
203
- 3.00 2.00 -1.00 | 1.00
208
+ 3.0 2.0 -1.0 | 0.0
204
- 0.00 0.33 2.33 | -1.00
209
+ 0.0 -1.0 4.3 | 3.9
205
- 0.00 0.00 18.00 | 2.00
210
+ 0.0 0.0 0.5 | 2.1

3

正しい答えに最も近づいたソースコードへの変更

2020/09/01 14:02

投稿

pocpocpoc
pocpocpoc

スコア26

title CHANGED
File without changes
body CHANGED
@@ -15,9 +15,9 @@
15
15
  printf("A, b = \n");
16
16
  for( i = 0 ; i < n ; i++){
17
17
  for( j = 0 ; j < n ; j++){
18
- printf("%6.2f\t", *(ai+j));
18
+ printf("%6.1f\t", *(ai+j));
19
19
  }
20
- printf("|%6.2f\n",*(b+i));
20
+ printf("|%6.1f\n",*(b+i));
21
21
  ai+=n;
22
22
  }
23
23
  }
@@ -31,21 +31,21 @@
31
31
  double mmax, temp;
32
32
  printf("N= ");
33
33
  scanf("%d",&N);
34
-
34
+
35
35
  /*行列Aのメモリ確保*/
36
36
  A = (double *)malloc(N*N*sizeof(double));
37
37
  if(A==NULL){
38
38
  printf("Can’t allocate memory. Y=n");
39
39
  exit(1);
40
40
  }
41
-
41
+
42
42
  /*ベクトルbのメモリ確保*/
43
43
  b = (double *)malloc(N*sizeof(double));
44
44
  if(b==NULL){
45
45
  printf("Can’t allocate memory. Y=n");
46
46
  exit(1);
47
47
  }
48
-
48
+
49
49
  /*行列・ベクトルの入力*/
50
50
  ai = A;
51
51
  ak = A;
@@ -56,13 +56,15 @@
56
56
  }
57
57
  ai+=N;
58
58
  }
59
-
59
+
60
60
  for(i = 0 ; i < N ; i++ ){
61
61
  printf("b[%d] = ", i );
62
62
  scanf("%lf", b+i);
63
63
  }
64
64
  printf("\n");
65
65
 
66
+ PrintAb(A, b, N);
67
+
66
68
  /* 前進消去（ピボット選択）*/
67
69
  for(k = 0; k < N-1; k++){ /* 第ｋステップ */
68
70
  m = k;
@@ -74,7 +76,7 @@
74
76
  }
75
77
  }
76
78
  printf("column %d: row %d is max\n", k, m);
77
- if(m != k){ /* 第ｋ行と第m行の交換　*/
79
+ if(m != k){ /* 第ｋ行と第m行を交換　*/
78
80
  for(i = k; i < N; i++){/* 係数行列　*/
79
81
  temp= *(ak+i);
80
82
  *(ak+i) = *(ak+m*N+i);
@@ -85,26 +87,22 @@
85
87
  b[k] = b[m];
86
88
  b[m] = temp;
87
89
  }
90
+ ai = A;
91
+ ak = A;
88
- PrintAb(A, b, N);
92
+ /* 第k列消去 */
93
+ for (k = 0; k < N - 1; k++){
94
+ for (i = k + 1; i < N; i++){
95
+ r = *(ai+i*N+k) / *(ak+k);
96
+ for (j = k; j < N; j++){
97
+ *(ai+i*N+j) -= *(ak+j) * r;
89
- }
98
+ }
90
-
91
- PrintAb(A, b, N);
92
-
93
- /*前進消去*/
94
- ai = A;
95
- ak = A;
96
- for (k = 0; k < N - 1; k++){
97
- for (i = k + 1; i < N; i++){
98
- r = *(ai+i*N+k) / *(ak+k);
99
- for (j = k; j < N; j++){
100
- *(ai+i*N+j) -= *(ak+j) * r;
99
+ b[i] -= b[k] * r;
101
100
  }
102
- b[i] -= b[k] * r;
101
+ ak+=N;
103
102
  }
104
103
  PrintAb(A, b, N);
105
- ak+=N;
106
104
  }
107
-
105
+
108
106
  /*後退代入*/
109
107
  ai = ak;
110
108
  for (i = N - 1; i >= 0; i--) {
@@ -115,11 +113,11 @@
115
113
  PrintAb(A, b,N);
116
114
  ai-=N;
117
115
  }
118
-
116
+
119
117
  /*メモリの解放*/
120
118
  free(A);
121
119
  free(b);
122
-
120
+
123
121
  return 0;
124
122
  }
125
123
  ```
@@ -184,37 +182,24 @@
184
182
  b[1] = -1
185
183
  b[2] = 5
186
184
 
187
- column 0: row 1 is max
188
185
  A, b =
186
+ 1.00 1.00 2.00 | 4.00
189
187
  3.00 2.00 -1.00 | -1.00
190
- 1.00 1.00 2.00 | 4.00
191
188
  2.00 -1.00 1.00 | 5.00
192
- column 1: row 2 is max
189
+ column 0: row 1 is max
193
190
  A, b =
194
- 3.00 -1.00 1.00 | -1.00
191
+ 3.00 2.00 -1.00 | -1.00
192
+ 0.00 0.33 2.33 | 4.33
195
- 1.00 1.00 2.00 | 5.00
193
+ 0.00 0.00 18.00 | 36.00
196
- 2.00 2.00 -1.00 | 4.00
197
194
  A, b =
198
- 3.00 -1.00 1.00 | -1.00
195
+ 3.00 2.00 -1.00 | -1.00
196
+ 0.00 0.33 2.33 | 4.33
199
- 1.00 1.00 2.00 | 5.00
197
+ 0.00 0.00 18.00 | 2.00
200
- 2.00 2.00 -1.00 | 4.00
201
198
  A, b =
202
- 3.00 -1.00 1.00 | -1.00
199
+ 3.00 2.00 -1.00 | -1.00
203
- 0.00 1.33 1.67 | 5.33
200
+ 0.00 0.33 2.33 | -1.00
204
- 0.00 2.67 -1.67 | 4.67
201
+ 0.00 0.00 18.00 | 2.00
205
202
  A, b =
206
- 3.00 -1.00 1.00 | -1.00
203
+ 3.00 2.00 -1.00 | 1.00
207
- 0.00 1.33 1.67 | 5.33
204
+ 0.00 0.33 2.33 | -1.00
208
- 0.00 0.00 -5.00 | -6.00
205
+ 0.00 0.00 18.00 | 2.00
209
- A, b =
210
- 3.00 -1.00 1.00 | -1.00
211
- 0.00 1.33 1.67 | 5.33
212
- 0.00 0.00 -5.00 | 1.20
213
- A, b =
214
- 3.00 -1.00 1.00 | -1.00
215
- 0.00 1.33 1.67 | 2.50
216
- 0.00 0.00 -5.00 | 1.20
217
- A, b =
218
- 3.00 -1.00 1.00 | 0.10
219
- 0.00 1.33 1.67 | 2.50
220
- 0.00 0.00 -5.00 | 1.20

2

ソースコードの修正

2020/09/01 13:40

投稿

pocpocpoc
pocpocpoc

スコア26

title CHANGED
@@ -1,1 +1,1 @@
1
- C言語 ピボット選択付ガウスの消去法
1
+ C言語 部分ピボット選択付ガウスの消去法
body CHANGED
@@ -81,9 +81,9 @@
81
81
  *(ak+m*N+i) = temp;
82
82
  }
83
83
  /* 既知ベクトル */
84
- m = b[k];
84
+ temp = b[k];
85
85
  b[k] = b[m];
86
- b[m] = m;
86
+ b[m] = temp;
87
87
  }
88
88
  PrintAb(A, b, N);
89
89
  }
@@ -186,35 +186,35 @@
186
186
 
187
187
  column 0: row 1 is max
188
188
  A, b =
189
- 3.00 2.00 -1.00 | 0.00
189
+ 3.00 2.00 -1.00 | -1.00
190
- 1.00 1.00 2.00 | -1.00
190
+ 1.00 1.00 2.00 | 4.00
191
191
  2.00 -1.00 1.00 | 5.00
192
192
  column 1: row 2 is max
193
193
  A, b =
194
- 3.00 -1.00 1.00 | 0.00
194
+ 3.00 -1.00 1.00 | -1.00
195
- 1.00 1.00 2.00 | 0.00
195
+ 1.00 1.00 2.00 | 5.00
196
- 2.00 2.00 -1.00 | 5.00
196
+ 2.00 2.00 -1.00 | 4.00
197
197
  A, b =
198
- 3.00 -1.00 1.00 | 0.00
198
+ 3.00 -1.00 1.00 | -1.00
199
- 1.00 1.00 2.00 | 0.00
199
+ 1.00 1.00 2.00 | 5.00
200
- 2.00 2.00 -1.00 | 5.00
200
+ 2.00 2.00 -1.00 | 4.00
201
201
  A, b =
202
- 3.00 -1.00 1.00 | 0.00
202
+ 3.00 -1.00 1.00 | -1.00
203
- 0.00 1.33 1.67 | 0.00
203
+ 0.00 1.33 1.67 | 5.33
204
- 0.00 2.67 -1.67 | 5.00
204
+ 0.00 2.67 -1.67 | 4.67
205
205
  A, b =
206
- 3.00 -1.00 1.00 | 0.00
206
+ 3.00 -1.00 1.00 | -1.00
207
- 0.00 1.33 1.67 | 0.00
207
+ 0.00 1.33 1.67 | 5.33
208
- 0.00 0.00 -5.00 | 5.00
208
+ 0.00 0.00 -5.00 | -6.00
209
209
  A, b =
210
- 3.00 -1.00 1.00 | 0.00
210
+ 3.00 -1.00 1.00 | -1.00
211
- 0.00 1.33 1.67 | 0.00
211
+ 0.00 1.33 1.67 | 5.33
212
- 0.00 0.00 -5.00 | -1.00
212
+ 0.00 0.00 -5.00 | 1.20
213
213
  A, b =
214
- 3.00 -1.00 1.00 | 0.00
214
+ 3.00 -1.00 1.00 | -1.00
215
- 0.00 1.33 1.67 | 1.25
215
+ 0.00 1.33 1.67 | 2.50
216
- 0.00 0.00 -5.00 | -1.00
216
+ 0.00 0.00 -5.00 | 1.20
217
217
  A, b =
218
- 3.00 -1.00 1.00 | 0.75
218
+ 3.00 -1.00 1.00 | 0.10
219
- 0.00 1.33 1.67 | 1.25
219
+ 0.00 1.33 1.67 | 2.50
220
- 0.00 0.00 -5.00 | -1.00
220
+ 0.00 0.00 -5.00 | 1.20

1

2020/09/01 13:30

投稿

pocpocpoc
pocpocpoc

スコア26

title CHANGED
@@ -1,1 +1,1 @@
1
- C言語 ピボットピボット選択付ガウスの消去法
1
+ C言語 ピボット選択付ガウスの消去法
body CHANGED
File without changes