顺序高斯消去法

Code

#普式化时将3作为n输入即可

arr = [[] for i in range(3)]

for i in range(3) :
    for j in range(4) :
        v = input()
        v = float(v)
        arr[i].append(v)
        
for i in range(3 - 1) :
    for j in range(i + 1, 3) :
        rate = arr[j][i] / arr[i][i]
        for k in range(0, i + 1) :
            arr[j][k] = 0
        for k in range(i + 1, 4) :
            arr[j][k] -= rate * arr[i][k]

#这里定义为10，通式时确保list足够大
x = [0] * 10

for i in range(3) :
    index = 3 - i
    sum = arr[index - 1][3]
    for j in range(i) :
        sum -= arr[index - 1][3 - 1 - j] * x[3 - j]
    x[index] = sum / arr[index - 1][index - 1]


for i in range(1, 3 + 1) :
    print(x[i])