From 3d1d16df101d2eb221d21c630468cac3c8d05f0b Mon Sep 17 00:00:00 2001 From: 123 <629825095@qq.com> Date: Wed, 18 Jun 2025 23:22:51 +0800 Subject: [PATCH] eee --- 227-4.py | 20 ++++++++++++++++---- 1 file changed, 16 insertions(+), 4 deletions(-) diff --git a/227-4.py b/227-4.py index ad66601..68b01af 100644 --- a/227-4.py +++ b/227-4.py @@ -22,12 +22,24 @@ def SqrtSolve(A,b): for k in range(j): L[getIndexFromDownMatrix(i,j)] -= L[getIndexFromDownMatrix(i,k)]*L[getIndexFromDownMatrix(j,k)] L[getIndexFromDownMatrix(i,j)] /= L[getIndexFromDownMatrix(j,j)] - print(L) + # 打印下三角矩阵 + print("下三角矩阵 L:") + for i in range(n): + L_row = [] + for j in range(n): + if j <= i: + L_row.append(L[getIndexFromDownMatrix(i,j)]) + else: + L_row.append(0) + print(L_row) + # print(L) for i in range(n): for k in range(i): b[i] -= L[getIndexFromDownMatrix(i,k)]*b[k] b[i] /= L[getIndexFromDownMatrix(i,i)] - + # 打印 b 向量 + print("y 向量:") + print(b) for i in range(n-1,-1,-1): for k in range(i+1,n): b[i] -= L[getIndexFromDownMatrix(k,i)]*b[k] @@ -42,5 +54,5 @@ if __name__ == "__main__": b = [10,5,4] - print("x:") - print(SqrtSolve(A,b)) + # print("x:") + print("x: \n",SqrtSolve(A,b))