mygauss_seidel.m

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

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

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

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