jacobi_method.m

来自「程序包中含有Lagrange插值、Newton插值、Hermite插值、jaco」· M 代码 · 共 30 行

function [x,m]=jacobi_method(A,b,kmax)
n=length(b);
x=zeros(n,1);
xk=zeros(n,1);
e=0.00001;
for k=1:kmax
    xk=x;
    for i=1:n
        sum=0;
        for j=1:n
            if j~=i
                sum=sum+A(i,j)*xk(j);
            end            
        end        
        x(i)=(b(i)-sum)/A(i,i);        
    end
    if norm(x-xk,1)<e
        m=k-1;
        break;
    end
end
if norm(x-xk,1)>e
    disp('iteraive times reach upper limit');
end
end