C语言高斯消元法与LU分解(3阶)实现

int main(){
   
	int k,i,j;
	float matrix[3][4] = {
   
		{
   -2,1.072,5.643,3},
		
		{
   -1,3.7212,4.623,2},
		{
   0.00001,2,3,1},
	};
	float  matrix1[3][4],matrix2[3][4];
	float n1,x3,x2,x1,n;
	for(i=0;i<3;i++)
		for(j=0;j<4;j++)
			matrix1[i][j]=matrix[i][j];
	//打印出原矩阵
	printf("原增广矩阵\n");
	for(i=0;i<3;i++){
   
		for(j=0;j<4;j++)
			printf("%lf\t",matrix1[i][j]);
		printf("\n");
	}

	printf("\n第一次消元\n");
	i=0;

		for(j=i;j<4;j++){
   
			for(k=i+1;k<3;k++){
   
				n1 = matrix[k][i]/matrix[i][i];
				matrix1[k][j] = matrix1[k][j] - matrix1[i][j]*n1;
			}
		}
		for(i=