雅可比迭代法MATLAB程序

function[y,n]=jacobi(A,b,x0,ep)ifnargin==3ep=1.0e-6;elseifnargin3errorreturnendendD=diag(diag(A));L=-tril(A,-1);U=-triu(A,10);B=D\(L+U);f=D\b;y=B*x0+f;n=1;whilenorm(y-x0)=epx0=y;y=B*x0+f;n=n+1;end