C语言实现方程组LU分解法和列主元LU分解法

根据老师所讲LU分解法和列主元LU分解法本质就是高斯消去法和列主元高斯消去法，

将L部分变为都变为0，则就化为高斯消去法，我们在方程求回代的过程中只用到了U和y，不需要用到L，那个L部分实际保存了每一次运算的倍数关系，详细可自己举数学例子验证便可知。

以例子解释规律，我们以（1，1）表示第一行第一列数,则（3，1）表示第三行第一列，那么第一行所有数保持不变，第一列数全部除上第一行第一列。 然后原来（2，2）-(2,1)*(1,2)即：5-0.7*7。 6-0.7*8，5-0.7*7而第二列为（6-0.8*7）/0.1,    (5-0.7*7)/0.1.先进性行运算，再进行列运算，第（3，3）则为2=（10-0.8*8-4*0.4），行则:3=(9-0.8*7-4*0.1),-6=(6-0.8*10-4*1）列则：1.5=（9-0.7*8-1*0.4）/2

大家应该能总结出规律了：所有数字更新完的上三角为U，下三角为L，L的对角线都为1.

代码实现：

#include<stdio.h>
# define M 100
int main()
{
	int