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