LU分解法求解线性方程组

LU分解法是求解线性方程组的一种算法。

先将系数矩阵分别转换成上三角矩阵U和下三角矩阵L，其中U(k,j)=a(k,j) - ∑L(k,m)*U(m,j)    m

求解Ly=b，b为因变量矩阵；

求解Ux=y；

c++代码如下：

#include

using namespace std;

#define  n  3    //矩阵阶数
int main()
{
	int i, j, k, m;
	float t, a[n][n], u[n][n], l[n][n], b[n], x[n], y[n];
	for (i = 0; i= 0; i--)
	{
		t = 0;
		for (k = i + 1; k