编辑代码

#include <stdio.h>
int main () {
    double a[4][4],b[4],x1,x2,x3,x4;
    int i,j;
    a[0][0]=1.1161;a[0][1]=0.1254;a[0][2]=0.1397;a[0][3]=0.1490;
    for(i=1;i<4;i++)
    {
        double t=-a[i][j-1]/a[i-1][j-1];
        b[i]+=b[i-1]*t;
        for(j=i;j<4;j++)
            a[i][j]+=a[i-1][j]*t;
    }
    x4=b[3]/a[3][3];
    x3=(b[2]-a[2][3]*x4)/a[2][2];
    x2=(b[1]-a[1][3]*x4-a[1][2]*x3)/a[1][1];
    x1=(b[0]-a[0][3]*x4-a[0][2]*x3-a[0][1]*x2)/a[0][0];
    printf("%lf %lf %lf %lf",x1,x2,x3,x4);
	return 0;
}