資源の無駄使いlu_decompの返り値が(転置されたL, U)なのであとから転置する。 # coding: utf8 # 正方行列のLU分解 # 再帰で書いたLU分解 def lu_decomp(A): if len(A)==1: return [ [1] ], [ [A[0][0]] ] else: li = [row[0] / A[0][0] for row in A] ui = …