%jacobi.m
function s=jacobi(A,b,X0,epx)
%A為系數(shù)矩陣，b即為Ax=b中的b,X0為迭代初值，eps為控制精度
D=diag(diag(A));  %求出對角陣
D=inv(D);  %求出對較陣的逆矩陣
L=tril(A,-1);  %求出嚴(yán)格下三角矩陣
U=triu(A,1);  %示出嚴(yán)格上三角矩陣
B=-D*(L+U);
f=D*b;
s=B*X0+f;
while norm(s-X0)>=eps
    X0=s;
    s=B*X0+f;  %迭代的矩陣表達形式
end
return