gauss_method.m

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

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