LU分解

由以上知识可以推导出：

程序代码：

#include <stdio.h>
int main(){
	double a[5][5]={{0,0,0,0,0},{0,4,2,1,5},{0,8,7,2,10},{0,4,8,3,6},{0,6,8,4,9}};
	double l[5][5]={0},u[5][5]={0};
	//输出a
	printf("a:\n"); 
	for(int i=1;i<5;i++){
		for(int j=1;j<5;j++){
			printf("%3.2f    ",a[i][j]);
		}
		printf("\n");
	}
	 
	for(int k=1;k<5;k++){
		//得到l的第k+1列 
		for(int i=k+1;i<5;i++){
			a[i][k] = a[i][k]/a[k][k];
		}
		//得到u的第k+1行 
		for(int i=k+1;i<5;i++)
		{
			for(int j=k+1;j<5;j++){
				a[i][j] = a[i][j] - a[i][k]*a[k][j];
			}
		}
	}
	//将改造好的a的值赋给l和u 
	for(int i=1;i<5;i++){
		for(int j=1;j<5;j++){
			if(j<i) l[i][j]=a[i][j];
			else	u[i][j]=a[i][j];
		}
		l[i][i]=1;
	}
	printf("l:\n");
	//输出l
	for(int i=1;i<5;i++){
		for(int j=1;j<5;j++){
			printf("%3.2f    ",l[i][j]);
		}
		printf("\n");
	}
	printf("u:\n");
	//输出u 
	for(int i=1;i<5;i++){
		for(int j=1;j<5;j++){
			printf("%3.2f    ",u[i][j]);
		}
		printf("\n");
	}
	return 0;
}

运行结果：