【文件属性】:
文件名称: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