Jacobiff雅可比

function x=Jacobiff(a,b,eps) n=length(a); d=diag(diag(a)); l=tril(a,-1); u=tril(a,1); x0=diag(zeros(n)); k=1 x=-inv(d)*(l+u)*x0+inv(d)*b while norm(x-x0)>eps k=k+1 p=x x=-inv(d)*(l+u)*x0+inv(d)*b x0=p end