5ちゃんねる ★スマホ版★ ■掲示板に戻る■ 全部 1- 最新50  

■ このスレッドは過去ログ倉庫に格納されています

Gauss-Jordanの掃出法

1 :名無しさん@1周年:2001/06/25(月) 17:01
で、下のCプログラムを書き換えてpivot=0でもちゃんと計算できる
ようにだったかな?とにかくそんなふうに下のプログラムを完全なものに
したいのですが・・・
どうやったらいいのですか?
ここをこうかえる、というように具体的に教えて下さい。
マジでお願いします。

#include <stdio.h>
#define IMAX 3
#define JMAX 4


double array[IMAX][JMAX]={

{ 2, 1,-1, 5 },

{-3, 3, 2, 1 },

{ 1,-2,-2,-1 },

};


void PrintMat(void){
int i,j;
for(i=0;i<IMAX;i++){
for(j=0;j<JMAX;j++){
printf("%5.2f ",array[i][j]);
}
putchar('\n');
}
putchar('\n');
}


void SweepOut(void){
int i,j,axis;
double pivot,aik;
for(axis=0;axis<IMAX;axis++){
pivot=array[axis][axis];
for(j=axis;j<JMAX;j++){
array[axis][j]/=pivot;
}
for(i=0;i<IMAX;i++){
if(i!=axis){
aik=array[i][axis];
for(j=0;j<JMAX;j++){
array[i][j]-=array[axis][j]*aik;
}
}
}
}
}


main(){
PrintMat();
SweepOut();
PrintMat();
}


2 :1:2001/06/25(月) 17:06
pivotが0になるのを防ぐように改良するんだったかも・・・

3 :仕様書無しさん:2001/06/25(月) 21:43
『C言語によるアルゴリズム事典』に載ってなかったっけ?

4 ::2001/06/25(月) 22:07
pivot=0や、pivot=とても小さい数になるのを防ぐように
改良するんだった。
具体的に猿でもわかるように教えてー

5 :1:2001/06/26(火) 09:50
頼むーー

6 :名無しさん@1周年:2001/06/26(火) 10:28
マルチポストウザすぎ。
ここで相手にされなかったからって他に書くな。
こういう子は放置よ!!!

7 :名無しさん@1周年:2001/06/26(火) 12:53
>>6
氏ね

8 :6:2001/06/26(火) 13:09
>>7
はーい。

9 :あほあほ東大生:2001/07/01(日) 03:30
きみレポートは自分でやるものです

10 : ◆UQdKaCeg :02/09/23 16:29
昔やったのに…駄目だな自分。関係ないですが後半の}×5行は}}}}}ってやっちゃ…
駄目ですかそうですか。空間嫌い。

11 :山崎渉:03/01/11 05:30
(^^)

12 :名無しさん@3周年:03/02/23 19:30
http://www.omaemona.com/

13 :n145025.ap.plala.or.jp:03/02/25 17:32


14 :山崎渉:03/04/17 09:29
(^^)

15 :山崎渉:03/04/20 03:57
   ∧_∧
  (  ^^ )< ぬるぽ(^^)

16 :山崎渉:03/05/21 21:55
━―━―━―━―━―━―━―━―━[JR山崎駅(^^)]━―━―━―━―━―━―━―━―━―

17 :山崎渉:03/05/21 23:33
━―━―━―━―━―━―━―━―━[JR山崎駅(^^)]━―━―━―━―━―━―━―━―━―

18 :山崎渉:03/05/28 14:26
     ∧_∧
ピュ.ー (  ^^ ) <これからも僕を応援して下さいね(^^)。
  =〔~∪ ̄ ̄〕
  = ◎――◎                      山崎渉

19 :山崎 渉:03/07/15 12:56

 __∧_∧_
 |(  ^^ )| <寝るぽ(^^)
 |\⌒⌒⌒\
 \ |⌒⌒⌒~|         山崎渉
   ~ ̄ ̄ ̄ ̄

20 :山崎 渉:03/08/02 02:26
   ∧_∧
  (  ^^ )< ぬるぽ(^^)

21 :ぼるじょあ ◆yBEncckFOU :03/08/02 03:01
     ∧_∧  ∧_∧
ピュ.ー (  ・3・) (  ^^ ) <これからも僕たちを応援して下さいね(^^)。
  =〔~∪ ̄ ̄ ̄∪ ̄ ̄〕
  = ◎――――――◎                      山崎渉&ぼるじょあ

22 :山崎 渉:03/08/15 18:43
    (⌒V⌒)
   │ ^ ^ │<これからも僕を応援して下さいね(^^)。
  ⊂|    |つ
   (_)(_)                      山崎パン

23 :名無しさん@3周年:03/11/26 04:27
G-J よりは、GAUSSの消去法(LU分解)を使うべし。

5 KB
■ このスレッドは過去ログ倉庫に格納されています

★スマホ版★ 掲示板に戻る 全部 前100 次100 最新50

read.cgi ver 05.04.00 2017/10/04 Walang Kapalit ★
FOX ★ DSO(Dynamic Shared Object)