Gauss-Jordanの掃出法
■ このスレッドは過去ログ倉庫に格納されています
0001デフォルトの名無しさん
NGNGで、下の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();
}
0002デフォルトの名無しさん
NGNG0003デフォルトの名無しさん
NGNG0004デフォルトの名無しさん
NGNG教えれーー
0005デフォルトの名無しさん
NGNGhttp://mentai.2ch.net/test/read.cgi?bbs=sim&key=993456114
0006デフォルトの名無しさん
NGNG0007デフォルトの名無しさん
NGNG0008デフォルトの名無しさん
NGNGおまえが教えてやれよ。
0009デフォルトの名無しさん
NGNGごめん。
みんなバカだから無理。
00101
NGNGしかたないですね・・・
馬鹿しかいないのですね、ここには
0011マジレス君
NGNG0012デフォルトの名無しさん
NGNGそうでーす。君みたいなレベルの高い人にはここは似合わないよ。
■ このスレッドは過去ログ倉庫に格納されています