【数值计算方法】线性方程组的迭代解法-数值实验


from formu_lib import *
import numpy as np

A=np.array([[-55,-5,12],        
            [21,36,-13],
            [24,7,47]])
b=np.array([41,52,12])

w=lambda t:0.1*t

xs,ys,ts=[],[],[]
for i in range(1,20):
    _,err=SORIter(A,b,w(i))
    xs.append(list(range(len(err+1))))
    ys.append(err)
    ts.append(f"SOR iter error with w={w(i)}")
    print(f"i={i}")