Jacobi Iterative Method

function apX=JacobiIterative(n,A,b,XO,TOL,N)
k=1;
x=zeros(n,1)

while k<=N
    for i=1:n
        x(i)=((-1)*(A(i,:)*XO)+b(i)+A(i,i)*x(i))/A(i,i)
    end
    k=k+1;
    for i=1:n
        XO(i)=x(i)
   end
end
XO
k

TOL部分未補，要先給上一個A的矩陣 n是矩陣大小 b是1Xn的向量 要找 Ax=b 的 solution.