목차
없음
본문내용
i=0;i<2;i++){
j=i+1;
while(M[i][i]==0){
for(k=0;k<4;k++){
t = M[i][k];
M[i][k] = M[j][k];
M[j][k] = t;
}
j++;
}
for(j=i+1;j<3;j++){
t=M[j][i];
for(k=i;k<4;k++){
M[j][k]=M[j][k]-M[i][k]*t/M[i][i];
}
}
}
x3 = M[2][3]/M[2][2];
x2 = (M[1][3]-M[1][2]*x3)/M[1][1];
x1 = (M[0][3]-M[0][1]*x2-M[0][2]*x3)/M[0][0];
printf("x1 = %f\n",x1);
printf("x2 = %f\n",x2);
printf("x3 = %f\n",x3);
}
/*
Gauss Elimination with complete pivoting for Exercise 6.3.4.a
*/
#include
main(){
int i,j,k;
float M[3][4]={{1,-5,1,7},{10,0,20,6},{5,0,-1,4}};
float t,x1,x2,x3;
for(i=0;i<2;i++){
for(j=i+1;j<3;j++){
if(M[i][0]
for(k=0;k<4;k++){
t = M[i][k];
M[i][k] = M[j][k];
M[j][k] = t;
}
}
}
}
for(i=0;i<2;i++){
for(j=i+1;j<3;j++){
if(M[0][i]
for(k=0;k<3;k++){
t = M[k][i];
M[k][i] = M[k][j];
M[k][j] = t;
}
}
}
}
for(i=0;i<2;i++){
j=i+1;
while(M[i][i]==0){
for(k=0;k<4;k++){
t = M[i][k];
M[i][k] = M[j][k];
M[j][k] = t;
}
j++;
}
for(j=i+1;j<3;j++){
t=M[j][i];
for(k=i;k<4;k++){
M[j][k]=M[j][k]-M[i][k]*t/M[i][i];
}
}
}
x3 = M[2][3]/M[2][2];
x2 = (M[1][3]-M[1][2]*x3)/M[1][1];
x1 = (M[0][3]-M[0][1]*x2-M[0][2]*x3)/M[0][0];
printf("x1 = %f\n",x1);
printf("x2 = %f\n",x2);
printf("x3 = %f\n",x3);
}
j=i+1;
while(M[i][i]==0){
for(k=0;k<4;k++){
t = M[i][k];
M[i][k] = M[j][k];
M[j][k] = t;
}
j++;
}
for(j=i+1;j<3;j++){
t=M[j][i];
for(k=i;k<4;k++){
M[j][k]=M[j][k]-M[i][k]*t/M[i][i];
}
}
}
x3 = M[2][3]/M[2][2];
x2 = (M[1][3]-M[1][2]*x3)/M[1][1];
x1 = (M[0][3]-M[0][1]*x2-M[0][2]*x3)/M[0][0];
printf("x1 = %f\n",x1);
printf("x2 = %f\n",x2);
printf("x3 = %f\n",x3);
}
/*
Gauss Elimination with complete pivoting for Exercise 6.3.4.a
*/
#include
main(){
int i,j,k;
float M[3][4]={{1,-5,1,7},{10,0,20,6},{5,0,-1,4}};
float t,x1,x2,x3;
for(i=0;i<2;i++){
for(j=i+1;j<3;j++){
if(M[i][0]
t = M[i][k];
M[i][k] = M[j][k];
M[j][k] = t;
}
}
}
}
for(i=0;i<2;i++){
for(j=i+1;j<3;j++){
if(M[0][i]
t = M[k][i];
M[k][i] = M[k][j];
M[k][j] = t;
}
}
}
}
for(i=0;i<2;i++){
j=i+1;
while(M[i][i]==0){
for(k=0;k<4;k++){
t = M[i][k];
M[i][k] = M[j][k];
M[j][k] = t;
}
j++;
}
for(j=i+1;j<3;j++){
t=M[j][i];
for(k=i;k<4;k++){
M[j][k]=M[j][k]-M[i][k]*t/M[i][i];
}
}
}
x3 = M[2][3]/M[2][2];
x2 = (M[1][3]-M[1][2]*x3)/M[1][1];
x1 = (M[0][3]-M[0][1]*x2-M[0][2]*x3)/M[0][0];
printf("x1 = %f\n",x1);
printf("x2 = %f\n",x2);
printf("x3 = %f\n",x3);
}
소개글