第27号专栏

function Gauss_Seidel(A,b,n,x0,E,N)k=1x=zeros(n,1);构造列向量 fprintf(/n**********start**************/n,k);for k=1:N    fprintf(%d/t,k);    %calculate x(1)    x(1)=b(1)/A(1,1);    for j=2:n        x(1

function Jacobi(A,b,n,x0,E,N)k=1x=zeros(n,1); fprintf(/n**********start**************/n,k);for k=1:N    fprintf(%d/t,k);    for i=1:n        x(i)=b(i)/A(i,i);        for j=1:n            if j~=i