rt

• n为未知数个数
• 在eqs输入矩阵（等式左边）
• res输入对应的常数项（等式右边）
• 结果存储在res

for i in range(n-1):
    for j in range(i+1,n):
        for k in range(i+1,n):
            eqs[j][k]-=eqs[i][k]/eqs[i][i]*eqs[j][i]
        res[j]-=res[i]/eqs[i][i]*eqs[j][i]
        eqs[j][i]=0
for i in range(n-1,0,-1):
    res[i]/=eqs[i][i]
    for j in range(0,i):
        res[j]-=eqs[j][i]*res[i]
        eqs[j][i]=0

实测有效