#include<stdio.h>
#include<math.h>
int main()
{
int i,j, k;
double A[3][4] = {{2,8,2,14},{1,6,-1,13},{2,-1,2,5}};
double x1,x2,x3;
double m;
for(i=0; i<2; i++)
{
for(j=i+1; i<3;j++)
{
m = A[j][i]/A[i][i];
for(k=0; k<4;k++)
{
A[j][k] = A[j][k] - m*A[i][k];
}
}
x3 = A[2][3]/A[2][2];
x2 = (A[1][3]-A[1][2]*x3)/A[1][1];
x1 = (A[0][3]-A[0][1]*x2-A[0][2]*x3)/A[0][0];
}
printf("x1=%lf x2=%lf x3=%lf \n",x1,x2,x3);
return 0;
}
I2luY2x1ZGU8c3RkaW8uaD4KI2luY2x1ZGU8bWF0aC5oPgppbnQgbWFpbigpCnsKICAgIGludCBpLGosIGs7Cglkb3VibGUgQVszXVs0XSA9IHt7Miw4LDIsMTR9LHsxLDYsLTEsMTN9LHsyLC0xLDIsNX19OwoJCglkb3VibGUgeDEseDIseDM7Cglkb3VibGUgbTsKCiAgICBmb3IoaT0wOyBpPDI7IGkrKykKCXsKCQlmb3Ioaj1pKzE7IGk8MztqKyspCgkJewoJCQltID0gQVtqXVtpXS9BW2ldW2ldOwoJCQlmb3Ioaz0wOyBrPDQ7aysrKQoJCQl7CgkJCQlBW2pdW2tdID0gQVtqXVtrXSAtIG0qQVtpXVtrXTsKCQkJfQoJCX0KCQkKCQl4MyA9IEFbMl1bM10vQVsyXVsyXTsKCQl4MiA9IChBWzFdWzNdLUFbMV1bMl0qeDMpL0FbMV1bMV07CgkJeDEgPSAoQVswXVszXS1BWzBdWzFdKngyLUFbMF1bMl0qeDMpL0FbMF1bMF07Cgl9CglwcmludGYoIngxPSVsZgl4Mj0lbGYJeDM9JWxmIFxuIix4MSx4Mix4Myk7CiAgICByZXR1cm4gMDsKfQ==