根据Doolittle LU分解，矩阵通常可以分解为A = LU。当A是非单数时，就存在这种分解。
对于以下示例，A = {{1,2,1}，{2,4,4}，{1,4,9}}。Example 它的行列式是-4（非零），因此存在矩阵L和U.
但是，当我按照算法对A进行分解时， pseudo code 得到u11 = 1，u12 = 2，u13 = 1，l21 = 2，l31 = 1，u22 = 0，u23 = 2，
然后，它在此步骤未能获得l32。 4 = l31 * u12 + l32 * u22
另外，如果u22 = 0，则表示det（U）= 0 =＆gt; DET（A）= DET（LU）= 0。它与det（A）= - 4 !!!
相矛盾任何人都可以指出我做错了什么？
