ガウスの消去法のプログラムを教えて下さい。

質問

3x+y+z=8
6x+2y+3z=4
9x+4y+5z=12
を計算する部分ピボット選択付きのガウスの消去法と、選択のないガウスの消去法のプログラムを教えて下さい＞＜。

何度も作り直しているのですがなかなか上手くいきません。よろしくお願いします。

回答No.1

DIooggooID

どのような処理の流れにしていますか？

　・繰り返し処理のところが、キーポイントになる感じがします。

補足

#include<stdio.h>
#include<stdlib.h>
#include<math.h>

int main(){
double a[3][3]={
{3,1,1},
{6,2,3},
{9,4,5}
};
double b[3]={8,4,12};
double max[3];
double maxp;
double temp;
double x[3];
double m;
int pivot=0;
int i,j,k;

for(k=0;k<3;k++){

for(j=k;j<3;j++){
if(maxp<fabs(a[j][k])){
pivot=j;
maxp=fabs(a[j][k]);
}
}

for(i=k;i<3;i++){
temp=a[k][i];
a[k][i]=a[pivot][i];
a[pivot][i]=temp;
}

temp=b[k];
b[k]=b[pivot];
b[pivot]=temp;

m=a[k][k];
for(i=k;i<3;i++){
a[k][i]=a[k][i]/m;
}
b[k]=b[k]/m;

for(j=k+1;j<3;j++){
m=a[i][k]/a[k][k];
for(i=k;i<3;i++){
a[j][i]-=a[k][i]*m;
}
b[j]-=b[k]*m;
}

for(j=2;j>=0;j--){
x[j]=b[i];
for(i=j+1;i<3;i++){
x[j]-=a[j][i]*x[i];
}
}

for(j=0;j<3;j++){
for(i=0;i<3;i++){
printf("%f",a[j][i]);
}
printf("\n");
printf("%f",x[j]);
}

return(0);
}
}

というプログラムを作ったのですが上手くいきません＞＜。