利用python实现高斯消元法（非列主元） 记录计算物理课程的课后作业

记录下本学期计算物理课程的课后作业。第一道题为高斯消去法

mt = [[2,-1,3,1],[4,2,5,4],[1,2,0,7]]

m = len(mt)
n = len(mt[0])

for i in range(0,m):
    for k in range(0,i):
        l = mt[i][k]/mt[k][k]
        for j in range(0,n):
            mt[i][j] = mt[i][j] - mt[k][j] * l


res = [0] * m
res[m-1] = mt[m-1][n-1]/mt[m-1][n-2]  #计算最后一个x的值

for i in range(m-2,-1,-1):       #回代nxl
    temp = 0
    for j in range(i+1,n-1):
        c = mt[i][j]*res[j]
        temp += c
    res[i] = (mt[i][n-1] - temp)/mt[i][i]

print(res)