Fixed Bug #2257

Author: Svn-Sp

arithmetic_analysis/lu_decomposition.py

@@ -13,15 +13,15 @@ def LUDecompose(table):
     if rows != columns:
         return []
     for i in range(columns):
-        for j in range(i - 1):
+        for j in range(i):
             sum = 0
-            for k in range(j - 1):
+            for k in range(j):
                 sum += L[i][k] * U[k][j]
             L[i][j] = (table[i][j] - sum) / U[j][j]
         L[i][i] = 1
-        for j in range(i - 1, columns):
+        for j in range(i, columns):
             sum1 = 0
-            for k in range(i - 1):
+            for k in range(i):
                 sum1 += L[i][k] * U[k][j]
             U[i][j] = table[i][j] - sum1
     return L, U