一、 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。

幂法求矩阵的特征值和特征向量 算法9.1 1． 输入矩阵A，初始向量X，误差限ε最大迭代次数N。 2． 置k=1，u=0。 3． 求xr=>λ，|xr|=max|xr|,i=1,….,n。 4． 计算 Y=X/λ，X=AY 5． 若|λ-u|〈ε，输出λ，X，停机；否则，转步骤6。 6. 若k<N，k+1=>k，λ=>u,转步骤3；否则，输出错误信息，停机

高斯赛德尔迭代实现，其中矩阵运算均是自写，纯原创-Gauss seidel iteration implementation, including matrix operations are from, pure original