#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=i; 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;
}
I2luY2x1ZGU8c3RkaW8uaD4KI2luY2x1ZGU8bWF0aC5oPgppbnQgbWFpbigpCnsKICAgIGludCBpLGosIGs7Cglkb3VibGUgQVszXVs0XSA9IHt7Miw4LDIsMTR9LHsxLDYsLTEsMTN9LHsyLC0xLDIsNX19OwoJZG91YmxlIHgxLHgyLHgzOwoJZG91YmxlIG07CgogICAgZm9yKGk9MDsgaTwyOyBpKyspCgl7CgkJZm9yKGo9aSsxOyBpPDM7aisrKQoJCXsKCQkJbSA9IEFbal1baV0vQVtpXVtpXTsKCQkJZm9yKGs9aTsgazw0O2srKykKCQkJewoJCQkJQVtqXVtrXSA9IEFbal1ba10gLSBtKkFbaV1ba107CgkJCX0KCQl9Cgl9CgoJeDMgPSBBWzJdWzNdL0FbMl1bMl07CgkJeDIgPSAoQVsxXVszXS1BWzFdWzJdKngzKS9BWzFdWzFdOwoJCXgxID0gKEFbMF1bM10tQVswXVsxXSp4Mi1BWzBdWzJdKngzKS9BWzBdWzBdOwoKCXByaW50ZigieDE9JWxmCXgyPSVsZgl4Mz0lbGYgXG4iLHgxLHgyLHgzKTsKICAgIHJldHVybiAwOwp9