一、　　Seidel迭代法

算法8.2：

1．　　输入矩阵A，右端项b ，初始迭代向量 X(0) ，维数n，容许误差ε ，最大容许迭代次数 N。

2．　　置k=1.

3．　　计算

X1=(b1-Σa1jxj(0))/a11 j=2,3,…….,n

Xi=(bi-Σaijxj-Σaijxj(0))/aii, i=2,3,……,n-1

Xn=(bn-Σanjxj)/ann j=1,2,…….,n-1

4．若||X-X(0)|| <ε，输出X , 停机；否则，转步骤5。

5．若k<N，k+1=>k,X=>X(0),转步骤3；否则输出失败信息，停机