jacobi.m

来自「用Jacobi叠待法解线性方程组 function Jacobi(A,b,n,」· M 代码 · 共 26 行

function Jacobi(A,b,n,x0,e,N)
k=1
x=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
                x(i)=x(i)-A(i,j)*x0(j)/A(i,i);
            end
        end
        fprintf('%f\t',x(i));
    end
    if norm(x-x0)<e
            break;
    else
        k=k+1;
        x0=x;
    end
    fprintf('\n');
end
 fprintf('\n***********end*************\n',k);
x0
x