首先从最简单的情况考虑，这也是我们做研究、做学术、做工程必须要时刻牢记心中的一点，很多人喜欢一上来就把所有问题、把最复杂的情况、把方方面面都给考虑到，然后再开始实现他的想法，我自己也有这个习惯，但是，这并不是一个好习惯，一上来就好高骛远、就想着高大上这本质上是一种急功近利的表现，那样的话你会陷入到各种各样的技术细节当中，你会想半天却仍然写不出半点实质性的东西出来，所以最好的办法是，先考虑最简单、最核心的情况，这样不仅大大降低问题的复杂度，同时也为将来进一步扩展程序、解决更复杂的情况打下了一个坚实的基础
1.Lu分解
U1j=a1j (j=1,2,n) li1=ai1/u11(i=2,3,,n) U的第一行和L的第一列元素。
矩阵A*B不等于矩阵B*A。
matlab里的LU分解为[L,U]=lu(B)
与科学计算与工程那本书得到的结果互为转置。
回答：如果U为单位上三角矩阵时，成为Crout分解，如果L为单位下三角矩阵时称为Doolittle分解，即L的对角线为1。
2. 求解
这个按照doolittle分解，
matlab的结果 [L,U,p]=lu(B)