myjacobi.m

来自「基于matlab实现的实质分析实例的计算 对基本的插值的计算 可运行」· M 代码 · 共 21 行

function x=myjacobi(a,b,x0,n,e)%线性方程组的Jacobi迭代法

m=length(b);
x=x0;
for  k=1:n
    k
    for i=1:m
        
    x(i)=(b(i)-sum(a(i,:).*x0)+a(i,i)*x0(i))/a(i,i);
end
v=max(abs(x-x0))
    if max(abs(x-x0))<e
        break
    end
    x0=x

end
if k==n
    '迭代失败'
end