有两种方法：
待定系数。直接设L,U的元素，然后计算L*U=A，解出L和U。虽然这种办法数值计算量大些，但是过程简单易理解。用在编程里更好
左乘行初等矩阵（初等行变化），一步步乘Pi，把A的对角线下面元素消去，然后剩下的就是U。Pn*......P2*P1*A=U，令P=Pn*P(n-1)*...P1，则有P*A=U，所以A=P^(-1)*U。这里P^(-1)是指P的逆。
建议用待定系数法，计算量不大的话。
