ガウス消去

質問

行列3*3で４列目に答えが入る、ガウス消去の前進消去のプログラムを考えたのですが、３行目まで正規化してしまいます。
前進消去したのを行列mに格納しています。

int main()
{
int i, j, k,n=3;
double a[3][4]={{ 3, 2,-3,-2},
{ 2,-2, 1, 7},
{-3, 6, 5,28}};
double m[3][4];

/* 前進消去 */
for( i = 0; i < n; i++ )
{
for( j = i+1; j <= n; j++ )
{
m[i][j] = a[i][j]/a[i][i];
m[i][i]=1;
for( k = i+1; k <= n+1; k++ )
{

a[k][j] = a[k][j] -a[k][i] * m[i][j];

m[j][i]=0;
}}}

回答No.1

BLUEPIXY · ベストアンサー

余計なアドバイスかもしれないけど、ｍ は、(1,0 になるとわかっている部分で構成されるから）別に要らないんじゃないかと思う
----------------------------------------------------------------
#include <stdio.h>

int main(){
int i, j, k, n=3;
double a[3][4]={
{ 3, 2,-3,-2},
{ 2,-2, 1, 7},
{-3, 6, 5,28},
};
//初期状態
for(i=0;i<n;i++){
for(j=0;j<=n;j++)
printf("%g ", a[i][j]);
printf("\n");
}
printf("\n");

for( i = 0; i < n; i++ ){
for( j = i+1; j <= n; j++ ){
a[i][j] /= a[i][i];
}
for( k = 0; k < n; k++ ){
if(k != i){
for(j = i+1;j <= n;j++){
a[k][j] -= a[k][i] * a[i][j];
}
}
}
}
//点検
for(i=0;i<n;i++){
printf("%lf\n", a[i][n]);
}

return 0;
}