|
r(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
|
|
메뉴펼치기
