연립 1차 방정식의 풀이를 위한 가우스 소거법(C++)

본문내용

[j][i])
{
pivot = j;
con = MatrixA[j][i];
}
}
for(j=0;j {
Sort[j] = MatrixA[pivot][j];
MatrixA[pivot][j] = MatrixA[i][j];
MatrixA[i][j]=Sort[j];
}
con = 0.0;
}
/////////////////// Gauss Elimination - forward elimination /////////////////////////////////
int row_in = row-1; //연산 상에서 배열의 번호 ; 0,1....
int column_in = column -1;// //연산 상에서 배열의 번호 ; 0,1....
for(k=0;k {
for(i=k+1;i {
factor = MatrixA[i][k]/MatrixA[k][k];
for(j=k+1;j {
MatrixA[i][j] = MatrixA[i][j]-factor*MatrixA[k][j];
}
MatrixA[i][column_in]=MatrixA[i][column_in]-factor*MatrixA[k][column_in];
}
}
Solution[row_in] = MatrixA[row_in][column_in]/MatrixA[row_in][row_in];
////////////////// Gauss Elimination - Back substitution //////////////////////////////////
for(i=row_in;i>0;i--)
{
sum = MatrixA[i-1][column_in];
for(j=i+1;j {
sum = sum - MatrixA[i-1][j-1]*Solution[j-1];
}
Solution[i-1] = sum/MatrixA[i-1][i-1];
}
cout<<"The Answer is\n";
for(i=0;i {
cout<<"X"< }
cout<<"\n";
////////////// 이차원 배열 동적 할당 해제////////////////////////////
for(k=0;k delete[] MatrixA[k];
}
delete[] MatrixA;
delete[] Sort;
delete[] Solution;
}

키워드

가우스,   가우스 소거법,   gauss,   C++,   동적할당,   연립방정식