解线性方程组的迭代方法(Jacobi，Gauss-Seidel)

function X=my(A,b)
%A为系数矩阵.
%b为右端常向量.
%X0为初始层，默认取0.
%mytol表示允许误差, 要求默认值为1e-6.
%nargin表示输入参数个数
T=diag(A);
D=diag(T);
L=-tril(A,-1);
U=-triu(A,1);
X1=ones(size(b));
n=1;  
%T是A主对角线元素，D是T的对角矩阵，X1是元素全1方阵
if nargin==2   %如果参数=2
    X0=zeros(size(b)); %X0是0矩阵
    mytol=1e-6; 
else
    error('输入参数数目不对')
end
[r,c]=size(A);
if r==c
    if abs(det(A))>=mytol  %det计算矩阵A的行列式的值，行列式不为0可逆
        if abs(det(D))>=mytol  %对角矩阵不为0，说明主对角线无0
            B=D^(-1)*(L+U);
            d=max(abs(eig(B))); %eig特征值绝对值的最大值
            if d<1  %谱半径<1,迭代法收敛
                disp('迭代法收敛')
                while max(abs(X0-X1))>=mytol  
                %终止循环,后一层和前一层结果的差的绝对值的最大值小于某个给定的小量mytol
                    X1=X0; %把X0赋给X1,相当于X1是前一层
                    X0=B*X0+D^(-1)*b;  %用前一层的计算得出后一层
                    n=n+1;
                end