整理文档很辛苦，赏杯茶钱您下走！
还剩 ... 页未读，继续阅读 >>
免费阅读已结束，点击下载阅读编辑剩下 ... 页
阅读已结束，您可以下载文档离线阅读编辑
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
