文件名称：Jacobiff雅可比
文件大小：20KB
文件格式：DOC
更新时间：2013-07-18 13:20:21
雅可比
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