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