学习数值分析这门课程时,对高斯消元法及主元素法部分代码的编写,希望能帮到大家,也希望能有大神对我做的不到位的地方进行指正。直接上代码。
def guess(A,B): #A-系数矩阵,B-对应方程的矩阵
n = len(A)
if n != len(A[0]): #确定系数矩阵符合未知数等于方程数
return False
for k in range(n-1):
# for l in range(k,n):
# if abs(A[l][k]) > abs(A[k][k]):
# A[k],A[l] = A[l],A[k]
# B[k],B[l] = B[l],B[k]
# 主元素法加上这几行进行排序操作
if A[k][k] == 0: #确定Akk不为零
ret