python——赛得尔迭代法求解方程组

方程组如下：

我们使用赛得尔迭代法求解改方程组，运行的误差在0.00001。代码如下：

---

import numpy as np
ef g(x,x0,n,d):#赛德尔迭代法
    x1=np.ones(n);y=np.ones(n);y[:-1]=-1;qd=0;x2=np.ones(n);x3=np.ones(n);x3[:]=x0[:]
    while np.any(abs(x1-x3)>d):
        qd+=1
        for u in range(n-1):
            x1[u]=x3[u]
        for j in range(n-1):
            tq=0
            for m in range(n):
                if m!=j:
                    tq+=(x[j,:]/x[j,j]*y)[m]*x3[m]
            x3[j]=tq;x2[j]=tq
        for re in range(n-1):
            x3[re]=x2[re]
    return x3,qd

---

运行结果：

---

print('赛德尔迭代法,最终的解X_i(角标从小到大):{},及迭代次数K:{}'.format(g(x,x0,n,d)[0][:-1],g(x,x0,n,d)[1]))



赛德尔迭代法,最终的解X_i(角标从小到大):[1.0000023  1.99999984 2.99998957 0.99998176],及迭代次数K:28

---